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ABSTRACT 


Investigation  of  the  Effect  of  Two-Dimensional  Cavities  on  Boundary  Layers  in 

an  Adverse  Pressure  Gradient 
by 

Richard  J.  Margason 


The  present  investigation  evaluated  one  aspect  of  the  feasibility  of  the  use 
of  multiple  cavities  as  an  airfoil  high-lift  device.  The  effects  of  cavities  on  the 
boundary  layer  characteristics  in  several  pressure  gradients  were  determined 
experimentally  and  computationally.  Experimentally,  it  was  found  that  up  to  four 
cavities  could  be  deployed  with  only  a  small  change  to  the  boundary  layer 
profiles  downstream  of  the  cavities  and  without  significantly  modifying  the 
resultant  streamwise  pressure  distribution.  From  the  computational  results  for 
both  of  the  wind  tunnel  test  section  lengths  used  in  the  experimental 
investigation,  it  was  found  that  a  grid  which  provided  a  converged  solution  in 
less  than  a  few  hundred  iterations  was  needed  before  a  reasonable  comparison 
with  experimental  data  could  obtained.  It  was  also  found  for  these  converged 
solutions  that  the  appropriate  grid  clustering  and  density  as  well  as  the  cell  size 
required  for  a  satisfactory  solution  was  not  always  apparent  before  comparing 
computational  results  with  experimental  data.  Overall,  the  investigation  results 
show  that  a  multiple  cavity  high-lift  concept  may  be  feasible. 
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I.  INTRODUCTION 


Aircraft  wings  are  usually  sized  by  their  cruise  requirements  which  occur  at 
a  relatively  low  lift  coefficient.  The  increased  lift  coefficient  needed  for  low- 
speed  flight,  including  take-off  and  landing,  usually  requires  the  deployment  of 
high-lift  devices.  These  devices  increase  either  wing  area  and/or  lift  coefficient  at 
a  given  angle-of-attack.  Several  examples  of  high-lift  devices  are  shown  in  Figure 
1.  Most  modem  aircraft  use  slats  and  flaps  (Figure  1(a))  which  extend  to  increase 
wing  area  and  deflect  to  increase  camber  (i.e.,  lift  coefficient).  Additionally,  less 
conventional  high-lift  concepts  include  trailing-edge  tabs  (Figure  1(b)),  variable- 
camber  airfoils  (Figure  1(c)),  and  spanwise  blowing.  The  pressure  coefficient  on 
airfoils  typically  has  a  favorable  gradient  from  flow  stagnation  location  near  the 
leading  edge  to  the  minimum  negative  pressure  coefficient  peak  which  is  usually 
located  several  percent  of  the  chord  length  from  the  leading  edge. 

For  high-lift  conditions,  the  pressure  coefficient  rises  from  the  negative 
pressure  peak  back  to  the  ambient  static  pressure  near  the  trailing  edge.  This  is  the 
adverse  pressure-gradient  region  (dp/dx  >  0).  If  the  pressure  gradient  is  too  large 
the  lift  is  reduced  because  the  flow  will  separate  on  the  wing  upper  surface  and  the 
suction  pressure  peak  is  reduced  and  moves  toward  the  leading  edge.  A  cavity  or 
multiple  cavities  could  be  useful  in  delaying  separation  to  a  more  adverse  pressure 
gradient.  The  usefulness  of  multiple  cavities  for  maintaining  attached  flow  in  an 
adverse  pressure  gradient  is  determined  by  their  effect  on  the  viscous  flow  in  the 
boundary  layer  downstream  of  the  cavities.  Maintenance  of  attached  flow  requires 
that  the  flow  over  the  cavities  and  downstream  of  them  retain  enough  momentum 
to  overcome  the  loss  due  to  an  adverse  pressure  gradient,  shear  layer  flow 
gradients,  and  viscous  dissipation.  The  effects  are  greater  near  the  wall  because  of 
the  reduced  velocity  near  the  surface.  If  the  velocity  at  the  wall  is  reduced  to  zero 
the  flow  separates  and  produces  flow  forward,  opposite  the  freestream  direction. 
As  a  result,  boundary  layer  assumptions  are  no  longer  valid. 

High-lift  aerodynamics  has  been  the  subject  of  numerous  investigations  as 
described  in  a  classic  review  by  A.  M.  O.  Smith  [Ref.  1]  about  twenty  years  ago. 
This  review  described  conventional  high-lift  devices  and  included  a  list  of  the  ten 
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most  important  basic  theoretical  problems  of  high-lift  aerodynamics.  The  list 
emphasized  development  of  computational  methods  capable  of  representing  three- 
dimensional  flow  with  boundary-layer  separation  and  merging  multi-element 
airfoil  boundary  layers.  There  have  been  many  experimental  investigations  and 
applications  of  computational  methods  to  high-lift  aerodynamics  during  the  past 
twenty  years.  As  a  consequence,  there  has  been  significant  progress  toward  the 
solution  of  the  basic  research  needs  listed  by  Smith.  Additionally,  the  review 
indicated  the  need  for  new  "inventions"  such  as  spanwise  blowing  and  trapped- 
vortex  concepts  to  provide  improved  lifting-surface  flow  control. 

There  have  been  several  attempts  to  develop  the  trapped-vortex  concept. 
Even  Leonardo  da  Vinci  observed  and  sketched  very  realistic  recirculating  eddies 
due  to  a  bluff  obstacle  (or  spoiler)  in  the  fifteenth  century.  Recently  it  was 
suggested  [Ref.  2]  that  several  spoilers  located  at  different  chordwise  locations 
could  be  deployed  to  form  several  cavities  in  the  chordwise  direction  (Figure  1(d). 
This  may  provide  a  means  of  increasing  upper  surface  camber  and  thereby 
increase  wing  lift.  The  spoilers  could  be  deployed  as  a  high-lift  device  for  steady 
flow  conditions.  In  other  scenarios,  the  spoilers  could  be  deployed  either  for 
transient  maneuvers  of  a  combat  aircraft  at  high  angles-of- attack  or  for  increasing 
the  lift  of  a  retreating  rotorcraft  blade.  In  still  another  application,  the  sudden 
opening  of  a  cavity  near  the  leading  edge  (Figure  1(e))  of  a  retreating  rotorcraft 
blade  may  suppress  or  delay  dynamic  stall.  Cavity  flows  may  also  be  caused  by 
finned  surfaces,  windows,  bomb  bays,  landing  gear  bays,  finned  heat  transfer 
surfaces,  and  other  surface  imperfections.  Some  cavity  flow  research  has  been 
conducted  specifically  to  study  the  more  general  problem  of  flow  separation. 

Detailed  cavity  experiments  were  first  reported  in  the  early  1950's.  Since 
then  a  large  number  of  investigations  of  cavity  flow  have  been  conducted.  This 
research  has  most  often  concentrated  on  a  single,  two-dimensional  rectangular 
cavity  in  a  uniform  freestream  flow  with  no  streamwise  static-pressure  gradient.  A 
survey  of  the  literature  was  made  to  locate  investigations  of  either  the  effect  of 
cavities  on  boundary-layer  characteristics,  especially  for  high-lift  conditions;  i.e., 
in  an  adverse  pressure-gradient.  No  experimental  data  were  found  in  the  literature 
search  which  documented  the  effect  of  either  adverse  pressure  gradients  or 
multiple  cavities  on  the  boundary-layer  development  or  on  the  pressure 
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distributions  in  or  near  the  cavities.  No  data  for  either  single  or  multiple  cavities 
were  found  which  documented  their  effect  on  the  local  flow  in  the  vicinity  of  the 
cavities.  These  effects  include  the  influence  of  a  door  or  spoiler  which  closes  the 
cavity  opening. 

A  basic  understanding  is  needed  to  determine  how  multiple  cavities  could 
be  used  on  an  airfoil  to  either  delay  flow  separation  or  to  increase  its  camber  and 
obtain  an  increased  lift  at  a  given  angle-of-attack  with  suitable  drag  characteristics. 
The  present  investigation  simulated  a  lifting  airfoil  using  a  wind  tunnel  with  test 
section  liners  to  represent  different  pressure-gradients  which  could  be  associated 
with  an  airfoil  at  selected  angles-of-attack.  This  investigation  concentrated  on  the 
effect  of  cavity  flow  on  the  attached  surface  boundary-layer  characteristics  in 
adverse  pressure-gradients.  Conditions  where  flow  separation  was  present  were 
not  considered.  The  purpose  of  the  present  investigation  was  to  experimentally 
and  computationally  evaluate  the  effect  of  1,  2,  or  4  two-dimensional  spanwise 
cavities  with  square  cross-sections  on  the  boundary-layer  characteristics  in  both  a 
constant  pressure  flow  and  for  two  adverse  pressure  gradient  flows.  The 
experimental  investigation  was  conducted  in  the  San  Jose  State  University  (SJSU) 
12"  by  12"  low  speed  tunnel.  Computational  results  were  also  obtained  using  a 
numerical  solution  to  the  incompressible  Navier-Stokes  equations. 
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(a)  Flap  and  slat  with  chord  extension. 


(b)  Tabs  on  airfoil  and  flap  trailing  edges. 


(d)  Multiple  spoilers  used  to  form  multiple  cavity  concept. 


(e)  Leading  edge  cavities  to  delay  dynamic  stall. 
Figure  1 .  High  lift  concepts. 


II.  PREVIOUS  CAVITY  FLOW  RESEARCH 


Two-dimensional  boundary  layers  in  subsonic  flow  are  well  defined  as 
described  in  a  very  complete  survey  paper  by  Ligrani  [Ref.  3]  and  textbooks  such 
as  Schlichting  [Ref.  4],  White  [Ref.  5],  and  Cebeci  and  Smith  [Ref.  6].  At  low 
Reynolds  numbers,  i.e.  short  run  lengths,  in  low-turbulence  flow  the  boundary 
layer  is  usually  laminar.  For  high-lift  applications,  the  boundary  layer  rapidly 
transitions  to  a  fully  turbulent  boundary  layer  whose  mean  velocity  profile  may  be 
represented  by  the  1/7  power  law  (u/Uoo  =  (y/d)l/7).  Additional  work  on 
boundary  layers  and  the  effects  of  skin  friction,  adverse  pressure  gradient,  surface 
roughness,  and  shear-layer  flow  are  described  in  References  17  through  31. 

Nearly  40  years  ago,  Roshko  [Ref.  32]  and  Krishnamurty  [Ref.  33] 
conducted  two  of  the  first  relatively-detailed  investigations  of  cavity  flow.  These 
experiments  did  a  good  job  of  identifying  the  important  flow  characteristics 
associated  with  a  single  cavity.  In  the  present  paper,  the  cavity  streamwise 
dimension,  length,  is  designated  as  L;  the  dimension  normal  to  the  freestream, 
depth,  is  D;  and  the  distance  across  the  tunnel  and  normal  to  the  flow  direction, 
width,  is  W.  In  general,  cavity  flows  exhibit  various  steady  and  unsteady 
phenomena.  The  upstream  boundary  layer  separates  at  the  cavity  lip  to  form  a 
shear  layer  over  the  cavity.  The  shear  layer  then  reattaches  (1)  either  on  the  cavity 
floor  or  (2)  on  the  downstream  cavity  wall  or  downstream  of  the  cavity.  The 
cavity  is  considered  to  be  closed  if  the  shear  layer  reattaches  on  the  cavity  floor 
and  then  recirculates  within  the  upstream  end  of  the  cavity.  Typically  this  flow 
can  occur  when  L/D  >  4.  The  cavity  is  considered  to  be  open  when  the  shear  layer 
reattaches  near  the  rear  lip  or  downstream  of  the  cavity.  This  flow  is  typical  of 
deeper  cavities  where  L/D  <  1.  For  L/D  values  between  1  and  4,  the  d2/L  strongly 

influences  whether  a  cavity  is  open  or  closed.  At  certain  conditions  disturbance 
waves  are  fed  back  to  the  upstream  cavity  lip,  the  original  disturbance  source. 
This  feedback  loop  can  amplify  the  disturbance  waves  and  create  large  oscillating 
pressure  waves  and  noise.  In  these  circumstances  the  unsteady  flow  can  dominate 
the  cavity  flowfield. 
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A.  MEAN  FLOW  CHARACTERISTICS 


Roshko  [Ref.  32]  used  a  cavity  length  L  of  4  inches  and  systematically 
varied  the  cavity  depth  D  to  study  the  time-averaged  flow  at  low  Mach  numbers 
(0.06  to  0.20)  for  D/L  ranging  from  0.02  to  2.50.  Pressure  coefficients  on  the 
cavity  walls  and  floor  were  measured  and  the  skin-friction  coefficients  were 
calculated.  The  structure  of  the  flow  in  the  cavity  was  also  observed  and  found  to 
be  a  function  of  the  cavity  depth-to-length  ratio  D/L.  For  very  shallow  cavities 
(D/L  <  0.1)  the  shear  layer  above  the  cavity  reattached  to  the  cavity  floor  forming 
a  recirculation  region  on  the  upstream  comer  of  the  cavity.  As  the  varying  cavity 
D/L  approached  1,  the  shear  layer  attachment  location  moved  aft  along  the  cavity 
floor  to  the  rear-wall  intersection  and  up  the  wall  to  near  its  intersection  with  the 
freestream  surface.  For  a  cavity  with  a  square  cross  section  (D/L  ~  1)  there  was 
steady  flow  in  a  single  predominant  vortex  with  secondary  vortices  in  the  comers 
of  the  cavity.  For  D/L  from  1  to  2.5,  a  single  vortex  continued  to  dominate  the 
cavity  flow. 

It  was  concluded  that  the  drag  due  to  the  cavity  was  almost  entirely  due  to 
the  pressure  on  the  cavity  walls.  The  drag  contribution  from  changes  in  the  skin 
friction  coefficient  on  the  cavity  surfaces  was  relatively  small.  Further,  the  friction 
forces  on  the  cavity  walls  were  found  to  be  small  compared  with  the  pressure 
forces.  The  drag  is  analyzed  in  detail  in  ssection  II.C.  While  the  friction  forces 
due  to  the  cavity  vortex  were  negligible  compared  to  the  pressure  drag,  they  do 
play  a  role  in  determining  the  vortex  stability.  An  indication  of  this  stability  is 
shown  in  Figure  2(a)  by  the  variation  of  the  pressure  coefficient  measured  at  the 
middle  of  the  cavity  floor  as  the  cavity  depth  was  increased  from  D/L  of  0.02  to 
2.50.  Steady  pressures  were  measured  when  D/L  was  less  than  0.50  and  when  D/L 
was  between  0.87  and  2.00.  For  D/L  from  0.50  to  0.87  and  D/L  > 
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(a)  Variation  of  pressure  coefficient  measured  at  the  middle  of  the  cavity 

floor  [Ref.  32]. 

Figure  2.  Effect  of  cavity  D/L  ratio  on  local  pressure-coefficients. 
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2.00  the  pressures  were  unsteady.  The  change  from  unstable  flow  to  stable  flow  as 
D/L  was  varied  through  0.87  was  quite  distinct.  There  were  two  stable  states  for 
this  range  of  D/L.  When  the  cavity  vortex  was  steady  (D/L  <  0.50  and  0.87  <  D/L 
<  2.00),  the  flow  over  the  cavity  was  also  steady. 

A  single,  stable  vortex  was  formed  by  the  deflection  of  a  portion  of  the 
shear  layer  at  the  downstream  cavity  edge  into  the  square  cross-section  cavity. 
The  relatively  high  pressure  on  the  cavity  wall  in  that  vicinity  accounts  for  most  of 
the  drag.  The  pitot  pressure  at  the  top  back  comer  is  shown  in  Figure  2(b).  This 
pressure  measurement  is  analogous  to  a  Preston-tube  measurement  and  it  is  related 
to  the  local  skin-friction  coefficient.  It  is  intended  to  give  a  measure  of  the 
pressure  near  the  top  of  the  rear  cavity  wall.  There  is  an  unsteadiness  shown  for 
D/L  between  0.50  and  0.87.  Above  D/L  of  1  there  is  a  hysteresis  which  shows  that 
the  pressure  coefficient  is  dependent  upon  whether  D  is  increasing  or  decreasing. 
For  the  square  cavity  (L/D  =1),  the  skin-friction  coefficient  was  calculated  from 
the  boundary-layer  profiles  with  the  following  results: 


Boundary-Layer  Profile  Location 

cf 

I  -  0.375  L  upstream  of  the  cavity  front  wall 

0.0015 

II  -  0.125  L  downstream  of  the  cavity  rear  wall 

0.0011 

III  -  0.375  L  downstream  of  the  cavity  rear  wall 

0.0012 

IV  -  0.375  L  downstream  of  the  cavity  rear  wall 
location  without  the  cavity  in  place 

0.0014 

Another  investigation  by  Maull  and  East  [Ref.  34]  studied  cavities  at 
low  subsonic  velocities  using  oil-flow  and  surface  static-pressure  distributions. 
They  found  that  the  flow  steadiness  depended  on  cavity  width  as  well  as 
streamwise  length  and  depth.  For  a  cavity  width-to-length  ratio  W/L  of  9,  a  non- 
uniform  spanwise  variation  of  pressure  coefficient  was  most  notable  below  D/L  of 
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0.85  and  between  D/L  of  1.5  to  over  2.  These  regions  are  roughly  consistent  with 
the  steadiness  of  Roshko's  pressure  coefficient  data  in  Figure  2.  There  were 
regions  of  uniform  pressure  distribution  near  D/L  of  1  and  2.5.  There  was  a  very 
rapid  switch  from  non-uniform  pressure-coefficient  variation  for  D/L  =  0.85  to  a 
nearly  uniform  pressure-coefficient  distribution  for  D/L  =  0.86.  These  results 
appeared  to  also  depend  on  the  ratio  of  boundary-layer  thickness  to  cavity  length 
5/L. 

Rossiter  [Ref.  35]  investigated  the  flow  over  rectangular  cavities  at 
subsonic  and  transonic  speeds.  The  investigation  was  conducted  in  a  2'  by  1.5' 
transonic  tunnel.  The  cavity  had  a  4"  width,  an  8"  length,  and  a  depth  which  was 
varied  from  0.8  to  8.0  inches.  Steady  and  unsteady  pressures  were  recorded  on  the 
cavity  floor  and  downstream  of  the  cavity.  While  the  flow  was  highly  unsteady,  it 
was  useful  to  briefly  discuss  the  nature  of  the  time-average  flow  over  cavities.  For 
very  shallow  cavities,  the  flow  over  the  front  and  rear  walls  may  be  considered 
independently  as  the  flow  down  and  up  a  step,  respectively.  The  airflow  will 
separate  from  the  front  edge  and  reattach  at  some  point  along  the  floor  of  the 
cavity.  The  pressure  in  the  separated  region  will  be  lower  than  the  freestream 
pressure.  This  occurs  because  the  freestream  flow  speeds  up  as  it  enters  the  cavity. 
Then  the  pressure  rises  at  the  attachment  point.  As  the  airflow  approaches  the  rear 
wall,  it  is  slowed  and  then  the  pressure  increases  until  a  position  is  reached  where 
the  boundary  layer  separates  to  form  the  boundary  layer  ahead  of  the  rear  wall. 
The  boundary  layer  will  usually  reattach  at  a  location  downstream  of  the  cavity. 
As  the  depth-to-length  ratio  of  the  cavity  increases,  the  attachment  and  separation 
points  on  its  floor  will  move  closer  together  until  a  reverse  flow  develops  between 
the  high  pressure  region  ahead  of  the  rear  wall  and  the  low  pressure  region  behind 
the  front  wall.  A  large  vortex  then  forms  within  the  cavity. 

As  shown  by  Roshko  [Ref.  32]  the  mean  flow  pattern  depends  on  the 
length-to-depth  ratio  L/D  of  the  cavity.  Rossiter  compared  his  data  with  Roshko's 
as  shown  in  Figure  3.  The  results  are  not  directly  comparable  between  these  tests 
because  of  different  length-to-depth  ratios  and  the  relatively  thicker  boundary 
layer  approaching  the  cavity  in  Rossiter's  investigation.  For  the  very  shallow 
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(b)  Pressure  coefficient  measured  at  the  top  back  corner  of  a  cavity  for  a 
range  of  depth-to-length  ratios  [Ref.  32]. 

Figure  2.  Concluded. 
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cavities  (L/D  >  8)  the  two  pressure  rises  associated  with  flow  attachment  and  with 
separation  on  the  floor  of  the  cavities  may  be  seen.  Immediately  behind  the  cavity 
front  wall  the  pressure  reaches  a  low  value.  The  pressure  increases  from  the  center 
of  the  closed  vortex  in  the  corner  at  the  front  wall  to  the  rear  end  of  the  vortex. 
The  pressure  levels  off  and  then  increases  again  as  the  downstream  wall  is- 
approached.  At  a  length-to-depth  ratio  of  8,  the  two  pressure  rises  have  merged 
and  extend  from  x/L  ~  0.4  to  0.9.  At  a  length-to-depth  ratio  of  6,  the  pressure  is 


Figure  3.  Comparison  of  Pressure  coefficients  measured  on  cavity  floors  and  on 
the  surface  downstream  of  the  cavity  by  Rossiter  [Ref.  35]  and  Roshko  [Ref.  32]. 


almost  constant  (Cp  ~  -0.02)  along  the  floor  of  the  cavity  indicating  that  the  two 

separations  have  combined.  For  the  deeper  cavities  (L/D  of  4  and  2)  there  is  a 
decrease  in  the  pressure  on  the  floor  associated  with  the  high  airspeeds  at  the 
periphery  of  the  cavity  vortex.  As  Mach  number  is  increased,  the  flow  attachment 
point  on  the  floor  of  the  shallower  cavities  moves  downstream  so  that  the  cavities 
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become  effectively  deeper.  As  a  result,  the  pressure  distributions  for  the  shallower 
cavities  show  a  large  variation  with  Mach  number  whereas  the  pressure 
distributions  in  the  deeper  cavities  are  comparatively  independent  of  Mach 
number. 

Other  investigations,  such  as  Charwat  et  al  [Ref.  36  and  37],  have  studied 
cavities  in  supersonic  flow  and  found  that  the  same  general  flows  exist  as  those 
found  at  subsonic  speeds.  Three-dimensional  cavities  which  relate  to  landing  gear 
or  weapons  bays  were  studied  by  Plentovich  [Ref.  38].  Greater  unsteadiness  was 
found  at  the  highest  Reynolds  number  and  the  pressure  distributions  were  found  to 
be  sensitive  to  the  thickness  of  the  boundary  layer  entering  the  cavity.  Another 
application  of  cavities  is  found  in  aircraft  where  telescopes  are  used  for  various 
observations.  The  objective  of  the  investigation  by  Buell  [Ref.  39]  was  to 
minimize  the  shear-flow  disturbance  using  an  antiresonance  device.  Devices  were 
developed  which  suppressed  cavity  shear-flow  resonance  and  produced  thinner 
shear  layers.  This  shear-layer  flow  suppression  enabled  better  light  transmission 
to  the  telescope.  However,  the  suppression  of  resonance  also  made  the  boundary 
layer  downstream  of  the  cavity  more  susceptible  to  flow  separation. 


B.  FLOW  INDUCED  OSCILLATIONS 


There  have  been  reviews  by  Rockwell  and  Naudascher  [Ref.  40]  and  by 
Komerath,  Ahufa,  and  Chambers  [Ref.  41]  which  emphasize  cavity  flow-induced 
oscillations.  Both  papers  used  the  classification  of  self-sustaining  flow  oscillation 
over  cavities  first  presented  by  Rockwell  and  Naudascher  and  reproduced  in 
Figure  4  from  reference  40.  Three  flow-interaction  categories  were  identified:  (1) 
fluid  dynamic,  (2)  fluid  resonant,  and  (3)  fluid  elastic.  In  many  situations  more 
than  one  of  these  interactions  may  be  involved.  The  fluid-dynamic  interactions 
involve  coupling  between  oscillations  of  the  shear  layer  over  the  cavity  with  the 
flow  inside  the  cavity.  Many  of  the  oscillating  flows  at  low  speeds  over  shallow 
cavities  (L/D  >1)  fall  under  this  category.  The  mechanisms  involved  are  believed 
to  arise  from  shear-layer  instability  and  vortex  shedding.  Large-scale  coherent 
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structures  present  in  the  shear  layer  are  known  to  play  the  major  role  in  such 
interactions.  Fluid-resonant  interactions  are  flow  oscillations  which  appear  to  be 
controlled  by  the  acoustic  modes  of  the  cavity.  These  are  usually  encountered  in 
cavities  which  have  large  depths  normal  to  the  flow  direction  (L/D  <  1).  These 
fluid-resonant  oscillations  are  observed  in  flows  over  cavities  at  high  Mach 
numbers.  Fluid-elastic  flows  involve  interactions  between  the  shear  layer  over  the 
cavity  and  the  elastic  boundaries  of  the  cavity.  These  interactions  may  cause 
sonic-fatique  problems. 


Figure  4.  Classification  of  cavity  flows  by  Rockwell  and  Naudascher  [Ref.  40]. 

Krishnamurty  [Ref.  33]  varied  the  rectangular  cavity  L/D  ratio  in  a  study  of 
the  sound  radiated  by  cavity  flow.  The  cavity  had  a  constant  depth  of  0.1"  and 
was  located  in  a  flat  plate.  A  trip  wire  and  the  flat-plate  angle-of-attack  were 
varied  to  change  the  boundary  layer  from  laminar  to  turbulent.  The  cavity  length 
was  varied  from  0  to  2  inches.  The  investigation  was  conducted  in  a  blowdown 
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wind  tunnel  over  a  Mach  number  range  from  0.4  to  1.5.  Data  included  hot  wire 
measurements  of  the  boundary  layer,  magnitude  and  directionality  of  the  radiated 
sound,  as  well  as,  schlieren  and  interferometer  flow  visualization.  Below  a 
minimum  cavity  length  it  was  found  that  the  shear  layer  flowed  over  the  cavity  and 
reattached  downstream  of  the  rear  cavity  wall.  For  these  conditions  there  were  no 
acoustic  oscillations. 

Above  this  cavity  length,  the  acoustic  frequency  was  inversely  proportional 
to  the  cavity  length  f  =  a  /  L  .  The  constant  of  proportionality  a  was  different 
for  laminar  and  turbulent  boundary  layers.  This  dependence  is  presented  in  Figure 
5  for  both  laminar  and  turbulent  boundary  layers.  While  there  was  a  dominant 
frequency  for  laminar  flow,  there  were  two  frequencies,  low  and  high,  for  the 
turbulent  case.  Krishnamurty  used  Strouhal  number  to  obtain  a  dimensionless 
frequency.  If  the  freestream  velocity  is  used  as  the  characteristic  velocity,  the 
Strouhal  number  becomes 


Ns  =  fL/U00=a/U00 


(1) 


At  low  Mach  numbers  there  was  little  sound  radiation  directionality;  at  increased 
Mach  numbers  the  sound  became  more  intense  and  directional.  The  radiation 
pattern  was  observed  using  a  schlieren  system.  It  was  shown  that  the  unsteady 
pressure  oscillations  were  greater  in  a  laminar  boundary  layer  than  in  a  turbulent 
boundary  layer. 

Rossiter  [Ref.  35]  found  that  the  unsteady  pressures  contain  both  random 
and  periodic  components.  The  random  component  predominates  in  the  shallower 
cavities  (L/D  >  4)  and  was  most  intense  near  the  rear  wall.  The  unsteady  pressure 
had  a  smooth  spectrum  over  a  broad  frequency  band  showing  the  random  character 
in  shallow  cavities.  For  very  shallow  cavities,  a  local  region  of  intense  pressure 
fluctuations  was  found  where  the  shear  layer  flow  attaches  to  the  floor  of  the 
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(a)  Proportionality  constant  a. 

Figure  5.  Effect  of  Mach  number  on  the  acoustic  field  due  to  a  cavity. 
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(b)  Strouhal  number  N$. 
Figure  5.  Concluded. 
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cavity.  There  were  usually  two  peaks  of  almost  equal  magnitude.  The  periodic 
component  predominates  in  the  deeper  cavities  (L/D  <  4)  and  may  form  standing 
wave  patterns  with  one  peak  whose  magnitude  is  much  larger  than  any  other.  It  is 
suggested  that  the  periodic  component  is  due  to  an  acoustic  resonance  within  the 
cavity  excited  by  a  phenomenon  similar  to  that  causing  edge-tones.  The  periodic 
pressure  fluctuations  may  be  very  large.  Mean  values  up  to  0.35  times  the 
freestream  dynamic-pressure  were  measured.  These  results  indicated  the 
predominance  of  the  periodic  component  over  the  random  oscillation.  Increasing 
Mach  number  caused  an  increased  periodic  component.  The  periodic  component 
was  reduced  in  the  thicker  boundary-layer.  This  suggests  that  the  lack  of  periodic 
pressure  fluctuations  in  shallow  cavities  is  partially  due  to  a  large  ratio  of 
boundary  layer  thickness  to  cavity  depth  5/D.  This  result  suggested  a  means  to 
reduce  the  magnitude  of  the  pressure  fluctuations.  The  suppression  was  achieved 
using  a  small  spoiler  located  upstream  of  the  cavity  to  effectively  increase  the 
boundary-layer  thickness. 

In  this  investigation,  Rossiter  identified  experimentally  the  principal 
characteristics  of  the  periodic  pressure  fluctuations  as  the  following: 

(1)  The  pressure  fluctuations  may  contain  a  number  of  periodic 
components  or  modes  which  occur  at  specific  frequencies  and  are 
designated  as  m  =  1,  2,  3,. . . 

(2)  The  frequency  of  any  component  is  inversely  proportional  to  the  cavity 
length  and  increases  with  freestream  velocity.  Dimensionally  this  suggests 

f(j)  =  g(M,Re)  (2) 

Over  the  Reynolds  number  range  investigated  it  was  found  that  the  effect  of 
cavity  length  was  small. 

(3)  Shadowgraphs  showed  that  the  pressure  fluctuations  are  accompanied 
by  the  periodic  shedding  of  vortices  from  the  front  lip  of  the  cavity  while 
the  principal  acoustic  source  is  close  to  the  rear  lip  of  the  cavity. 
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(4)  Under  some  conditions,  standing  wave  patterns  occur  in  the  cavity 
which  are  probably  due  to  an  acoustic  resonance  within  the  cavity. 


The  regular  shedding  of  vortices  found  in  flow  over  cavities  is  a  feature  of 
the  edge  tone  phenomenon  as  well  as  the  strong  acoustic-radiation.  A  connection 
was  assumed  between  these  two  acoustic  features.  Further  it  was  assumed  that  the 
acoustic  radiation  initiates  the  vortex  shedding  and  that  the  passage  of  the  vortices 
over  the  rear  lip  of  the  cavity  is  responsible  for  the  acoustic  radiation.  If  the 
average  speed  of  the  vortices  over  the  cavity  is  K  times  the  freestream  velocity  and 
sound  waves  travel  upstream  in  the  cavity  at  a  mean  velocity  c  then 


f_KUTO_  c 


(3) 


where  Ay  was  assumed  to  equal  Aa.  The  phase  relation  between  the  two  motions  is 
unknown. 


A  solution  was  found  by  identifying  the  particular  phase  of  the  acoustic 

radiation  which  reaches  the  front  lip  of  the  cavity  at  the  instant  that  a  vortex  is 
shed.  It  is  assumed  that  a  vortex  is  y  Ay  behind  the  rear  lip  when  this  particular 

phase  of  the  acoustic  radiation  leaves  the  source  at  the  rear  lip.  It  was  assumed 
that  there  were  mv  complete  wavelengths  of  the  vortex  motion  and  ma  complete 

wavelengths  of  acoustic  radiation.  Further  it  was  assumed  that  there  is  a  time 
interval  t'  between  arrivals  of  an  identified  phase  of  the  acoustic  radiation  at  the 

front  lip  just  as  a  vortex  is  shed.  Therefore  the  vortex  pattern  moves  downstream  a 
distance  KU^t'  in  this  time  interval  so  that 

mvAv  =L  +  yAv  +  KU^t'  (4) 

Also  in  time  t'  the  internal  wave  system  has  moved  a  distance  ct'  so  that 

L  =  maAa  +  ct'  (5) 

Then  t'  is  eliminated  between  equations  4  and  5.  The  freestream  velocity  was 

replaced  by  Ma;  further,  it  was  assumed  the  acoustic  velocities  internal,  c,  and 
external,  a,  to  the  cavity  were  equal;  and  that  the  total  number  of  wavelengths,  mv 
+  ma  ,  equal  the  mode  number,  m. 
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These  assumptions  led  to  a  physical  model  of  the  flow  which  may  be 
represented  by  the  empirical  frequency  equation 


f  = 


u.  (m-y) 

L  (Um 

u  ^ 


(6) 


where  K  equals  0.57,  m  is  an  integer,  and  y  equals  0.25.  When  the  frequency  of 
one  of  the  components  (m-y)  is  close  to  the  natural  frequency  of  the  volume  of  air 
in  the  cavity,  resonance  occurs.  Both  cavity  length  and  depth  influence  which 
mode  number  is  dominant  for  a  particular  cavity  geometry. 


Additional  investigations  of  flow  induced  cavity  pressure 
oscillations  and  acoustic  resonance  to  refine  Rossiter's  work  were  conducted  by 
East  [Ref.  42],  Heller  et  al  [Ref.  43],  Bilanin  and  Covert  [Ref.  44],  Miles  [Ref.  45], 
Block  [Ref. 46],  Sarohia  [Ref.  47],  Yu  [Ref. 48],  and  Gharib  and  Roshko  [Ref.  49]. 
The  investigation  by  Tam  and  Block  [Ref.  50]  obtained  experimental  data  for  the 
tone-frequency  characteristics  for  flow  Mach  numbers  less  than  0.4.  It  was 
observed  that  there  was  a  transition  from  the  normal  mode  resonance  mechanism 
to  a  feedback  instability  mechanism  for  disctete  tone  generation  as  flow  Mach 
number  increases.  A  mathematical  model  of  acoustic-feedback  oscillations  was 
developed  which,  in  contrast  to  the  Bilanin  and  Covert  model,  accounted  for  the 
shear-layer  thickness.  Compared  with  the  Rossiter  model,  the  Tam  and  Block 

model  added  two  parameters;  (1)  the  cavity  length-to-depth  ratio  L/D;  and  (2)  the 
ratio  of  the  momentum  thickness  of  the  shear  layer  to  the  cavity  length  82/L. 


The  acoustic -wave  generation  process  sketched  in  Figure  6  (from  Ref.  51) 
for  a  supersonic  freestream  was  adopted.  Flow-induced  cavity  oscillations  are 
caused  by  the  interaction  of  the  free  shear  layer  and  the  complex  internal  cavity 
wavetrains.  The  shear  layer  oscillated  up  and  down  near  the  trailing  edge  of  the 
cavity.  During  the  upward  motion  of  the  cycle,  the  fluid  of  the  shear  layer  shields 
the  trailing  edge  of  the  cavity  from  the  external  flow  and  the  predominant  flow  is 
over  the  cavity  with  no  pressure  waves.  When  the  shear  layer  deflects  downward , 
there  is  flow  of  the  freestream  into  the  cavity  which  causes  a  transient  high 
pressure  region  near  the  cavity  trailing  edge,  which  forces  propagation  of  a 
compression  wave  in  all  directions.  The  convection  effect  of  the  freesteam 
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modifies  the  shape  of  the  wave  front  as  it  radiates  away  from  the  cavity  trailing- 
edge. 


In  Figure  6,  the  essential  features  of  a  typical  oscillation  cycle  are  divided 
into  six  time-sequential  parts,  identified  by  the  letters  A  through  F.  The  external- 
flow  Mach  waves  represent  a  frees tream  Mach  number  of  1.5.  Each  part  sketched 
in  Figure  6  should  be  viewed  as  a  typical  phase  of  the  oscillation  cycle.  The 
direction  of  the  motion  of  each  wave  is  indicated  by  an  arrow.  The  exact  wave 
representation  depends  on  the  cavity  geometry  (L,  D,  and  W),  the  external  Mach 
number,  the  boundary  layer  characteristics  at  the  upstream  end  of  the  cavity,  and 
the  freestream  turbulence.  The  starting  point  for  the  cycle  is  arbitrarily  selected.  It 
is  helpful  to  review  the  entire  cycle  to  gain  a  good  understanding  of  cavity  flow 
phenomena.  This  cycle  is  now  discussed  starting  at  A: 

(A)  The  pressure  wave  from  the  previous  trailing-edge  disturbance  reaches 
the  front  of  the  cavity  and  reflects.  Another  wave,  moving  from  the  front  wall 
approaches  the  rear  wall.  The  shear  layer  is  above  the  cavity  trailing  edge,  so  the 
external  flow  does  not  produce  a  disturbance  at  the  cavity  trailing  edge.  Some 
fluid  leaves  the  cavity  at  the  rear. 

(B)  The  shear-layer  waveform  travels  rearward  and  reduces  the  height  of 
the  shear  layer  above  the  trailing  edge.  A  new  compression  wave  begins  to  form 
at  the  rear  as  the  flow  interacts  with  the  trailing  edge  and  fluid  is  added  to  the 
cavity.  The  front  compression  wave  has  reflected  off  the  front  wall  and  moves 
rearward  nearly  in  phase  with  the  shear  layer  displacement.  The  previous 
rearward  wave  has  reached  the  trailing  edge. 

(C)  The  wave  reflected  off  the  front  wall  continues  to  move  rearward  in 
phase  with  the  shear-layer  displacement.  The  shear  layer,  which  is  now  below  the 
trailing  edge  at  the  rear  of  the  cavity,  forms  a  new  forward-traveling  compression 
wave  as  the  external  flow  impinges  on  the  back  of  the  cavity. 

(D)  The  newly-generated,  forward-traveling  compression  wave  and  the 
reflected,  rearward-traveling  compression  wave  meet  and  interact  near  the  cavity 
center. 
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Figure  6.  Typical  oscillation  cycle  for  cavity  flow  [from  Ref.  51]. 
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(E)  After  interaction,  the  waves  continue  in  their  respective  directions.  The 
external  part  of  the  forward-traveling  wave  moves  into  the  supersonic  flow,  thus 
causing  it  to  be  tipped  more  than  the  external  flow  angle.  The  rearward  wave 
moves  in  the  same  direction  as  the  external  flow  and  travels  at  a  subsonic  speed 
relative  to  it.  This  subsonic  relative  speed  explains  why  the  rearward-traveling 
wave  stops  at  the  shear  layer.  At  the  rear,  the  shear  layer  reaches  the  trailing  edge 
height. 

(F)  The  shear  layer  is  now  above  the  trailing-edge  height.  The  wave 
generated  at  the  trailing  edge  approaches  the  front  of  the  cavity,  and  the  reflected 
wave  nears  the  rear  of  the  cavity.  The  next  step  is  the  same  as  (A),  and  the 
oscillation  cycle  repeats. 

For  subsonic  flow,  the  process  is  essentially  the  same,  particularly  as  regards  the 
internal  wave  structure.  The  forward-traveling  wave  will  still  be  supersonic  with 
respect  to  the  external  flow.  The  external  wave  structure  will  usually  be  non¬ 
existent.  Instead,  the  shear  layer  will  tend  to  roll  into  transverse  vortices  with  the 
number  dependent  primarily  on  the  cavity  geometry  (L/D))  and  on  the  freestream 
Mach  number.  Typically  there  will  be  either  two  (mode  II)  or  three  (mode  III) 
vortical  structures  present. 

For  the  purpose  of  computing  the  phases  of  the  acoustic  waves  generated  at 
the  trailing  edge  of  the  cavity,  Tam  and  Block  [Ref.  50]  made  the  size  of  the  noise 
source  very  small.  A  periodic  line  source  was  used  to  produce  this  flow  pattern. 
Inside  the  cavity  the  mean  velocity  was  zero.  The  mathematical  model  of  the 
effect  of  the  interaction  of  the  acoustic- wave  field  on  the  instabilities  of  the  shear 
layer  was  assumed  to  be  convective.  A  mean  shear-layer  momentum  thickness 
was  used  in  the  model.  The  resultant  equations,  which  are  fairly  complex,  are 
presented  in  Reference  50.  This  model  provided  good  agreement  between 
predicted  discrete-tone  frequencies  and  Rossiter's  data  (0.4  <  M  <  1.2)  and  Tam  & 
Block's  data  (M  <  0.2).  For  the  very  low  Mach  numbers,  M  <  0.2,  it  appeared  that 
the  tones  were  generated  by  the  normal-mode  resonance  mechanism.  The 
transition  between  the  feedback  mechanism  and  the  normal-mode  resonance  was 
rather  gradual.  This  suggested  that  the  Tam  and  Block  model  may  provide  the 
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basis  for  a  unified  model  of  the  flow-induced  tone  phenomenon  to  include  cavity- 
tone  frequencies  throughout  the  Mach  number  range. 

Ziada  and  Rockwell  [Ref.  52]  found  that  impingement  of  mixing  layers  on 
solid  boundaries  enhanced  the  organization  of  all  harmonic  modes  in  a  shear  layer. 
This  effect  was  investigated  by  Rockwell  and  Knisely  [Ref.  53}  specifically  for 
cavity  flow  using  a  laser  velocimeter.  As  shown  in  Figure  7,  inserting  a 
downstream  cavity-impingement  edge  makes  a  large  change  in  the  unsteady  nature 
of  the  shear  flow.  Autocorrelations  were  taken  along  the  edge  of  the  shear  layer 
where  u/Uqo  =  0.95  using  an  average  of  six  sequential  samples.  These  results  are 

presented  in  Figure  7;  note  that  the  ordinate  scale  changes  for  each  curve.  Both 
data  sets  were  taken  at  the  same  shear  layer  location  of  x/80  =  134.  For  the 

backward  facing  cavity,  there  were  large  variations  in  the  autocorrelations.  The 
typical  period  of  oscillation  is  quite  irregular  and  tends  to  increase  with 
downstream  location.  For  the  cavity  flow,  the  enhanced  organization  of  the  flow 
is  quite  dramatic.  In  addition,  there  are  no  apparent  variations  in  the  oscillation 
frequency.  This  enhancement  of  the  shear  layer  organization  extended  along  the 
entire  length  of  the  shear  layer.  This  finding  substantiates  the  concept  of 
disturbance  feedback.  The  perturbations  propagated  upstream  from  the 
impingement  surface  affected  the  amplitude  of  the  locally-induced  pressure 
fluctuations  and  the  resultant  shear  layer  deformation. 

A  few  examples  of  the  many  additional  cavity-flow  investigations  are  cited 
below  to  illustrate  some  results  pertinent  to  the  present  investigation.  Franke  and 
Carr  [Ref.  54]  investigated  the  effect  of  geometry  on  modification  of  flow-induced 
oscillations  for  open  cavity  flow  in  the  Mach  number  range  from  0.6  to  3.3.  It  was 
found  that  ramps  at  the  leading  and  trailing  edges  of  the  rectangular  cavities  could 
be  effective  in  reducing  pressure  amplitudes  under  some  conditions.  Some  of  the 
cavity  configurations  studied  by  Franke  and  Carr  included  two  cavities  in  tandem. 
Tandem  cavities  were  also  investigated  by  Zhang  and  Edwards  [Ref.  55]  in 
supersonic  flow.  They  examined  both  the  mean  and  unsteady  flowfields.  It  was 
found  that  two  cavities  whose  L/D  was  either  1  or  3  had  little  effect  on  each  other. 
In  contrast,  a  L/D  =  3  cavity  completely  altered  the  flow  in  a  downstream  L/D  =  1 
cavity. 
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Betts  [Ref.  56]  experimentally  evaluated  the  effect  of  slotted  walls  on  flow 
in  a  rectangular  water-channel.  It  was  found  that  the  effect  of  geometry  on  flow 
oscillations  could  be  modeled  by  empirical  relations.  Sarno  and  Franke  [Ref.  57] 
evaluated  several  devices  located  upstream  of  a  cavity  which  attempted  to  reduce 
flow  oscillation.  Like  Buell  [Ref.  39]  it  was  found  that  stationary  fences  at  the 
cavity  leading  edge  were  the  most  effecive  flow-suppression  device  since  they 
effectively  reduced  L/8<> 
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Figure  7.  Comparison  of  autocorrelations  taken  with  and  without  the  cavity 
impingement  edge  at  corresponding  locations  in  the  cavity  shear-layer  [Ref.  62]. 
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C.  CAVITY  DRAG 


The  drag  of  cavities,  holes,  and  gaps  were  discussed  by  Hoerner  [Ref.  13] 

in  his  book  on  fluid-dynamic  drag.  Cavity  drag  is  relatively  small  compared  to  the 
overall  aircraft  drag.  The  cavity  drag  coefficient  Cp  was  typically  less  than  0.01 

based  on  a  reference  area  equal  to  the  cavity  opening,  LW.  For  comparison, 
complete  aircraft  drag  coefficient  is  usually  based  on  wing  area;  the  cavity  opening 
area  LW  of  landing  gear  or  a  stores  cavity  is  usually  at  least  one  or  two  orders  of 
magnitude  less  than  wing  area  so  that  the  Cp  contribution  to  aircraft  drag 

coefficient  is  of  the  order  of  0.001  to  0.0001  or  less.  The  shape  of  the  edges  of  the 

cavities  have  a  major  influence  on  the  magnitude  of  the  drag  as  shown  in  Figure  8. 
For  square  cavities,  L/D  =  1  such  as  those  used  in  the  present  investigation,  Cp 
was  0.0083.  When  the  slope  on  the  top  of  the  rear  wall  was  changed  to  -5°,  Cp 
was  reduced  to  0.0072.  Rounding  the  rear  slope  reduced  Cp  to  0.0060.  The 
other  sketches  in  Figure  8  show  a  cavity  with  the  same  opening  area  but  a  larger 
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Figure  8.  Cavity  drag  coefficient  Cp  measured  for  several  edge  shapes  [Ref.  13]. 
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cavity  volume  with  variations  in  the  slopes  into  and  out  of  the  cavity.  It  is  shown 
that  the  value  of  Cp  can  be  varied  from  0.0070  to  as  high  as  0.0400  by  changing 

these  slopes.  Cavities  with  front  slopes  which  direct  freestream  flow  into  the 
cavity  and  with  rear  wall  slopes  which  inhibit  flow  out  of  the  cavity  greatly 
increase  the  drag.  For  minimum  drag,  the  upstream  edge  should  be  nearly  parallel 
with  the  freestream  and  the  downstream  edge  should  be  bent  at  a  small  angle  into 
the  cavity. 

The  maximum  shear  stress  increases  with  the  cavity  L/D.  Gharib  and 
Roshko  [Ref.  49]  measured  the  maximum  Cf  as  0.013.  This  Cf  value  compares 

closely  with  0.0125  measured  by  Liepmann  and  Laufer  [Ref.  58]  for  a  two- 

dimensional  shear  layer  and  0.012  measured  by  Kistler  and  Tan  [Ref.  59]  for  a 
two-dimensional  cavity  shear  layer.  The  maxima  of  cavity  Cf  profiles  did  not 

have  a  constant  value.  In  contrast,  Cf  was  found  to  be  constant  in  two- 
dimensional,  self-similar  free  turbulent  shear  layers. 

Cavity  drag  is  defined  as  the  net  force  in  the  flow  direction  experienced  by 
the  cavity.  The  cavity  control  volume  which  can  be  used  for  the  momentum 

balance  is  shown  in  Figure  9(a).  Surface  1  is  the  open  interface  between  the  cavity 
and  freestream  flows  where  drag  equals  the  momentum  flux  drag  Djyj  occurs.  On 
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(a)  Schematic  of  the  cavity  and  its  control  volume. 

Figure  9.  Cavity  drag  and  effect  of  cavity  length  on  maximum  Cp. 
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On  the  upstream  and  downstream  cavity  walls  (surfaces  4  and  2,  respectively)  the 
pressures  are  integrated  to  give  a  pressure  drag  Dp.  On  the  cavity  floor  (surface 

3)  the  shear  stress  gives  a  third  drag  force  contribution  Dg.  These  terms  are 

summed  to  zero  for  the  control  volume.  The  drag  force  on  walls  2,  3,  and  4  is 
given  by 

D  =  Dp  +  Dg  =  Dm  (7) 


The  cavity  drag  is  expressed  as  either  the  direct  forces  on  the  cavity  walls  and 
floor  or  the  cavity  momentum  flux.  The  drag  terms  in  equation  7  may  be 
expanded  to  provide  the  cavity  momentum-drag  balance  in  integral  form  as 


D  =  J[2]  PdA2  -  Jf4l  PdA4  +  Jf3l  V*A3  =  |m  [  ~PU  V  ~  PU'  V1  +p 
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Roshko  [Ref.  32]  found  that  the  dynamic  pressure  in  the  recirculating  cavity 
flow  is  low.  On  their  axisymmetric  experimental  cavity  model,  which  had  no 
cavity  floor  at  the  centerline,  there  was  a  standing  axisymmetric  vortex.  At  the 
body  centerline  the  vortex  velocity  components  tended  to  cancel  and  the  resultant 

velocity  was  nearly  zero.  As  a  consequence,  the  resultant  drag  due  to  shear  stress, 
D§  ,  was  negligible  compared  to  the  pressure  drag.  Dp.  This  simplified  the  drag 

equation  to 

D  =  Dp  =  DM  (9) 


This  drag  force  can  be  nondimensionalized  as  a  drag  coefficient  using  the 
freestream  dynamic  pressure  and  an  arbitrary  reference  area  so  that 
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Alternatively,  the  reference  area  may  be  the  cavity  opening  Sref  =  LW  to  obtain 
an  equivalent  average  friction  coefficient  due  to  the  cavity  (or  cavity  drag 
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coefficient).  Both  coefficients,  Cp  and  Cp,  give  the  same  drag  value.  Gharib  and 
Roshko  [Ref.  59]  choose  to  associate  Cp)  with  the  forces  on  the  cavity  walls  and 
Cp  with  cavity  momentum  flux.  For  cavities  in  the  non-oscillation  mode  Cp  was 
near  0.0001;  in  the  self-sustained  mode,  the  Cp  was  less  than  0.01;  and  while  in 
the  wake  mode  (L/82  >  155  or  L/D  >  1.25)  the  Cp  was  between  0.06  and  0.08. 
For  comparison,  the  estimated  Cp  based  on  boundary  layer  friction  in  the  absence 
of  a  cavity  was  about  0.005. 

Measurement  of  the  pressure  distribution  on  the  solid  surfaces  of  the  cavity 

provides  an  understanding  of  drag  generation.  In  short  cavities  with  non¬ 
oscillating  flow,  the  pressure  coefficient  Cp  on  both  the  upstream  and  downstream 

walls  tends  to  be  slightly  negative  and  uniform  indicating  essentially  no  drag 
contribution.  In  cavities  with  flow  oscillations  (84  <  L/82  <  155  or  0.68  <  L/D  < 

1.25)  the  Cp  on  both  upstream  and  downstream  walls  approaches  -0.1,  except 
near  the  freestream  edge  of  the  rear  wall.  The  Cp  on  the  rear  wall  in  the  region 

from  0.1  D  to  the  freestream  surface  increases  to  a  positive  value  equal  to  a  few 
tenths.  This  small  region  accounts  for  the  entire  pressure-drag  force.  For  cavities 
in  the  wake  mode,  shear-layer  impingement  occurs  further  below  the  edge  of  the 
rear  wall.  As  a  result  the  positive  pressure  coefficient  extends  much  further  into 
the  cavity.  For  large  L/D,  impingement  even  occurs  on  the  cavity  floor  leading  to 
positive  Cp  values.  On  the  upstream  wall,  the  negative  Cp  is  nearly  twice  that 

for  self-sustaining  oscillation  modes. 

Examples  of  selected  cp  values  from  Gharib  and  Roshko  [Ref.  49]  are 

presented  in  Figure  9(b),  The  plot  shows  the  variation  of  maximum  pressure 
coefficient  with  increasing  cavity  length.  The  cp  on  the  upstream  comer  of  the 

cavity  is  nearly  constant  (  ~  -0.07)  in  the  self-sustained  oscillation  mode  II  and  III 

regimes  and  are  more  negative  (  ~  -0.10  to  -0.15)  in  the  wake  regime.  In  the  self- 
sustaining  oscillation  regimes,  the  maximum  Cp  occurs  at  the  downstream  comer 

as  shown  by  the  overlap  of  the  plus  and  solid-circle  symbols.  As  cavity  length 
increases  in  the  wake-mode  flow  regime,  the  peak  Cp  on  the  downstream  moves 

into  the  cavity  from  the  downstream  comer.  This  is  consistent  with  shear-layer 
impingement  inside  the  cavity,  as  discussed  earlier. 
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Gharib  and  Roshko  [Ref.  49]  used  a  laser  velocimeter  to  measure  both  the 
turbulent  and  mean  shear  layer  velocity  components.  These  proved  to  be  difficult 
to  obtain  especially  in  the  wake  mode.  The  integral  of  uv  across  the  cavity 

opening  represents  a  relatively  small  difference  between  two  large  contributions. 
It  was  found  that  the  variation  of  Cp  while  showing  a  lot  of  scatter,  was 

generally  consistent  with  Cp  p.  In  conclusion,  it  was  found  that  the  oscillating 
shear  layer  in  the  cavity  is  fundamentally  different  from  a  free  shear  layer.  This 
difference  is  attributable  to  coupling  between  the  upstream  and  downstream  edges 
of  the  cavity. 
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(b)  Variation  of  maximum  pressure  coefficient  with  cavity  length  L/82  [Ref.  49]. 

Figure  9.  Concluded. 

In  their  study  of  coaxial  flow  over  a  combination  of  a  disk  and  cylinder 
separated  by  a  gap,  Koenig  and  Roshko  [Ref.  60]  defined  two  flow  regimes  based 
on  the  magnitude  of  the  drag  of  the  combined  body.  These  regimes  were 
identified  as  low  and  medium  drag.  The  medium  drag  is  significantly  lower  than 
the  drag  that  exists  when  the  bodies  are  well  separated  and  the  downstream  body 
has  little  influence  on  the  upstream  body,  but  in  the  low-drag  regime  drag  is  about 
an  order  of  magnitude  lower.  Their  study  suggests  that  the  low-drag  flows  appear 
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to  be  more  general  examples  of  the  classical  cavity  flow.  From  the  Gharib  and 
Roshko  [Ref.  49]  work,  it  is  now  clear  that  the  low-drag  regime  of  Ref.  60 
corresponds  to  the  regime  for  self-sustained  cavity  flow  oscillations  while  the 
medium-drag  regime  corresponds  to  the  cavity  wake  mode.  The  flow 
visualizations  and  flowfield  measurements  of  Ref.  49  show  that  in  the  low-drag 
regime  the  cavity  shear  layer  always  stagnates  at  the  downstream  comer.  Only  in 
the  wake  mode,  the  increased  drag  case,  does  it  stagnate  inside  the  downstream 
comer.  The  stagnation  location  moves  further  into  the  cavity  as  the  cavity  length 
increases. 

In  summary,  the  mechanism  for  a  stable  equilibrium  of  the  shear  layer  may 
be  explained.  An  inward  displacement  of  the  shear  layer  at  the  rear  comer  lowers 
the  feedback  signal  and  the  resultant  Reynolds  stress.  This  decreases  the 
entrainment  in  the  initial  portion  of  the  shear  layer  and  reduces  the  positive 
outflow.  Correspondingly,  the  inflow  at  the  downstream  end  must  be  reduced  and 
the  shear  layer  returns  to  its  initial  position.  In  the  self-sustained  flow  oscillation 
modes  the  flow  is  very  stable  and  resistant  to  moderate  external  disturbances  until 
the  gap  becomes  so  large  that  the  wake-mode  instability  takes  over.  These 
phenomena  have  been  shown  to  determine  that  a  cavity  in  the  self-sustained  flow 
oscillation  modes  has  low  drag  and  that  a  cavity  in  the  wake  mode  has  high  drag. 


D.  NUMERICAL  PREDICTIONS 


The  earliest  separated-eddy  and  cavity-flow  computations  were  reviewed  in 
detail  by  Burggraf  [Ref.  61].  These  computations  ranged  from  incompressible, 
inviscid  such  as  Rossow  [Ref.  62]  to  solutions  of  full  Navier-Stokes  equations  at 
low  Reynolds  numbers.  The  Burggraf  paper  included  a  numerical  solution  to  the 
full  Navier-Stokes  equations  using  a  stream  function  and  vorticity  formulation,  for 
the  case  of  a  square  cavity  in  the  Reynolds  number  (UL/o)  range  from  0  to  400. 
From  all  of  the  studies  reviewed,  a  fairly  clear  description  was  obtained  for  the 
overall  flow  characteristics.  However,  the  detailed  structure  of  secondary  vortices 
was  poorly  represented,  even  for  these  low  Reynolds  number  flows.  Pan  and 
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Acrivos  [Ref.  63]  calculated  steady  flow  in  rectangular  cavities  (L/D  from  0.25  to 
5)  where  the  flow  was  driven  by  a  uniform  translation  of  the  top  wall  of  the  cavity. 
These  creeping  flow  solutions  were  obtained  using  a  relaxation  technique  in  a 
formulation  similar  to  that  of  Burggraf.  Computed  streamline  contours  of  the 
secondary  comer  vortices  were  an  improvement  upon  earlier  solutions.  However, 
a  minimum  grid  size  of  0.0 1L  was  too  coarse  to  represent  detailed  streamline 
patterns  inside  the  comer  eddies.  It  should  be  noted  that  these  comer  vortices 
occupied  only  about  0.5  percent  of  the  total  area  of  the  cavity  cross-sectional  area. 
Unfortunately,  the  numerical  method  had  serious  instabilities  for  Reynolds 
numbers  greater  than  400,  which  made  those  results  invalid. 

The  flow  in  a  two-dimensional  channel  with  a  rectangular  cavity  was 
studied  numerically  by  Metha  and  Lavan  [Ref.  64],  For  simplification,  the  length 
of  the  channel  was  taken  to  be  infinite  and  the  upper  wall  of  the  channel  was 
moved  with  a  constant  velocity.  This  kept  the  flow  approaching  the  cavity 
identical  in  all  cases.  Steady  laminar,  incompressible  flow  in  two-dimensional 
channels  with  a  rectangular  cutout  cavity  were  computed  using  an  explicit 
numerical  solution  of  the  Navier-Stokes  equations  in  a  stream  function  and 
vorticity  formulation.  Solutions  were  obtained  for  cavity  L/D  of  0.5,  1 .0,  and  2.0 
and  for  Reynolds  numbers  of  1,  10,  100,  and  500.  These  cavity  L/D  ratios  were 
chosen  to  give  reattachment  of  the  flow  over  the  cavity  and  not  inside.  One  vortex 
was  observed  for  cavities  with  L/D  of  1  and  2.  For  a  L/D  of  0.5  there  were  two 
vortices  present,  one  on  top  of  the  other.  The  streamline  dividing  the  external  flow 
and  the  cavity  flow  was  concave  for  very  low  Reynolds  numbers  and  convex  for 
the  higher  Reynolds  numbers.  As  the  Reynolds  was  increased,  the  vortex  center 
moved  downstream  and  upward,  creating  a  thin  shear  layer.  The  shear  layer  on 
top  of  the  cavity  and  along  the  cavity  wall  was  not  very  thin,  even  at  the  larger 
Reynolds  number.  This  result  suggested  that  a  relatively  coarse  grid  may  be 
adequate  to  define  the  flow  features. 

Borland  [Ref.  65]  obtained  numerical  solutions  for  the  oscillating  flowfield 
in  an  open  cavity  exposed  to  a  high  subsonic  freestream  for  comparison  with  both 
wind  tunnel  and  flight  test  data.  Two-dimensional  Euler  equations  for  time- 
dependent  inviscid  compressible  flow  were  formulated  and  solved  using  two 
different  numerical  algorithms.  A  first-order,  fluid-in-cell  method  (FLIC)  which 
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used  an  entirely  Eulerian  scheme  provided  the  fundamental-frequency  pressure- 
oscillation  mode  results.  The  grid  cells  were  control  volumes  where  the  equations 
for  conservation  of  mass,  momentum,  and  energy  were  solved  subject  to 
appropriate  boundary  conditions.  A  second-order,  fluid-in-cell  MacCormack 
predictor-corrector  method  was  also  used  in  an  attempt  to  predict  some  of  the 
higher-frequency  oscillation  modes.  These  finite-difference  equations  were  used 
for  cell  boundaries  which  were  within  the  computational  region.  Upstream  or 
inflow  boundary  conditions  (b.c.'s)  used  the  freestream  variables.  Downstream  or 
outflow  b.c.'s  assumed  that  the  gradients  of  the  basic  flow  variables  vanish.  This 
is  a  more  complex  b.c.’s  and  sometimes  produces  less  reliable  results.  The  solid 
wall  b.c.'s  used  cell  boundaries  which  were  coincident  with  the  physical 
boundaries.  For  inviscid  flow,  reflective  conditions  were  used  where  the  normal 
velocities  at  the  solid  walls  are  cancelled  by  assuming  a  fictitious  grid  cell  with  an 
opposing  velocity  just  beyond  the  physical  boundary.  All  of  the  other  primitive 
flow  variables  are  also  duplicated  by  the  fictitious  grid  cell.  The  first-order  FLIC 
results  gave  a  good  correlation  with  the  fundamental  oscillation  frequency  and  the 
magnitude  of  the  pressure  oscillations  observed  experimentally.  The  second-order 
predictor-corrector  method  showed  a  capability  to  predict  higher-order  harmonics 
as  well  as  the  fundamental.  It  was  found  that  this  computation  was  not  stable  and 
the  oscillatory  flow  field  damped  out  with  increasing  numbers  of  iterations  and 
was  not  observed  late  in  the  computation.  The  outflow  b.c.  was  suspected  to  be 
the  cause  of  this  computational  instability. 

Hankey  and  Shang  [Ref.  66]  analyzed  pressure  oscillations  in  an  open 
cavity  using  time-dependent  numerical  computations  of  supersonic  flow  over  an 
open  cavity  with  a  L/D  equal  to  2.25.  The  unsteady  Navier-Stokes  equations  were 
solved  by  the  MacCormack  finite- difference  explicit  method  for  a  freestream 
Mach  number  of  1.50,  and  a  Reynolds  number  of  2.6  x  107.  The  most  obvious 
feature  was  that  the  flowfield  within  the  cavity  is  subsonic,  except  in  the  vicinity 
of  the  cavity  opening.  The  orderly  development  of  the  shear  layer  above  the 
cavity  was  shown.  An  attached  turbulent  boundary  layer  upstream  of  the  cavity 
separates  at  the  front  wall  to  form  the  shear  layer  over  the  cavity.  The  shear  layer 
then  reattaches  downstream  of  the  cavity.  The  pressure  oscillation  was  also 
predicted,  and  both  the  fundamental  frequency  and  magnitude  were  in  agreement 
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with  experimental  data  [Ref.  51].  However,  the  numerical  solution  would  be 
needed  for  a  longer  period  to  completely  verify  the  experimental  spectral  analysis. 
This  paper  provided  the  first  complete  viscous  solution  of  the  pressure  oscillating 
cavity. 

Ghia,  Ghia,  and  Shin  [Ref.  67]  used  a  multigrid  method  to  solve  the  Navier-Stokes 
equations  for  incompressible  flow  in  a  cavity.  The  laminar  incompressible  flow 
was  computed  in  a  square  cavity  whose  top  wall  moved  with  a  uniform  velocity  in 
its  own  plane  for  moderately  high  values  of  Reynolds  number  (100  to  10,000). 
The  objective  of  the  paper  was  to  achieve  the  solutions  in  a  greatly  reduced 
computational  time.  The  two-dimensional  cavity  flow  was  represented 
mathematically  by  a  stream  function  and  vorticity  form  of  the  Navier-Stokes 
equations.  The  non-slip  boundary  conditions  were  applied  at  the  cavity  walls  by 
requiring  zero  normal  velocities  at  all  of  the  boundaries.  The  two-dimensional 
cavity  space  was  discretized  by  a  uniform  mesh.  Second-order  accurate  central 
finite-difference  approximations  were  employed  for  all  of  the  second-order 
derivatives.  The  convective  terms  were  represented  by  a  first-order  upwind 
difference  scheme.  In  the  multigrid  method,  the  mesh  density  is  cycled  from  low 
to  high,  iteratively,  to  smooth  the  results  by  damping  the  high  frequency  errors. 
The  solutions  used  the  coupled,  strongly-implicit  procedure  of  Rubin  and  Khosla 
[Ref.  68], 

An  example  of  the  computational  results  from  Reference  67  is  presented  in 
Figure  33  which  shows  the  streamline  contours  in  Figure  33  for  a  Reynolds 
number  of  10,000  using  a  uniform  grid  (257  x  257).  The  center  of  the  primary 
vortex  was  offset  towards  the  top  right  corner  at  a  Reynolds  number  of  100  (x  = 
0.6172  and  y  =  0.7344).  It  moves  toward  the  geometric  center  of  the  cavity  with 
increase  of  the  Reynolds  number.  Above  a  Reynolds  number  of  5000  the 
movement  ceases  and  the  primary  vortex  was  located  at  x  =  0.5117  and  y  = 
0.5233.  As  the  Reynolds  number  increased  from  100  to  10,000  the  vorticity 
strength  at  the  vortex  center  decreased  from  3.17  to  1.88. 

This  solution  shows  the  presence  of  additional  counter-rotating  vortices  in 

or  near  the  cavity  comers.  The  notation  shown  in  Figure  10  uses  the  letters  T,  B, 
L,  and  R  to  denote  top,  bottom,  left,  and  right,  respectively.  For  example,  BR2 
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refers  to  the  second  in  the  sequence  of  secondary  vortices  that  occur  in  the  bottom 
right  comer  of  the  cavity.  At  a  Reynolds  number  of  100  there  were  only  two 
secondary  vortices;  as  the  Reynolds  number  was  incresed  to  10,000  the  number  of 
secondary  vortices  increased  to  the  six  shown  in  Figure  10.  Initially,  all  of  the 
secondary  vortices  are  located  very  near  the  cavity  corners.  With  increasing 
Reynolds  number  the  secondary  vortices  move,  very  slowly,  toward  the  cavity 
center.  Figure  33  includes  magnified  views  of  the  contours  of  the  secondary 
vortices.  The  values  of  the  stream  function  for  each  contour  on  the  figure  are 
listed  below: 


contour  letter 

¥ 

contour  number 

¥ 

a 

-1  x  10-10 

0 

X 

o 

1 

00 

b 

-1  x  10-7 

1 

1  X  10-7 

c 

-1  x  10-5 

2 

1  x  10-6 

d 

-1  x  10-4 

3 

1  x  10-5 

e 

-0.0100 

4 

5  x  10-5 

f 

-0.0300 

5 

1  x  10-4 

g 

-0.0500 

6 

2.5  x  10-4 

h 

-0.0700 

7 

5  x  10-4 

i 

-0.0900 

8 

1  x  10-3 

j 

-0.1000 

9 

1.5x10-3 

k 

-0.1100 

10 

3  x  10-3 

1 

-0.1150 

m 

-0.1175 

33 


Figure  10.  Streamline  pattern  for  primary,  secondary,  and  additional  comer 
vortices  for  a  Reynolds  number  of  10,000  using  an  uniform  grid  (257  x  257). 


In  the  last  ten  years  there  have  been  many  additional  computational 
solutions  of  the  cavity  flow  problem  [Ref.  69  through  78].  These  efforts  covered  a 
wide  range  of  flow  conditions.  In  spite  of  this  large  amount  of  computational 
activity  there  are  still  some  major  deficiencies.  The  extensive  regime  of  unsteady 
and  oscillatory  flow  over  and  inside  cavities  still  contains  many  unresolved 
problems.  Prediction  of  the  possible  modes,  amplitudes,  and  frequencies  have 
been  very  difficult  to  obtain  accurately.  Many  of  these  prediction  efforts  have  not 
adequately  modeled  the  effect  of  the  initial  boundary  layer  thickness  on  the  shear 
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layer  interaction  with  the  cavity  flow  oscillations.  Attempts  to  get  around  the 
details  of  the  shear  layer  have  usually  succeeded  only  in  situations  where  the  shear 
layer  is  not  the  controlling  factor,  such  as  for  deep  cavities  and  large  enclosures 
with  small  openings.  Improvements  upon  the  these  deficiencies  will  require  a 
significant  development  effort. 


E.  MULTIPLE  CAVITIES 


A  two-dimensional  inviscid,  incompressible  computational  investigation  of 
the  aerodynamics  of  an  airfoil  with  a  vortex  trapped  by  two  spanwise  fences  was 
conducted  by  Rosso w  [Ref.  62].  To  evaluate  the  concept,  computed  results  were 
obtained  for  a  Clark-Y  airfoil  using  a  variety  of  fence  geometries  at  several  angles- 
of-attack  to  determine  the  aerodynamic  characteristics  of  these  trapped- vortex, 
high-lift  devices.  These  computational  results  suggested  that  two  spanwise  fences 
should  be  used  to  enclose  a  trapped  vortex.  The  two  fence  heights  were  adjusted 
so  that  an  equilibrium  condition  could  more  easily  be  achieved  than  using  only  a 
single  fence.  Little  or  no  mass  removal  from  the  core  region  of  the  vortex  was 
required  to  stabilize  the  trapped  vortex.  The  lift  increase  was  found  to  be  inversely 
proportional  to  the  chordwise  spacing  of  the  fences.  It  appeared  that  for  two 
fences  there  would  be  no  profile  drag  penalty,  that  the  vortex  would  be  easier  to 
form  and  would  be  more  stable  than  for  a  single  fence.  The  results  also  showed 
that  the  vortex  bubble  could  be  moved  fore  and  aft  on  the  airfoil  to  control  the 
pitching  moment.  It  was  suggested  that  an  extension  of  the  single  trapped-vortex 
geometry  to  two  or  more  trapped  vortices  on  the  upper  surface  of  an  airfoil  may 
provide  the  same  lift  with  less  cumbersome  equipment. 

The  progression  in  complexity  from  the  two-dimensional  configurations 
studied  above  to  the  full  three-dimensional  flowfield  of  a  wing  requires 
examination  of  several  additional  factors.  Such  devices  were  recently  studied  in 
an  unpublished  experimental  investigation  by  Rossow  and  Ross  [Ref.  79]  using  a 
low-aspect-ratio,  highly-swept-back  (A  =  60°)  wing  (Figure  11)  with  a  semispan 
of  24.76  inches  and  a  chord  of  18.98  inches.  There  were  differences  between  the 
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actual  wing  area  and  the  resultant  aspect  ratio  and  the  values  used  in  the 
aerodynamic  coefficients  as  the  reference  geometry  as  shown  below: 


actual  geometry 
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S,  in^ 

783.90 

467.66 

AR 

3.13 

5.24 

This  test  was  conducted  in  the  NASA  Ames  Research  Center  7-  by  10-  Foot 
Tunnel  using  a  semispan  model.  The  model  was  mounted  on  a  metric  balance 
frame.  Force  data  were  obtained  from  an  external  floating-frame  system  of  beam 
balances  located  below  the  test-section  floor. 


Figure  11.  Semispan  wing  (NACA  0012  airfoil  normal  to  the  leading  edge)  with 
two  part-span  spoilers  tested  by  Rossow  and  Ross  [Ref.  79]. 

A  variety  of  spoiler  configurations  were  tested.  As  a  sample  of  results  from 
this  investigation,  the  following  data  were  obtained  using  a  pair  of  spoilers  (0.08c 
high  at  0.15c  and  0.10c  high  at  0.45c)  which  were  used  to  form  a  single  cavity 
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swept  along  the  local  chord  line.  The  spoilers  extended  from  the  root  chord  to 
approximately  1/3  of  the  wing  span.  The  intent  was  to  take  advantage  of  the 
inherent  spanwise  flow  to  form  a  standing  vortex  and  to  avoid  active  mass  removal 
from  the  core  region  of  the  vortex.  If  the  spoilers  are  effective,  then  the  flow  over 
the  upper  surface  would  be  displaced  passively  to  increase  the  wing  camber.  This 
experimental  investigation  found  that  below  8°  angle-of-attack  there  was  little  or 
no  lift  increase.  However,  there  was  a  lift  benefit  as  high  as  30%  at  an  angle-of- 
attack  greater  than  10°  (Figure  12(a)). 


(a)  Lift  curve. 

Figure  12.  Rossow-Ross  experimental  data  [Ref.  79]. 

However,  there  was  an  unexpected  and  undesired  drag  increase  over  a  wider 
angle-of-attack  range  as  shown  by  the  polar  plot  in  Figure  12(b).  For  reference, 
the  ideal  polar  (l/7t  AR  e)  is  shown  by  a  solid  line,  where  the  span  efficiency 
factor  e  is  1,  which  corresponds  to  an  elliptic  spanload  distribution.  An  efficient 
transport  wing  will  typically  have  a  span  efficiency  factor  greater  than  0.9. 
Lower-aspect-ratio  fighter  aircraft  have  a  span  efficiency  which  drops  into  the  0.7 
or  0.8  range.  Even  lower  span-efficiency  factors  indicate  either  an  inefficient 
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spanload  distribution  or  the  presence  of  flow  separation.The  polar  for  the  wing 
without  spoilers  shows  a  drag  increase  above  the  ideal  drag  polar  for  Cp  >  0.3  (or 

a  >  5°).  The  drag  for  the  spoiler-deployed  configuration  is  increased  even  further. 
In  an  attempt  to  quantify  the  drag  increase  the  experimentally  determined  span 
efficiency  factor,  e,  is  tabulated  in  the  following  table: 


(b)  Lift-drag  polar. 
Figure  12.  Concluded. 


CL 

spoilers  deployed 

none 

two 

0.4 

0.57 

0.31 

0.6 

0.56 

0.33 

0.8 

0.58 

0.38 

1.0 

0.42 
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These  adverse  drag  results  suggest  the  presence  of  flow  separation  and  indicate  the 
need  for  a  better  understanding  of  the  flow  phenomena  involved. 

In  the  early  1960's,  Migay  [Ref.  80]  conducted  a  series  of  investigations 
which  evaluated  the  effect  of  multiple  transverse  cavities  in  a  diffuser  (Figure 
13(a)).  The  diffuser  surface  had  a  constant  radius  of  11  centimeters.  The  first 
cavity  was  located  at  the  entrance  to  the  diffuser  (  a  =  0°).  Up  to  12  transverse 
cavities  distributed  along  the  flow  direction  through  the  diffuser  (at  a  >  0°)  were 
tested.  Pressure  taps  were  distributed  along  the  diffuser  in  the  direction  of  the 
flow.  The  pressures  were  non-dimensionalized  using  a  pressure  coefficient 
referenced  to  atmospheric  pressure 


c 


_  Pi  Patm 


(12) 


where  U|a=0  is  the  velocity  in  the  test  section  upstream  of  the  diffuser.  The 

pressure  data  were  presented  as  the  ratio  of  the  pressure  coefficient  measured  at  a 
particular  tap  location  to  the  pressure  coefficient  measured  at  a  =  0.  The  measured 
pressure  distribution  is  presented  in  Figure  13  (b)  and  it  is  seen  that  the  largest 


curved  diffuser  with 


Figure  13.  Results  of  Migay  investigation  [Ref.  80]. 
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pressure  reduction  occured  near  a  =  10°.  It  appeared  that  the  surface  flow 
separated  downstream  of  the  a  =  10°  location.  It  was  found  that  the  cavities  were 
effective  in  preventing  flow  separation  as  long  as  the  first  cavity  was  located  at  an 
oc  between  0°  and  10°.  When  the  first  cavity  was  located  further  downstream,  the 
cavities  had  no  beneficial  effect  on  the  location  of  flow  separation.  The  effect  of 
the  number  of  cavities  is  presented  in  Figure  13  (c)  where  the  ratio  of  pressure 
coefficients  is  shown  as  a  function  of  the  number  of  cavities.  No  effect  was  found 
for  up  to  three  cavities.  Then  there  was  an  increasing  pressure  coefficient 
reduction  from  four  to  eight  cavities  with  no  additional  change  for  more  than  eight 
cavities.  These  results  indicate  that  the  cavity  benefit  can  be  achieved  with  a 
limited  number  of  cavities  if  they  are  properly  located  in  the  diffuser. 

An  investigation  by  Lin  et  al  [Ref.  81-86]  evaluated  in  some  detail  the 
potential  of  multiple  cavities  as  a  base-drag  reduction  device.  An  example  of  these 
results  is  presented  in  Figure  14  for  an  axisymmetric  body  with  a  modified  bluff 
base  which  had  up  to  13  circumferential  grooves.  The  transition  from  the 
axisymmetric  body  to  the  tapered  base  used  a  radius  on  the  shoulder  at  the 
intersection  of  the  body  and  the  base.  For  the  results  in  Figure  14,  the  shoulder 
radius  was  one-half  of  the  body  radius.  The  downstream  edge  of  the  first  groove 
was  located  at  the  start  of  the  shoulder  radius.  The  grooves  were  numbered  from  1 
to  13  as  their  location  varied  from  the  shoulder  to  more  downstream  positions. 
The  bluff-body  base  had  a  30°  angle;  a  streamline  body  with  a  base  angle  of  10° 
was  also  tested.  Both  bodies  had  a  diameter  of  2.39  inches;  the  bluff  body  was  3 
diameters  long  while  the  streamline  body  was  nearly  5  diameters  long. 

The  data  in  Figure  14  present  the  variation  of  drag  coefficient  with 
Reynolds  number  for  axisymmetric-body  configurations  with  and  without 
circumferential  grooves.  They  show  that  the  streamline  body  had  the  lowest  drag 
and  that  the  ungrooved  bluff  body  had  one  of  the  highest  drag  coefficients.  The 
effect  of  the  grooves  was  determined  by  opening  the  grooves  one  at  a  time. 
Groove  number  1  was  opened  at  the  shoulder  location.  Then  groove  number  2 
was  also  opened  at  the  next  downstream  location.  This  process  was  continued 
until  all  13  grooves  were  open.  When  one  groove  was  opened  the  drag  increased 
above  that  for  the  ungrooved  body.  When  a  second  groove  was  opened  there  was 
a  drag  reduction.  However,  for  1,  2,  or  3  grooves  open,  the  drag  was  greater  than 
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for  the  ungrooved  body.  As  the  number  of  grooves  was  increased  from  4  to  10  the 
drag  was  increasingly  reduced  below  that  for  the  ungrooved  body. 

The  largest  drag  reduction  was  achieved  with  either  10  or  13  grooves  open. 
For  these  two  configurations  the  lowest  drag  occurred  in  the  Reynolds  number 
range  from  about  80,000  to  100,000  and  then  drag  tended  to  increase  slowly  with 
further  Reynolds  number  increase.  It  was  found  that  the  boundary  layer 
transitioned  on  the  base  from  laminar  to  turbulent  in  this  Reynolds  number  range. 
As  a  result,  these  circumferential  grooves  on  the  body  with  a  shoulder  radius 
appear  to  only  be  effective  for  a  laminar  boundary  layer  approaching  the  shoulder. 
The  reason  for  this  Reynolds-number  sensitivity  was  not  identified.  This 
investigation  also  evaluated  many  other  devices  and  found  several  which  were  as 
good  or  better  than  multiple  cavities  without  being  limited  to  a  particular 
Reynolds-number  range.  This  study  also  indicated  a  need  to  improve  the 
understanding  of  the  flow  changes  due  to  the  cavities. 


GROOVE  NO.  1.  2.  3.  4.  ETC. 


O.M  0.08  0.12  0.16  0.20  *  106 


Figure  14.  Effect  of  circumferencial  grooves  on  the  afterbody  on  the 
axisymmetric  bluff-body  drag;  shoulder  radius  is  0.2  body  radius  [Ref.  86]. 
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III.  EXPERIMENTAL  PROCEDURE  AND  RESULTS 


The  present  investigation  included  a  wind  tunnel  test  where  the  effect  of 
single  and  multiple  cavities  on  the  boundary  layer  development  and  the  pressure 
distributions  were  measured  in  both  a  zero  and  two  adverse  pressure  gradients. 
The  boundary  layer  profiles  were  measured  at  four  or  five  longitudinal  stations  in 
the  vicinity  of  the  cavity(s).  One  station  was  located  immediately  upstream  of  the 
first  cavity  and  the  rest  were  either  between  the  multiple  cavities  and/or 
downstream  of  the  cavities. 


A.  WIND  TUNNEL  AND  APPARATUS 


The  experimental  investigation  was  conducted  in  the  San  Jose  State 
University  (SJSU)  Low  Speed  Wind  Tunnel.  The  test  apparatus  included  a  three- 
dimensional  traverse  from  the  Naval  Postgraduate  School  as  well  as  a  computer 
and  data  aquisition  system  from  the  NASA  Ames  Research  Center.  The  test 
hardware  for  the  cavity  configurations  and  the  test  section  diffusers  were 
specifically  fabricated  for  this  investigation  in  the  SJSU  tunnel. 


1.  Wind  Tunnel 


A  sketch  of  the  SJSU  tunnel  is  presented  in  Figure  15.  This  is  a  continuous 
flow  tunnel  with  an  Eiffel-type  return.  The  flow  enters  the  30"  square  entrance 
section,  shown  on  the  left  in  Figure  16,  and  then  goes  through  a  honeycomb  and 
two  screens  into  the  36"  long  contraction.  The  contraction  ratio  is  6.25.  The  test 
section  is  12"  high,  12"  wide,  and  24"  long.  The  contraction,  test  section,  and 
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diffuser  are  mounted  on  a  fixed  frame.  The  drive  motor-fan  assembly  is  located 
on  a  movable  frame  located  on  the  right  side  of  the  sketch.  There  is  also  a 
honeycomb  upstream  of  the  motor.  Between  the  downstream  end  of  the  diffuser 
and  the  honeycomb  there  is  a  variable  gap  which  is  used  to  control  the  tunnel 
speed.  As  the  gap  was  varied  from  6"  to  0"  the  test  section  speed  Uoo  increased 
from  40  feet  per  second  (fps)  (q  =  1.9  psf)  to  the  maximum  tunnel  speed,  142  fps 
(q  =  24  psf).  The  present  test  was  conducted  using  a  speed  of  130  feet  per  second. 
A  96"  long  test  section  extension  was  constructed  for  this  investigation  to  increase 
the  boundary  layer  thickness  for  a  portion  of  the  tests.  Longitudinal  position  in  the 
test  section  was  measured  from  the  test  section  entrance:  station  0"  through  24"  for 
the  24"  test  section;  and  stations  0"  through  120"  for  the  120"  test  section.  In  the 
longer  test  section  the  original  24"  long  test  section  was  located  between  the  96" 
and  120"  stations.  Photographs  of  the  SJSU  tunnel  with  each  test  section  length 
are  presented  in  Figure  17. 

To  achieve  adverse  pressure  gradients  in  the  24"  long  test  section,  the  first 
6"  was  reduced  in  height.  Then  the  increased  contraction  ratio  was  faired  using  an 
insert  on  the  tunnel  bottom  in  the  16"  immediately  upstream  of  the  test  section 
(x<0").  A  test  section  diffuser  was  installed  between  6"  and  24"  stations  on  the 
floor  of  the  test  section.  One  diffuser  had  a  constant  7°  angle  and  the  second  test 
section  diffuser  had  a  constant  14°  angle.  The  test  section  extension  height  was 
variable  to  accomodate  the  reduced  height  upstream  of  the  same  test  section 
diffusers  used  in  the  24"  test  section.  The  contraction  ratios  for  the  test  section 
with  the  7°  and  14°  test  section  diffusers  were  9.38  and  10.13,  respectively. 
Coordinates  are  presented  in  Figure  36  for  the  original  test  section  (solid  lines)  and 
for  the  two  test  section  diffusers  (dashed  lines).  The  x  and  y  scales  exagerate  the  y 
dimensions  relative  to  the  x  dimensions.  Tests  were  conducted  in  each  of  the  three 
configurations  shown;  (1)  original  test  section,  y  =  6"  (top)  and  -6"(  bottom);  (2) 
7°  test  section  diffuser,  y  =  6"  (top)  and  varying  from  -4"  to  -6"  (bottom);  and  14° 
test  section  diffuser,  y  =  6"  (top)  and  varying  from  -1.4"  to  -6"  (bottom).  The 
same  test  section  top  was  used  for  each  of  the  basic  tunnel  configurations.  For  the 
120"  test  section  the  same  test  section  diffusers  were  used.  As  mentioned  in  the 
previous  paragraph  their  x  locations  were  increased  by  96",  the  length  of  the  test 
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(b)  Extended  120"  long  test  section  configuration. 
Figure  17.  Concluded 
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(a)  Coordinates  of  the  24  "  test  section  and  for  the  7°  and  14°  test  section  diffusers. 

Figure  18.  Test  section  diffusers  used  to  generate  an  adverse  dp/dx. 

section  extension  which  was  added  to  the  front  of  the  original  test  section  and 
downstream  of  the  modified  contraction  ratio  (x  <  0"). 

Photographs  of  the  SJSU  tunnel  are  also  presented  in  Figure  18  which 
shows  the  7°  test  section  diffuser  with  four  one  inch  cavities  (Figure  18(b))  and  the 
14°  test  section  diffuser  with  one  two  inch  cavity  (Figure  18(c)).  The  test  section 
diffusers  were  fabricated  using  plywood  frames  mounted  in  the  streamwise 
direction  and  covered  by  0.25"  thick  clear  acrylic  sheets.  The  photograph  for  the 
7°  test  section  diffuser  also  shows  the  four  cavity  configuration  mounted  on  the 
top  of  the  test  section.  The  photograph  for  the  14°  test  section  diffuser  also  shows 
the  single  two  inch  cavity  configuration  mounted  on  the  top  of  the  test  section. 
Some  of  the  clutter  in  the  background  of  the  photographs  is  the  tygon  tubing  used 
for  the  pressure  distribution  measurements  on  the  tunnel  top  and  bottom.  The 
vertical  portion  of  the  pitot-static  probe  used  for  tunnel  reference  dynamic  pressure 
measurement  is  visible  above  the  circular  access  window.  The  boundary  layer 
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(b)  Photograph  of  the  7°  test  section  diffuser  with  flow  from  the  right. 

Figure  18.  Continued. 


(c)  Photograph  of  the  14°  test  section  diffuser  with  flow  from  the  right. 

Figure  18.  Concluded 
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total  pressure  probe  is  near  the  downstream  cavity  in  Figure  18(b)  or  near  the  two 
inch  cavity  in  Figure  18(c)  adjacent  to  the  tunnel  ceiling  in  both  photographs. 


2.  Test  Hardware 


The  test  hardware  built  assembled  specifically  for  the  present  investigation 
consisted  of  the  test  section  diffuser,  the  cavity  configurations  (1,  2,  and  4 
cavities),  the  three-dimensional  traverse  mechanism,  and  the  computer  with  its 
data  acquisition  program. 


a.  Single  and  Multiple  Cavity  Hardware 


There  were  four  cavity  configurations  built  for  the  investigation  of 
the  effect  of  the  number  of  cavities  on  the  boundary  layer  profiles.  As  shown  by 
the  photographs  in  Figure  18  the  cavities  were  mounted  in  the  test  section  ceiling 
in  the  adverse  pressure  gradient  portion  of  the  test  section  diffuser.  For 
nomenclature  purposes,  the  edge  of  the  cavity  adjacent  to  the  freestream  is  called 
the  top  and  the  cavity  horizontal  surface  parallel  to  the  freestream  is  called  the 
floor.  For  each  cavity  the  upstream  vertical  surface  is  called  the  front  wall  and  the 
downstream  vertical  surface  is  called  the  rear  wall.  Each  cavity  had  a  nominal 
length  to  depth  ratio  L/D  ofl.  Each  of  these  cavity  configurations  was  equipped 
with  a  distribution  of  surface  pressure  taps  both  in  the  cavities  and  on  the  ceiling 
surface  both  upstream  and  downstream  of  the  cavities.  The  coordinates  of  all  of 
these  pressure  taps  are  presented  in  Appendix  B,  Pressure  Tap  Locations.  The 
cavity  test  hardware  which  was  fabricated  for  this  investigation  is  summarized  in 
the  following  table: 
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number  of 

cavities 

cavity  number 
(from 
upstream) 

cavity  depth, 
D,  in. 

cavity  length, 
L,  in. 

longitudinal 
location  of 
cavity  front 
wall,  in. 

1 

1 

2.00 

1.94 

7.50 

1 

1 

0.94 

1.00 

6.50 

2 

1 

1.00 

0.94 

6.00 

2 

2 

1.00 

0.94 

8.00 

4 

1 

0.94 

0.94 

6.00 

4 

2 

0.99 

0.94 

7.95 

4 

3 

1.00 

0.94 

9.94 

4 

4 

0.99 

0.95 

11.95 

b.  Traverse 


The  three-dimensional  traverse  assembly  from  the  Naval 
Postgraduate  School  is  shown  in  Figure  19.  It  was  used  to  support  and  traverse  a 
pitot  probe  used  to  measure  the  boundary  layer  profile.  A  0.375  inch  circular  hole 
in  the  ceiling  provided  clearance  for  the  0.25  inch  diameter  probe  support.  The 
probe  was  used  to  survey  the  boundary  layer  profile  adjacent  to  the  test  section 
ceiling.  The  traverse  was  mounted  on  top  of  a  frame  which  wrapped  around  the 
tunnel  test  section  and  did  not  touch  the  wind  tunnel.  Three  electric  motors 
powered  the  gear-driven  screws  which  moved  the  traverse  independently  in  each 
Cartesian  direction.  The  following  traverse  ranges  were  available:  streamwise. 
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Figure  19.  Photograph  of  the  traverse  mechanism  mounted  at  the  SJSU  tunnel. 


20";  laterally,  22.5";  and  vertically,  12".  These  movements  were  executed  by 
manually  entering  the  desired  movement  direction  (X,  Y,  or  Z)  and  the  movement 
distance.  The  repeatibility  of  probe  positioning  was  determined  by  moving  the 
probe  ten  inches  and  then  moving  back  to  the  starting  point.  There  was  less  than 
0.002  inch  error  in  returning  to  the  original  position.  Significantly  no  backlash 
due  to  reversing  the  movement  direction  was  observed.  To  assure  the  best 
positioning  accuracy,  all  of  the  boundary  layer  traverses  were  obtained  by  moving 
in  one  direction  from  the  ceiling  into  the  freestream  to  avoid  any  possibility  of 
backlash.  For  the  next  run,  the  probe  was  then  returned  to  a  position  hard  against 
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the  test  section  ceiling.  The  direction  was  then  reversed  to  a  position  where  the 
probe  just  touched  the  ceiling  prior  to  the  next  traverse. 


c.  Boundary  Layer  Probe 


The  boundary  layer  total  pressure  was  measured  using  a  total 
pressure  probe  which  was  made  by  flattening  a  circular  stainless  steel  tube.  Using 
a  0.060  inch  diameter  tube,  a  probe  tip  height  of  0.022  inch  was  used  for  runs  101 
through  254.  This  probe  was  damaged  and  replaced  by  another  probe  with  a  tip 
height  of  0.034  inch  for  runs  255  through  380.  The  probe  location  for  plotting 
boundary  layer  profiles  was  taken  to  be  one-half  of  the  probe  tip  height  when  it 
touched  the  wall.  For  the  tests,  this  half  height  was  identified  as  zero  on  the  probe 
traverse  position  data  display. 

The  pressure  was  measured  using  a  0.5  psid  transducer.  The  NASA 
Ames  calibration  laboratory  found  that  this  particular  transducer  had  a  standard 
deviation  equal  to  0.141%  of  full  scale.  This  translated  to  a  0.5%  (0.102  psf)  error 
for  differential  pressures  measured  at  a  freestream  dynamic  pressure  of  20  psf. 
After  the  probe  was  installed,  the  calibration  was  checked  using  an  electronic 
pressure  calibration  device.  Regression  analyzes  of  repeated  calibrations  gave  a 
correlation  coefficient  r^  of  between  0.9986  and  0.9997.  This  represented  a 
maximum  repeatibility  error  of  0.07%  of  full  scale  or  0.05  psf.  A  third  error 
source  was  the  electronic  zero  shift  between  the  beginning  and  end  of  a  run.  In 
154  runs  this  shift  averaged  0.03  psf.  If  these  errors  were  summed  the  maximum 
error  becomes  0.18  psf  or  0.9%  of  the  freestream  dynamic  pressure. 

A  final  error  source  was  the  effect  of  the  wall  proximity  or  boundary 
layer  displacement  effect.  The  work  of  Coles  [Ref.  9]  was  applied  to  the  present 
data.  It  was  found  that  this  error  was  negligible  when  used  to  correct  the 
displacement  thickness  Si  and  the  momentum  thickness  §2-  Application  of  this 
correction  to  the  measured  velocity  adjacent  to  the  wall  indicated  a  velocity  error 


55 


AV/V  of  0.0145.  At  a  freestream  dynamic  pressure  of  20  psf,  this  was  an  error  of 
0.004  psf  or  0.02%.  This  was  considered  to  be  negligible. 


d.  Computer  and  Data  Acquisition  System 


The  digital  data  acquisition  system  was  controlled  by  a  desktop 
computer  through  a  graphical  user  interface  (GUI)  computer  program  and  two 
dedicated  data  acquisition  circuit  boards.  Each  data  channel  had  its  own  amplifier 
which  conditioned  the  data  signal  and  amplified  its  voltage  so  that  it  could  be 
sampled  with  the  greatest  possible  accuracy  and  then  multiplexed  to  the  computer. 
The  computer  had  8  megabytes  (Mb)  of  random  access  memory  (RAM),  an  40  Mb 
internal  hard  disk,  and  a  removable  44  Mb  hard  disk  which  was  used  to  store  the 
experimental  data. 

The  application  program  for  the  GUI  was  a  set  of  virtual  instruments 
(VI)  which  each  consistd  of  three  parts:  (1)  the  front  panel;  (2)  the  block  diagram; 
and  (3)  the  icon/connector.  The  first  two  parts  were  visible  on  the  computer 
screen  as  windows.  The  front  panel  was  the  user  interface  to  the  VI  and  was  used 
during  data  acquisition  by  the  program  operator  to  interact  with  it.  It  included 
input  objects  called  controls  and  output  objects  called  indicators.  The  objects  were 
portrayed  on  the  computer  screen  as  command  buttons  and  status  indicators.  Data 
was  also  be  displayed  in  either  tabular  and/or  grapical  forms.  The  block  diagram 
was  a  VI  source  code  which  was  created  using  a  graphical  programming  language. 
This  was  a  virtual  wiring  diagram  for  connecting  the  program  modules  which  were 
called  nodes  to  the  terminals  for  the  controls  and  indicators.  The  icon  was  a 
representation  of  the  VI.  The  connector  actd  as  the  port  through  which  the  data 
passed.  The  icon/connector  together  represented  the  VI  in  a  manner  analogous  to 
a  subroutine  call  statement  when  the  VI  was  used  as  a  sub VI  in  another  Vi's  block 
diagram.  Many  layers  of  sub  Vis  could  be  used  in  higher  order  Vis. 

The  present  investigation  GUI  used  a  main  panel  (Figure  20(a)),  a 
boundary  layer  measurements  panel  (Figure  20(b)),  and  a  pressure  distribution 
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measurements  panel  (Figure  20(c)).  The  data  acquisition/reduction  program  was 
executed  by  first  loading  the  GUI  program  application  code  and  then  launching  the 
cavity  flow  experiment  main  panel  VI  and  all  of  the  related  sub  Vis.  At  the 
completion  of  these  steps  the  computer  screen  displayed  the  main  panel  shown  in 
Figure  20(a).  There  were  nine  command  buttons  shown  on  the  top  portion  of  this 
panel  and  status  indicators  for  each  data  channel.  The  data  acquisition  was 
initiated  by  sequentially  using  the  computer's  mouse  to  select  the  clear  acquisition 
button,  then  select  the  initialize  acquisiton  buttom,  and  finally  select  the  take  zeros 
button.  Then  there  were  five  choices  available:  (1)  data  acquisition  setup  button 
which  was  used  to  add  or  remove  data  channels;  (2)  on-line  monitoring  button  to 
view  data  output  for  each  channel;  (3)  boundary  layer  measurements  button;  (4) 
pressure  distribution  measurements  button;  and  (5)  quit  button  to  terminate  and 
exit  the  program. 

In  this  investigation  four  groups  of  data  measurements  were  taken 
and  recorded  in  data  files:  (a)  tunnel  conditions;  (b)  boundary  layer  profiles;  and 
(c)  static  pressure  distributions.  An  example  of  the  data  recorded  for  a  boundary 
layer  profile  run  is  presented  in  Figure  21.  This  particular  run  was  identified  as 
R0105BL  which  stands  for  run  R,  number  0105  and  boundary  layer  profile  BL, 
respectively.  Subsequent  boundary  layer  profile  runs  have  different  run  numbers. 
As  shown  in  Figure  21,  the  data  file  generated  by  the  GUI  program  started  with  the 
run  number,  the  date  and  time  of  the  run,  and  the  tunnel  code  (see  Appendix  A  for 
a  detailed  description).  Then  there  was  an  additional  line  of  group  1  data  which 
consisted  of  the  freestream  conditions  measured  by  a  pitot  probe  located  at  x  =  1" 
and  y  =  2.5"  and  then  averaged  over  the  entire  run.  These  data  include  the  tunnel 
reference  values  for  pitot  probe  static  Ps,tun  and  total  pressures  Pt,tun  ,  as  well 
as  their  difference,  the  calculated  dynamic  pressure  qc.  For  an  independent 
measurement,  the  dynamic  pressure  was  also  measured  by  a  separate  pressure 
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(a)  Cavity  flow  experimental  main  panel. 

Figure  20.  Front  panels  used  as  the  GUI  in  the  present  investigation. 
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(b)  Boundary  layer  instruments  panel 
Figure  20.  Continued. 
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(c)  Pressure  distribution  measurement  panel. 
Figure  20.  Concluded 
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transducer  and  displayed  as  q.  For  pressure  coefficients  this  value  of  dynamic 
pressure  was  used  to  non-dimensionalize  the  static  pressures.  Finally  the  tunnel 
total  temperature  Tt,tun  and  atmospheric  pressure  P  atm  were  listed.  The  tunnel 
conditions  were  always  measured  along  with  one  of  the  other  three  groups  of  data 
described  above. 

Shown  in  Figure  21  is  a  group  of  boundary  layer  probe  data, 
identified  as  group  2  data.  The  number  of  data  points  varied  from  run  to  run.  In 
this  example  from  run  105,  15  points  were  recorded.  Then  seven  columns  of  data 
measurements  were  displayed  for  each  boundary  layer  height:  (1)  the  height  from 
the  bottom  of  the  boundary  layer  total  pressure  probe  to  the  tunnel  ceiling  (y  = 
6.00")  Y  bl;  (2)  the  boundary  layer  total  pressure  Pt,p;  (3)  the  ratio  of  local  u 
velocity  to  estimated  local  freestream  velocity  u/U;  (4)  pitot  probe  dynamic 
pressure  q,tun;  (5)  tunnel  total  temperature  Tt,tun;  (6)  pitot  probe  total  pressure 
Pt,tun;  and  (7)  pitot  probe  static  pressure  Ps,p.  The  height  of  the  boundary  layer 
measurement  was  the  sum  of  the  probe  height  Y  bl  plus  one-half  of  the  probe 
vertical  thickness  Ayp/2  which  was  0.011  inch  for  runs  100-257  and  0.017  inch 

for  runs  258-380.  The  four  measured  pressures  and  the  measured  temperature 
were  tabulated  in  the  engineering  units  shown  in  Figure  21.  The  boundary  layer 
profile  in  run  105  was  taken  at  x  =  9",  which  is  downstream  of  the  location  of  the 
single  cavity  configurations.  The  boundary  layer  profile  is  expressed  as  the  ratio 

of  boundary  layer  velocity  computed  from  the  total  pressure  measured  by  the 
boundary  layer  probe  pt?p  at  the  height  Y  bl  +  Ayp/2  divided  by  the  freestream 

velocity  at  the  x  station  where  the  profile  was  measured.  The  velocity  within  the 
boundary  layer  at  a  given  height  was  computed  from  the  Bernoulli  equation  as 


where  the  static  pressure  measured  at  a  static  pressure  orifice  located  on  the  ceiling 
at  the  x  station  of  the  boundary  layer  profile  was  psurf.  The  quantity  Pt,p-Psurf 

was  measured  as  the  pressure  difference  across  a  pressure  transducer. 

During  the  on-line  measurements  the  local  freestream  velocity  was 
estimated  by  the  GUI  program.  For  tunnel  configurations  without  the  diffuser 
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inserts,  the  local  freestream  velocity  was  assumed  to  be  constant  along  the  length 
of  the  test  section.  For  tunnel  configurations  with  a  diffuser  the  local  velocity  U 
or  Ux  was  computed  using  the  incompressible  continuity  equation 

U„  =  U„^  (14) 

Ax 

where  either;  (1)  Ax  =  12*y@x=5  for  or  for  Ax  where  x  <  6  or;  (2) 
Ax  =12(y@x=5  +  (x-6)sin(a))  for  Ax  when  x  >  6.  The  test  section  diffuser  had 

a  constant  height  from  x  =  0"  to  x  =  6".  From  x  =  6"  to  x  =  24"  the  diffuser  had  a 
constant  slope  of  either  7°  or  14°.  For  the  final  processed  data,  Ux  was  replaced 
by  the  velocity  at  the  edge  of  the  boundary  layer  Ue  which  was  computed  from 
the  boundary  layer  probe  measurements  at  y  >  d  where  the  maximum  value  of  pt  ^ 
-  pSUIf  was  measured. 

An  example  of  group  3  data  from  run  173  is  presented  in  Figure  22. 
This  particular  run  was  identified  as  R0173CP  which  stands  for  run  R,  run 
number  0173  and  pressure  coefficient  data  CP,  respectively.  The  run  heading, 
run-date-time-tunnel  code,  and  group  1  data  are  displayed  the  same  way  they  were 
shown  in  Figure  21.  The  group  3  data,  pressure  distributions,  are  grouped  by  the 
scanivalve  module  that  measured  a  group  of  pressure  taps.  Each  section  of  these 
data  is  headed  by  the  scanivalve  number  (1,  2,  or  3)  and  the  number  of  pressures 
recorded.  Notice  that  the  number  of  pressures  is  one  greater  than  the  number  of 
taps  for  each  scanivalve.  The  extra  pressure  is  a  reference  pressure  which 
represents  zero  return  after  the  group  of  taps  were  sampled  and  indicates  whether 
there  is  any  electrical  drift  in  the  pressure  data  during  the  data  sample  period. 

Scanivalve  1  (Figure  22(a))  measured  the  16  taps  located  on  the  test 
section  diffuser  between  x  ~  2"  and  x  ~  21".  When  there  was  no  test  section 
diffuser  installed,  scanivalve  1  was  not  sampled;  only  the  test  section  ceiling 
pressures  were  measured.  For  run  173  scanivalve  2  (Figure  22(b))  measured  the 
17  taps  located  on  the  test  section  ceiling  upstream  (from  x  ~  2")  and  downstream 
of  the  two  inch  cavity  (to  x  «  22").  Notice  that  pressure  tap  13  measured  a 
pressure  nearly  equal  to  zero  for  this  run.  Since  this  was  not  seen  in  measurements 
of  this  tap  in  other  runs,  this  measurement  was  interpreted  as  faulty  data  for  this 
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run  and  omitted  from  the  final  data  plots.  Scanivalve  3  (Figure  22(c))  measured 
the  21  taps  located  in  and  near  the  two  inch  cavity.  As  seen  in  the  data  listing,  the 
upstream  cavity  wall  was  located  at  x  =  7"  and  the  downstream  wall  was  located  at 
x  =  9".  The  floor  of  the  cavity  was  located  at  z  =  8"  (for  the  2"  cavity).  The  data 
listing  for  each  scanivalve  includes  the  following  information:  (1)  tap  number;  (2) 
tap  x  location;  (3)  tap  y  location;  (4)  the  surface  static  pressure  Psurf;  (5)  the 
pressure  coefficient  Cp  which  is  defined  as 

c  Psurf -Ps.tun  (15) 

q 

The  complete  list  of  runs  completed  in  this  investigation  is  presented  in  Appendix 
A.  The  surface  tap  locations  for  all  of  the  test  configurations  are  listed  in 
Appendix  B.  The  cavity  pressure  data  are  presented  in  Appendix  C  for  all  of  the 
measured  data.  The  boundary  layer  profile  data  are  presented  in  Appendix  D  for 
all  of  the  measured  data. 
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(a)  Group  1  tunnel  conditions  data  and  group  3  scanivalve  (SV)  1  pressure  data. 
Figure  22.  Data  file  generated  by  the  pressure  distribution  measurement  panel. 
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(b)  Group  3  scanivalve  (SV)  2  pressure  data. 
Figure  22.  Continued. 
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(c)  Group  3  scanivalve  (SV)  3  pressure  data 
Figure  22.  Concluded. 
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B.  WIND  TUNNEL  RESULTS 


The  experimental  wind  tunnel  data  described  in  the  previous  section  and 
presented  in  the  appendices  are  analyzed  in  this  section  in  three  parts:  (1) 
boundary  layer  profile  analyses;  (2)  static  surface  and  cavity  pressure  data 
analyses;  and  (3)  dynamic  flow  due  to  cavity  opening.  The  entrance  velocity 
profile  was  surveyed  in  the  24  inch  test  section  at  the  x  =  4"  station  using  the 
boundary  layer  total  pressure  probe.  It  was  found  that  the  flow  dynamic  pressure 
was  uniform  (less  than  0.2  percent  variation)  with  an  increase  of  about  0.8  percent 
in  dynamic  pressure  at  the  edge  of  the  ceiling  and  floor  boundary  layers.  A  single 
hot-wire  probe  was  located  at  the  tunnel  centerline  of  the  24"  test  section  at  the  x  = 
4"  station  to  measure  the  streamwise  flow  turbulence.  These  data  found  that  the 
Au/u  was  approximately  0.003. 


1 .  Test  Section  Boundary  Layer 


For  the  plain  test  section  ceiling,  the  boundary  layers  were  measured  at  four 
longitudinal  locations  in  the  original  24"  test  section  (x  =  5.5",  9.0",  12.0",  and 
16")  and  in  the  extended  120"  test  section  (x  =  101.5",  105.0",  108.0",  and  112.0"). 
The  test-section  configurations  included  the  basic  12"  by  12"  cross  section  test 
section  with  essentially  zero  pressure  gradient  (dp/dx  «  0)  and  the  two  test  section 
diffuser  configurations  with  adverse  pressure  gradients  (dp/dx  >  0).  These  data  are 
analyzed  in  the  following  section  to  evaluate  how  representative  the  present  basic 
tunnel  data  are  of  a  two-dimensional  boundary  layer  over  a  flat  plate. 


a.  Zero  Longitudinal  Pressure  Gradient 

For  a  two-dimensional,  turbulent  boundary  layer  in  a  zero  pressure 
gradient,  the  velocity  profile  may  be  estimated  using  the  1/7  power  law 
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u/U^  =  (y/S)(1/7).  The  1/7  power  law  may  then  be  integrated  [Ref.  4]  to  obtain  the 
above  boundary  layer  thickness  quantities  as  8/x  =  0.37  R^1/5) 
81/x  =  0.37  R~^  /  (1  +  n)  where  n  =  7,  and  82/x  =  0.036  Rx-^^.  In  the 
present  investigation  these  relations  only  apply  to  the  cases  where  the  velocity  is 
constant  with  increasing  run  length  x  ;  that  is  for  zero  pressure  gradient,  the  tunnel 
configurations  without  the  7°  and  14°  diffusers  in  the  test  section. 

These  three  boundary  layer  thickness  parameters  (8,  8j,  and  82) 

were  computed  using  both  the  1/7  power  law  relations  and  the  experimental  data 
from  the  SJSU  tunnel  for  both  the  basic  24"  long  test  section  (at  x  =  5.5",  9",  12", 
and  16")  and  the  120"  long  test  section  (at  x  =  101.5",  105",  108",  and  112").  The 
results  are  compared  in  Figure  23.  In  the  24"  test  section  data  all  three 
experimental  thicknesses  were  found  to  be  slightly  greater  than  the  1/7  power  law 
thicknesses.  This  could  be  due  to  the  fact  that  x  was  measured  from  the  start  of  the 
test  section.  Ahead  of  the  test  section  in  the  contraction  the  pressure  was 
decreasing  (dp/dx  <  0)  so  that  the  boundary  layer  thicknesses  were  decreasing. 
However,  they  almost  certainly  do  not  go  to  zero  at  the  x  =  0  station.  If  the  initial 
station  were  moved  upstream  3  or  4  inches  to  a  virtual  origin  for  application  of  the 
1/7  power  law  (to  more  adequately  represent  the  equivalent  zero  thickness 
location)  then  the  experimental  data  would  agree  with  the  1/7  power  law  relations. 
The  experimental  boundary  layer  profiles  for  all  four  stations  where  the  boundary 
layer  was  measured  in  the  24"  test  section  agree  with  the  1/7  power  law  profile. 
As  the  Reynolds  number  increases  the  boundary  layer  profiles  are  more  nearly 
one-eighth  or  one-ninth  or  a  lower  power  profile.  In  the  120"  test  section  the 

boundary  layer  thicknesses  are  well  below  those  obtained  from  the  1/7  power  law. 
If  a  one-eleventh  power  law  is  used  to  compute  the  momentum  thickness  8j 

agreement  with  the  experimental  data  can  be  achieved.  For  the  thicker  boundary 
layer  in  the  120"  test  section,  the  presence  of  four  walls  may  be  the  reason  that  the 
boundary  layer  thickness  does  not  grow  as  large  as  a  two-dimensional  boundary 
layer.  While  all  of  the  reasons  for  the  difference  between  the  24"  and  120"  test 
sections  have  not  been  identified,  these  effects,  especially  for  the  longer  test 
section,  will  be  evaluated  numerically  in  a  later  section  using  the  incompressible 
Navier-Stokes  equations. 
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(a)  Boundary-layer  thickness,  8. 

Figure  23.  Comparison  of  experimental  data  from  present  investigation  using  the 
24"  and  120"  test  sections  with  the  1/7  power-law  boundary-layer  thicknesses. 
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(b)  Boundary-layer-displacement  thickness,  8j. 
Figure  23.  Continued. 
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(c)  Boundary-layer-momentum  thickness,  82- 
Figure  23.  Concluded. 

Rotta's  analysis  [Ref.  9]  introduced  universal  similarity  parameters 
which  were  later  demonstrated  by  Clauser  [Ref.  8].  Clauser  plotted  data  from 
three  separate  investigations  as  u/U  versus  y/5  which  appeared  to  lack  similarity. 
He  then  replotted  these  data  using  a  universal  similarity-parameter  plot  ((u-U)/u* 
versus  y/8).  This  change  greatly  reduced  the  data  scatter  and  convincingly 
demonstrated  similarity  independent  of  both  Reynolds  number  and  roughness.  A 
similar  correlation  was  done  using  the  data  from  the  present  investigation  where 
the  24"  test  section  data  (Figure  Dl)  deviated  from  the  1/7  power  law  and  the  120" 
test  section  data  (Figure  D2)  tended  to  agree  with  the  1/7  power  law.  The 
universal  similarity  parameter  form  of  these  data  is  presented  in  Figure  24(a).  The 
data  provide  a  good  correlation  for  most  of  the  boundary-layer  profile  with  only  a 
limited  dispersion  at  heights  below  0.15  for  the  shortest  run  length  data  (x  =  5.5 
inches).  This  correlation  provides  confidence  that  these  boundary  layer  profiles 
are  similar  and  representative  of  those  measured  in  a  two-dimensional  boundary- 
layer. 
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(a)  Experimental  data  from  present  investigation. 
Figure  24.  Universal  velocity  distribution  law  comparisons. 


In  Schlichting  [Ref.  4]  the  velocity  defect  laws  of  both  von  Karman 


and  Prandtl 


u  k  vy/5  J 


(16) 


(17) 


compared  favorably  with  experimental  data  from  both  smooth  and  rough  pipes.  It 
is  stated  that  both  laws  were  obtained  for  two-dimensional  flow  in  a  channel. 
They  were  found  to  agree  with  both  two-dimensional  and  axisymmetric 
experimental  data.  Both  of  these  equations  are  compared  in  Figure  24(b)  with  zero 
pressure  gradient,  flat  plate  experimental  data  both  from  the  present  investigation 
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(Figure  24(a))  and  from  the  experimental  data  correlation  done  by  Clauser  [Ref. 
8].  The  two  sets  of  experimental  data  correlate.  In  contrast  to  Ref.  4,  there  is  a 
difference  between  the  two  velocity  defect  laws  (equations  16  and  17)  and  the  two 
sets  of  experimental  data  in  Figure  24  (b). 
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(b)  Prandtl  and  von  Karman  equations  compared  with  present  investigation  and 
with  expermental  data  correlated  by  Clauser  [Ref.  17]. 

Figure  24.  Concluded. 


relation 


The  data  from  the  SJSU  low  speed  tunnel  is  compared  with  log-law 


(18) 


in  Figure  25  (a)  for  the  standard  24"  long  test  section  (where  6  =  0.23  at  x  =  5.5" 
and  8  =  0.45"  at  x  =  16")  and  in  Figure  25  (b)  for  the  extended  120"  long  test 
section  (8  ~  1.17"  at  the  locations  where  the  measurements  were  taken).  The 
comparison  is  typical  of  the  expected  agreement  in  the  log-law  region.  There  is  a 
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boundary  layer  probe  displacement  effect  at  y+  for  low  heights  in  the  inner  region; 
and  there  is  a  wake  effect  at  y+  for  the  greater  heights.  The  wake  is  characterized 
by  a  region  at  heights  approaching  the  edge  of  the  boundary  layer  region  where  the 
measured  velocity  is  greater  than  that  given  by  the  log  law.  The  specific  locations 
for  these  effects  are  given  below: 


test  section  length,  inches 

24 

120 

boundary  layer  probe  displacement  effect  below  y+ 

100 

100 

boundary  layer  wake  effect  at  y+  greater  than 

300 

1000 

In  Figure  25  (a)  for  the  24"  test  section  there  is  a  larger  wake  effect 
at  the  x  =  12"  location  than  at  the  other  locations.  It  should  be  noted  that  the 
boundary  layer  run  lengths  are  relatively  short  and  so  it  should  not  be  surprising 
that  there  are  differences  where  the  boundary  layer  is  in  its  initial  development.  In 
contrast,  as  shown  in  Figure  25  (b)  for  the  120"  test  section,  where  there  is  a  fully 
developed  turbulent  boundary  layer,  there  are  no  significant  differences  among  the 
wakes.  As  a  consequence,  for  the  y+  regions  both  below  the  boundary  layer 
displacement  effect  limit  and  above  the  wake  effect  limit  the  experimental  U+ 
data  is  greater  than  the  log  law  relation.  In  between  these  y+  limits  where 

equation  13  applies,  the  data  in  the  log  law  region  agrees  well  with  the  above  log- 
law  equation  when  the  constants  are  k  =  0.41  and  Cj  =  5.0. 

Computation  of  the  turbulent  boundary  layer  was  the  subject  of  the  1968 
AFSOR-IFP-Stanford  Conference  on  Computation  of  Turbulent  Boundary  Layers 
[Ref.  9].  Volume  II  of  reference  9  included  a  compilation  of  experimental 
boundary  layer  data  by  Coles  and  Hirst  and  a  paper  by  Coles  which  was  a  guide  to 
the  experimental  data  compiled  for  use  at  the  conference  by  the  various  predictive 
methods  which  were  being  evaluated.  In  this  paper,  Coles  discussed  two  problems 
associated  with  evaluation  of  experimental  data:  (1)  the  boundary  layer 
displacement  effect  where,  near  the  wall,  the  data  are  in  error  because  of  the 
effects  of  probe  wall  interference  and  local  strong  turbulence;  and  (2)  the  wake 
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effect  near  the  edge  of  the  boundary  layer  where  the  data  are  correct  but  the  law- 
of-the-wake  equations  are  not  correct  as  discussed  earlier  in  this  section.  To 
account  for  the  near  wall  effects,  Coles  used  a  standard  sublayer  profile  to  obtain 
the  following  standard  functions  to  reproduce  the  tabulated  thicknesses: 

Jo5°U+dy+ =540.6  and  Jo50(U+)2dy+ =6546,  (19) 

The  experimental  boundary  layer  displacement  thickness  and  the  momentum 
thickness  have  been  recomputed  using  these  relations  (equations  19)  and  are 
shown  in  Figures  26(a)  and  26(b),  respectively.  Here  the  1/7  power  law  results  are 
compared  with  the  experimental  data  from  the  SJSU  tunnel  for  both  the  integrated 

uncorrected  data  and  these  data  modified  by  the  above  correction  procedure.  It  is 
shown  that  the  corrections  made  only  a  small  change  to  and  82-  Since  these 

changes  are  so  small,  for  all  of  the  rest  of  the  boundary  layer  data  in  the  present 
paper  only  the  integrated  experimental  data  without  the  Coles  correction  will  be 
used. 


x,  in. 


(a)  Boundary-layer  displacement-thickness. 

Figure  26.  Comparison  of  1/7  power  law  boundary-layer  thicknesses  with 
experimental  data  both  with  and  without  Coles  near-wall  correction. 
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(b)  Boundary-layer  momentum-thickness. 

Figure  26.  Concluded. 

Relations  for  turbulent  wall-friction  were  determined  using  the 
experimental  velocity  profile  data.  Experimentally,  the  skin  friction  coefficient 
may  be  found  using  a  variety  of  skin-friction  strain  gauge-devices,  by  an  obstacle 
block,  by  a  heat  transfer  analogy  [Ref.  16],  by  surface  flow  visualization,  by  a 
razor  blade  technique,  by  a  sublayer  fence  [Ref.  14  through  17],  by  a  Preston  tube 
[Ref.  18],  by  empirical  correlation  [Ref.  17  and  19]  or  by  analysis  of  the  boundary 
layer  profile  [Ref.  8,  17,  and  20].  In  the  present  report  the  skin  friction  was 
determined  from  the  boundary  layer  profile  using  a  Clauser  plot  [Ref.  8].  The  log- 
law  equation  (equation  8)  may  be  re-arranged  to  become  an  implicit  equation  for 
skin  friction  coefficient,  giving 


+  Ci 


(20) 


An  example  of  the  application  of  this  equation  to  the  analysis  of  the  data  from  the 
present  investigation  is  presented  in  Figure  27  for  station  112  inches  in  the  120 
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•  Figure  27.  Application  of  the  Clauser  plot  for  the  determination  of  skin-friction 
coefficient  for  data  from  SJSU  120"  test  section  at  station  1 12  inches. 

inch  long  test  section  where  the  skin-friction  coefficient  was  found  to  be  0.0030. 
Typically,  repeat  runs  had  differences  in  cf  no  larger  than  0.0002.  These  results 
show  that  while  there  is  a  small  reduction  in  skin  friction  due  to  the  cavity,  the 
differences  from  the  plain  wall  configuration  were  nearly  within  measurement 
accuracy.  This  procedure  was  used  to  obtain  the  skin-friction  coefficient  from  all 
of  the  data  described  in  the  test  run  schedule  (Appendix  A). 

There  are  a  number  of  correlations  which  may  be  used  to  estimate  the  skin- 
friction  coefficient  on  a  flat  plate.  The  equation  based  on  the  1/7  power  law  has 
already  been  given  as  Cf  =0.0592(RX)-0'2  .  Another  correlation  often  used  is 
the  empirically  developed  skin-friction  equation  by  Schultz-Grunow  [Ref.  19] 

cf  =  0.37(logRx)-2'584  (21) 

The  skin  friction  coefficients  obtained  from  the  Clauser  plots  for  the  conditions 
with  zero  pressure  gradient  are  compared  in  Figure  28  (a)  with  the  corresponding 
values  obtained  from  the  1/7  power  law  and  with  the  Schultz-Grunow  equation. 
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The  correlations  and  the  experimental  data  agree  for  the  24"  test  section.  But  for 
the  120"  test  section  both  Cf  correlations  are  about  0.0005  less  than  the 

experimental  data.  Whether  this  difference  may  be  due  to  surface  roughness  in  the 
96"  test  section  extension  will  be  evaluated  in  the  fourth  subsection  of  the  general 
background  section. 

The  following  correlation  by  Ludwieg-Tillmann  [Ref.  17]  was  developed  to 
account  for  the  effect  of  pressure  gradients 

cf  =  0. 246(R52  )"0'268io-°-678H  (22) 

This  correlation  is  compared  with  data  from  both  test  section  lengths  in  Figure  28 
(b).  In  contrast  to  the  1/7  power  law  and  Schultz-Grunow  correlations  in  Figure 
28  (a),  the  Ludwieg-Tillmann  correlation  demonstrates  good  agreement  with  the 
Clauser  method  results.  There  is  scatter  for  the  data  from  the  24"  test  section. 
There  is  agreement  with  the  experimental  data  from  the  120"  test  section.  This 
latter  comparison  for  x  >  100"  represents  the  major  difference  between  the 


(a)  Comparison  of  Clauser  analysis  of  present  data  from  both  the  14"  and  120" 
test  sections  with  correlations  from  1/7  power  law  and  by  Schultz-Grunow. 

Figure  28.  Comparisons  of  skin  friction  from  several  methods. 
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(b)  Comparison  of  Clauser  analysis  with  the  Ludwieg-Tillmann  correlation. 

Figure  28.  Concluded. 

Ludwieg-Tillmann  correlation  and  both  the  1/7  power  law  and  Schultz-Grunow 
correlation. 


Initially,  the  pressure  distributions  for  the  24"  and  120"  long  test 
sections  were  measured.  An  example  of  the  measured  pressure  distribution  on  the 

short  test  section  ceiling  is  presented  in  Figure  29  using  open  circle  symbols. 
These  data  show  a  reduction  in  Cp  with  increasing  x  distance  from  the  test  section 

dCp 

entrance.  The  measured  pressure  gradient  was  — -  =  -0.000643.  This  cD 

dx  F 

decrease  was  caused  by  the  small  increase  in  local  velocity  due  to  the  expected 

increase  in  displacement  thickness.  At  an  arbitrary  x  location  the  effective  test 

section  cross-section  area  for  the  flow  may  be  expressed  as  the  physical  dimension 

reduced  by  the  boundary  layer  displacement  thickness  on  each  surface,  giving 

A  =  (12.0-  2  Sj)2  (23) 

and  the  local  velocity  using  the  incompressible  continuity  equation  becomes 
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(24) 


U2  —  UjAj  /  A2 


Finally,  the  computed  Cp  is  given  by 


1- 


vU.y 


\2 


(25) 


The  Cp  change  due  to  the  increasing  is  called  the  buoyancy  correction.  For  the 

four  stations  where  the  boundary-layer  profile  was  measured,  the  buoyancy- 
induced  Cp  was  computed  and  is  presented  in  Figure  29  as  solid  square  symbols. 

The  two  curves  shown  provide  reasonable  agreement  and  illustrate  why  the 
pressure  gradient  was  negative. 
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Figure  29.  Effect  of  buoyancy  on  pressure  distribution  on  the  ceiling  of  the  24" 

long  test  section. 
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b.  Adverse  Pressure  Gradient 


The  pressure  distributions  measured  in  the  original  test  section  and 
in  the  test  section  with  the  diffuser  are  presented  in  Figure  30.  When  either  test 
section  diffuser  (7°  or  14°)  was  installed  the  reference  velocity  (or  dynamic 
pressure)  was  measured  by  a  pitot-static  probe  located  at  station  1".  The  surface 
pressure  data  from  both  test  section  lengths  showed  a  difference  between  the  test 
section  ceiling  flat  surface  and  the  test  section  floor  with  a  deflected  surface.  On 
the  test  section  floor  at  station  6"  (or  102"),  the  surface  had  an  obtuse  angle  where 
there  was  an  abrupt  change  of  flow  direction.  As  discussed  in  Milne-Thompson 
[Ref.  87]  for  potential  flow,  the  velocity  is  locally  infinite  at  an  angle  apex.  The 
effects  of  a  sharp  corner  in  viscous  flow  also  show  an  increased,  but  finite, 
velocity.  At  the  angle  apex,  measured  pressure  coefficients  were  of  the  order  of 
minus  one  and  increased  rapidly  a  small  distance  from  the  comer  to  a  few  tenths. 
These  flows  are  discussed  in  greater  detail  by  Mason  [Ref.  88].  The  data  from  the 
present  investigation  confirms  the  trends  discussed  above.  The  floor  pressure 
distribution  had  a  discontinuity  near  the  angle  apex  at  station  6"  (or  102")  where 
the  local  pressure  coefficient  peaked  as  low  as  -0.19  at  x  =  7.4"  for  the  24"  test 
section  and  at  about  -0.1  at  x  =  5"  for  the  120"  test  section.  These  pressure 
coefficients  correspond  to  10%  or  5%  velocity  increases,  respectively. 

The  inviscid  pressure  distribution  for  both  test  section  diffusers  was 
computed  using  the  low-order  panel  method  PMARC  [Ref.  89]  and  the  results  for 
the  7°  diffuser  are  presented  in  Figure  31.  Without  viscosity  the  pressure  peak 
near  the  bottom  apex  angle  is  seen  to  peak  at  about  -0.10  for  the  7°  diffuser  which 
is  similar  to  the  experimental  data  and  at  about  -1.40  for  the  14°  diffuser  which  is 
much  lower  than  the  experimental  data  in  Figure  3 1 .  This  indicates  that  while  the 
flow  for  the  7°  diffuser  was  attached,  the  flow  for  the  14°  diffuser  may  have 
separated  near  the  angle  apex  and  as  a  result  generate  a  reduced  pressure  peak. 
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Figure  31.  Test-section  pressure  distribution  computed  by  the  inviscid  panel- 
method  PMARC  [Ref.  90]  for  the  24"  test-section  with  the  7°  test-section  diffuser. 

Since  the  velocity  was  subsonic,  the  streamlines  above  the  angle 
apex  should  quickly  smooth  out,  and  the  disturbance  die  down  exponentially  with 
increasing  distance  [Ref.  88].  The  ceiling  pressure-distributions  for  all  six  test 
section  configurations  show  the  expected  and  desired  smooth  profile  in  the  vicinity 
of  the  cavity  location(s).  The  14°  diffuser  data  still  has  a  slope  change  near  station 
6"  (or  102")  and  shows  the  effect  on  the  ceiling  pressure-distribution  of  the  change 
in  pressure  gradient  slope  at  station  6"  (or  102").  The  comparison  between  the  two 

test  sections  shows  that  the  thicker  boundary  layer  (see  Figure  30)  in  the  120"  test 
section  was  more  effective  in  smoothing  the  Cp  distribution  on  the  ceiling 

pressure-distribution.  These  data  also  show  that  there  was  no  significant  flow 
separation  apparent  in  these  experimentally  measured  profiles. 

The  experimental  pressure  distributions  in  the  0°,  7°,  and  14° 

diffuser  constant-angle  portions  of  both  the  24"  and  120"  long  test  sections  are 
presented  in  Figure  32.  These  data  show  a  linear  variation  of  Cp  with  x  location. 

Repeatibility  is  shown  by  the  data  for  the  14°  diffuser. 
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x,  in. 

(a)  24"  long  test  section. 

Figure  32.  Effect  of  diffuser  angle  on  the  test  section  ceiling  pressure  distribution. 


x,  in. 

(b)  120"  long  test  section. 
Figure  32.  Concluded. 
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The  resultant  pressure  gradients  are  tabulated  below  and  plotted  in  Figure  33. 


Test  Section 
Length,  in. 

Diffuser  Angle, 
deg. 

height  at  test 
section  diffuser 
entrance(x=6"),  in. 

dcn 

p  ,  1/in. 
dx 

24. 

0 

12.0 

-0.00064 

24. 

7 

10.0 

0.01103 

24. 

14 

7.4 

0.02443 

120 

0 

12.0 

-0.00036 

120 

7 

10.0 

0.00944 

120 

14 

7.4 

0.01653 

0  5  10  15 

diffuser  angle 


Figure  33.  Test  section  pressure  gradient  as  a  function  of  the  diffuser  angle. 
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Also  shown  in  Figure  33  is  the  pressure  coefficient  for  inviscid  flow. 
This  curve  was  computed  using  equation  24  with  Uj  =  U^,  the  velocity  at  the  test 

section  entrance,  and  equation  25.  The  effect  of  boundary  layer  growth  is  to 
increase  local  velocity  above  that  computed  inviscidly,  as  described  above  for 
Figure  29.  The  slopes  shown  in  Figure  33  are  listed  below: 


Case 

(dx)Ai/deg- 

inviscid 

0.00297 

24"  long  test  section 

0.00175 

120"  long  test  section 

0.00121 

These  data  show  a  41  %  and  a  60%  reduction  for  viscous  effects  on  the  pressure 
gradient  per  degree  of  diffuser  angle  in  the  24"  and  120"  long  test  sections, 
respectively. 


Boundary-layer  profiles  were  measured  on  the  ceiling  of  the  test 
section  at  four  longitudinal  stations  in  each  test  section:  x  =  5.5",  9",  12",  and  16" 
in  the  24"  test  section;  or  x  =  101.5",  105",  108",  and  112"  in  the  120"  test  section. 
These  boundary-layer  profiles  are  presented  in  Appendix  D  in  the  following 
figures: 


Configuration 

Test  Section 
length,  in 

Diffuser  Angle, 
deg. 

Figure 

0000 

24 

0 

D1 

8000 

120 

7 

D2 

0700 

24 

14 

D3  | 

8700 

120 

0 

D4 

87 


Configuration 

Test  Section 
length,  in 

Diffuser  Angle, 
deg. 

Figure 

1400 

24 

7 

D5 

8400 

120 

14 

D6 

These  profiles  were  used  to  compute  the  boundary-layer  thickness  8,  displacement 
thickness  5  j,  momentum  thickness  82,  pressure-gradient  parameter  [3  ,  the  ratio 

of  Clauser  defect  thickness  to  boundary-layer  thickness  A/8,  and  shape  factor  G 
which  are  presented  in  the  following  figures: 


Parameter 

8 

52 

P 

A/8 

G 

Figure 

34 

35 

36 

37 

38 

39 

The  data  from  the  24"  test  section  are  presented  in  part  (a)  and  from  the  120"  test 
section  are  presented  in  part  (b)  of  each  of  the  above  figures. 

The  boundary-layer  thickness,  boundary-layer  displacement 
thickness,  and  the  boundary-layer  momentum  thickness  data  for  the  basic  test 
sections  (configurations  0000  and  8000)  without  pressure  gradients  were 
compared  in  Figure  23  with  results  from  application  of  the  1/7  power  law.  While 
the  24"  test  section  data  were  found  to  agree  with  the  1/7  power  law,  the  higher 
Reynolds  number  120"  test  section  data  tended  to  agree  with  a  lower  order  1/11 
power  law.  For  the  24"  test  section  the  addition  of  the  7°  and  14°  diffusers 
reduced  the  values  of  8,  8j,  and  82  at  corresponding  stations  in  the  divergent 

section  (x  =  9",  12",  and  16").  In  contrast,  for  the  120"  test  section  the  addition  of 
the  7°  and  14°  test-section  diffusers  increased  the  values  of  8,  8j,  and  82  at  the 

corresponding  stations  in  the  divergent  section  (x  =  105",  108",  and  112"). 

The  Clauser  equilibrium  parameter  [3  is  often  used  to  characterize 
the  magnitude  of  a  pressure  gradient.  The  dominant  term  in  this  parameter  is  the 
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pressure  gradient,  dp/dx.  This  term  is  multiplied  by  the  ratio  of  displacement 
thickness  to  skin  friction.  In  Clauser's  paper  [Ref.  8],  constant  values  of  (3  were 
used  to  characterize  equilibrium  boundary  layers.  The  results  in  Figure  37  show 
that  the  present  investigation  achieved  equilibrium  for  configurations  0000,  8000, 
and  0700.  The  Clauser  equilibrium  parameter  for  the  24"  test  section  presented  in 
Figure  37(a)  shows  a  variation  similar  to  the  actual  pressures  shown  in  Figure 
32(a).  The  Clauser  equilibrium  parameter  for  the  120"  test  section  presented  in 
Figure  37(b)  also  shows  a  variation  similar  to  the  actual  pressures  shown  in  Figure 
32(b)  except  for  the  7°  diffuser.  Inspection  of  the  values  that  make  up  the  Clauser 
equilibrium  parameter  show  that  the  skin  friction  coefficients  and  displacement 
thicknesses  were  nearly  constant,  which  lead  to  a  small  variation  of  (3. 

Another  measure  of  the  pressure  gradient  is  the  Clauser  shape  factor 
[Ref.  8],  which  is  the  ratio  of  the  second  to  the  first  moment  about  the  axis  (U- 
u)/u*  =  0.  Another  useful  parameter  is  the  Clauser  defect  thickness  A.  The 
relation  of  G  to  the  universal  parameter  A/6  was  originally  shown  by  Clauser 
[Ref.  8].  Equivalent  information  from  the  present  investigation  are  shown  in 
Figure  38.  These  data  indicate  that  the  test  section  diffusers  induced  a  relatively 
modest  adverse  pressure  gradient. 

The  conventional  shape  factor  H  is  often  used  as  an  indication  of 
flow  separation.  Schlichting  [Ref.  4,  page  630]  indicates  that  separation  occurs 
when  H  =  1.8  to  2.4.  The  data  in  Figure  39  show  that  all  of  the  test  sections  gave 
shape  factors  which  were  well  below  1.8.  It  is  interesting  to  observe  that  the  no- 
diffuser  configuration  for  the  24"  test  section  had  shape  factors  greater  than  either 
of  the  two  diffuser  configurations.  In  contrast,  the  opposite  relation  was  measured 
in  the  120"  test  section.  All  of  these  data  show  that  the  goal  of  obtaining  adverse 
pressure  gradients  without  flow  separation  was  achieved  in  the  present 
investigation. 

The  log-law  relation  (equation  18)  for  the  data  from  the  SJSU  low 
speed  tunnel  was  presented  in  Figure  25  and  showed  the  expected  agreement.  The 
corresponding  data  for  the  7°  and  14°  test  section  diffusers  are  presented  in  Figure 
40.  The  data  for  the  24"  test  section  (Figures  40(a)  and  (b))  show  reasonable 
agreement  with  the  law-of-the  wall  equation  using  the  generally-accepted 
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empirical  constants.  In  contrast,  the  120"  test-section  data  fell  below  the  law-of- 
the  wall  equation.  This  is  consistent  with  the  fact  that  the  boundary  layer 
thickness  8  was  less  than  that  given  by  the  1/7  power  law.  It  appears  that  this 
boundary  layer  thickness  is  characteristic  of  a  higher  Reynolds  number  flow 
consistent  with  a  1/9  to  as  low  as  a  1/1 1  power-law  correlation  depending  on  the  x 
station  of  the  boundary-layer  survey. 

Clauser  [Ref.  8]  used  a  universal  plot  of  equilibrium  turbulent 
velocity  profiles,  (U-u)/u*  versus  y/A,  to  display  the  effect  of  boundary-layer  run 
length  on  the  boundary-layer  shape.  The  increasing  x  distance  was  shown  in 
Figure  37(a)  to  correspond  to  the  greatest  increase  in  the  pressure-gradient 
parameter  P  achieved  in  this  investigation.  The  resultant  effect  of  the  pressure 
gradient  on  the  velocity  profile  is  shown  in  Figure  41(b)  using  data  from  the  24" 
test  section  with  the  14°  diffuser.  The  changing  velocity  profile  from  station  5.5" 
to  station  16"  due  to  increasing  P  is  consistent  with  the  changes  shown  by 
Clauser. 


Since  these  profiles  seemed  to  depend  more  on  x  location  than 
specific  diffuser  geometry,  the  data  from  Figure  41  was  replotted  at  a  given  x 
location  for  the  24"  test  section  in  Figure  42  and  for  the  120"  test  section  in  Figure 
43.  The  data  at  station  5.5"  (Figure  42(a))  showed  an  anomaly  in  that  an 
apparently  large  effect  of  pressure  gradient  was  present.  For  this  station  the  7° 
data  was  apparently  invalid  and  the  14°  data  showed  a  large  profile  change.  The 
other  three  stations  from  the  24"  test  section  showed  only  small  profile  changes 
consistent  with  the  data  correlation  shown  in  Figure  38.  For  the  120"  test  section 
the  expected  larger  profile  changes  due  to  the  larger  pressure  gradient  parameter 
variation  are  seen.  At  x  =  101.5"  there  was  no  pressure  gradient,  so  the  boundary 
layer  profiles  are  seen  to  be  in  equilibrium.  That  is,  all  three  surveys  are  the  same. 
At  the  next  three  stations  there  was  an  increasing  gradient  in  the  pressure  variation 
with  increasing  x  which  is  consistent  with  increasing  pressure  gradient  parameter, 
P. 
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(a)  24"  test  section. 

Figure  34.  Effect  of  test-section  diffuser  angle  on  the  variation  of  boundary-layer 
thickness  as  a  function  of  longitudinal  location. 


(a)  24"  test  section. 

Figure  35.  Effect  of  test-section  diffuser  angle  on  the  variation  of  boundary-layer 
displacement  thickness  as  a  function  of  longitudinal  location. 


(a)  24"  test  section. 

Figure  36.  Effect  of  test-section  diffuser  angle  on  the  variation  of  boundary-layer 
momentum  thickness  as  a  function  of  longitudinal  location. 


(b)  120"  test  section. 
Figure  36.  Concluded. 
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(a)  24"  test  section. 

Figure  37.  Effect  of  test-section  diffuser  angle  on  the  variation  of  freestream 
pressure-gradient  parameter  as  a  function  of  longitudinal  location. 


(b)  120"  test  section. 
Figure  37.  Concluded. 
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(a)  24"  test  section. 

Figure  38.  Effect  of  test-section  diffuser  angle  on  the  ratio  of  the  Clauser  defect 
thickness  to  the  boundary-layer  thickness  as  a  function  of  the  Clauser  shape  factor. 
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(b)  120"  test  section. 
Figure  38.  Concluded. 
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(a)  24  test  section  with  7°  diffuser. 

Figure  40.  Effect  of  test-section  diffuser  angle  on  the  variation  of  boundary-layer 
thickness  as  a  function  of  longitudinal  location. 
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(b)  24  test  section  with  14°  diffuser. 
Figure  40.  Continued. 
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(c)  120"  test-section  with  7°  diffuser. 
Figure  41.  Continued. 
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(d)  120"  test-section  with  14°  diffuser. 
Figure  41.  Concluded. 
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(U-u)/u* 


(a)  Station  5.5"  in  the  24"  test  section. 

Figure  42.  Effect  of  test-section  diffuser  angle  on  the  universal  velocity  profiles. 


(U-u)/u* 


(b)  Station  9"  in  the  24"  test  section. 
Figure  42.  Continued. 
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(c)  Station  12"  in  the  24"  test  section. 

Figure  42.  Continued. 


(d)  Station  16"  in  the  24"  test  section. 
Figure  42.  Concluded. 
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(a)  Station  101.5"  in  the  120"  test  section. 

Figure  43.  Effect  of  test-section  diffuser  angle  on  the  universal  velocity  profiles. 


(b)  Station  105"  in  the  120"  test  section. 
Figure  43.  Continued. 


103 


10 


(U-u)/u* 


(U-u)/u* 


il 

leg. 

I 

--•c 

a 

-  < 

>-  14 

- 

■ 

L . 

■ 

n 

< ) 

0  0.2  0.4  0.6  0.8 

y/A 


(c)  Station  108"  in  the  120"  test  section. 
Figure  43.  Continued. 


(d)  Station  112"  in  the  120"  test  section. 
Figure  43.  Concluded. 
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2.  Open  Cavity  Configurations 


The  purpose  of  the  present  investigation  was  to  experimentally  and 
computationally  evaluate  the  effect  of  1,  2,  or  4  two-dimensional  spanwise  cavities 
with  square  cross-sections  on  the  boundary-layer  characteristics  in  both  a  constant 
pressure  flow  and  two  adverse  pressure-gradient  flows.  The  experimental  results 
presented  in  the  present  section  are  based  on  the  boundary-layer  profiles  presented 
in  Appendix  D  and  the  surface/cavity  pressure  distributions  presented  in  Appendix 
C.  These  data  will  be  compared  and  analyzed  in  the  following  two  subsections. 

a.  Effects  on  the  Boundary-Layer  Characteristics 


The  effects  of  the  cavity  geometry  on  the  boundary-layer 
characteristics  are  analyzed  in  this  subsection  to  identify  how  they  are  influenced 
by  the  flowfield.  The  boundary  layer  profiles  measured  on  the  tunnel  ceiling  in 
the  vicinity  of  the  cavities  were  used  to  compute  the  boundary  layer  thickness  8, 
displacement  thickness  8j  (equation  4),  momentum  thickness  82,  shape  factor  H, 
and  skin-friction  coefficient  Cf  which  are  presented  in  the  following  figures: 


Parameter 

8 

§1 

82 

H 

Cf 

44 

45 

46 

47 

48 

The  data  from  the  120"  test  section  are  presented  in  parts  (a),  (b),  and  (c)  (test 
section  diffuser  angles  of  0°,  14°,  and  7°,  respectively)  and  from  the  24"  test 
section  are  presented  in  parts  (d),  (e),  and  (f)  (test  section  diffuser  angles  of  0°, 
14°,  and  7°,  respectively)  of  each  of  the  above  figures. 

The  longer  120"  test  section  provided  a  fully-developed  turbulent 
boundary  layer  with  only  small  increases  in  the  x  =  100"  to  112"  measurement 
region  where  the  Reynolds  number  was  large  (Re  >  5x10^).  The  boundary-layer 
thicknesses  (8,  8j,  and  82)  for  the  120"  test  section  with  either  the  0°  or  the  7°  test 
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section  diffusers  did  not  change  significantly  with  variation  in  the  cavity 
configuration  (parts  (a)  and  (c)  of  Figures  44  through  46).  Depending  on  the 
selected  comparison,  these  data  were  within  ±  2%  to  ±  12%  for  a  given  plot.  Only 
the  four-cavity  configuration  with  the  14°  diffuser  shows  a  noticeable  increase  (up 
to  50%)  in  these  thicknesses  (part  (b)  of  Figures  44  through  46). 

The  24"  test  section  had  a  much  thinner  boundary  layer  which  is 
initially  barely  turbulent  since  the  Reynolds  number  ranged  from  350,000  to  about 
1,100,000.  The  boundary  layer  profiles  presented  in  Appendix  D  tend  to  agree 
with  the  1/7  power  law  which  is  typical  of  a  two-dimensional  turbulent  boundary 
layer.  The  boundary  layer  thicknesses  (8,  8j,  and  82)  for  the  24"  test  section  with 

either  the  0°  or  7°  test  section  diffuser  are  typically  within  ±  10%.  In  the  14°  test 
section  diffuser,  while  there  was  more  data  scatter,  there  tended  to  be  an  increase 

in  thickness  with  the  increasing  diffuser  angle.  The  shape  factor  H  (Figure  47) 
and  skin  friction  coefficient  Cf  (Figure  48)  do  not  show  any  noticeable  effect  of 

the  presence  of  the  cavities. 

The  boundary  layer  log-law  profiles  were  compared  (Figures  49  to 
51)  both  upstream  and  downstream  of  the  cavity  (or  cavities)  in  an  attempt  to 

identify  their  effect  on  the  flow.  It  is  shown  that  all  of  the  data  obtained  without  a 
diffuser  tend  to  agree  with  the  log  law  (equation  18)  where  k  =  0.41  and  C]  =  5.0. 

The  data  with  the  7°  and  14°  test- section  diffusers  asymptote  below  the  log-law 
defined  by  the  above  constants.  It  was  found,  as  shown  in  Figure  52,  that  if  Cj  is 

changed  to  3.5  for  the  7°  diffuser  and  to  3.0  for  the  14°  diffuser  that  these  data  for 
the  test-section  diffuser  geometries  agree  with  the  log-law  equation. 

The  effect  of  variation  of  the  cavity  configuration  (0,  1,  2,  or  4 
cavities)  on  the  boundary-layer  profiles  at  a  given  x  station  is  presented  in  Figures 
53,  54,  and  55.  The  test-section  entrance  profiles  (Figure  53)  demonstrate 
boundary  layer  similarity  for  the  basic  test  section  and  for  the  7°  test  section 
diffuser.  These  profiles  are  in  agreement  with  the  1/7  power  law  profile.  There  is 
a  significant  boundary  layer  profile  variation  from  the  1/7  power  law  for  the  14° 
test-section  diffuser.  The  data  from  Figure  54,  which  are  for  the  downstream  x 
station,  are  presented  in  Figure  55  using  the  universal  velocity  profile.  This  data 
format  demonstrates  a  modest  variation  in  the  boundary-layer  profiles  which  are 
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typical  of  an  adverse  pressure  gradient.  The  single  and  double  cavities  show  only 
a  small  deviation  from  the  profile  measured  with  no  cavity.  The  four-cavity 
configuration  profile  in  most  cases  represents  a  more  adverse  boundary-layer 
profile  than  any  other  configuration. 


Figure 

test  section  length 

x  station 

test  section 
diffuser  angle,  deg. 

53  a 

120 

-  101 

0 

53b 
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14 

53  c 
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53  e 

24 

~  5 
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54  a 
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16 
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54  f 
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b.  Effects  on  the  Surface  Pressure  Distribution 


The  effect  of  the  test-section  diffuser  on  the  surface-pressure 
distribution  was  discussed  on  page  82.  The  focus  of  the  present  section  is  the 
evaluation  of  the  effect  of  cavities  on  the  pressure  distribution  on  the  wind-tunnel 
ceiling.  Several  configurations  were  measured  several  times  to  determine  the  data 
repeatibility.  An  example  is  presented  in  Figure  56  where  the  ceiling  pressure 
distributions  from  five  runs  of  configuration  8431c  (120"  test  section  with  a  2" 
deep  cavity  covered  by  a  door  located  between  x  =  102"  and  104")  are  plotted 
along  with  their  average.  Immediately  downstream  of  the  cavity  the  pressure 
coefficient  is  seen  to  be  depressed.  This  was  probably  caused  by  a  misalignment 
of  the  door  with  the  ceiling.  Most  of  the  pressure  coefficients  are  within  ±  0.02  of 
the  average  pressure  coefficient  curve  shown  on  Figure  56. 

The  effect  of  the  cavities  in  the  24"  test  section  is  presented  in 
Figures  57  and  58.  The  basic  test  section  shows  a  pressure  coefficient  magnitude 
of  approximately  -0.02;  ahead  of  the  cavities  there  is  a  similar  pressure  coefficient. 
Downstream  of  the  cavity  the  pressure  coefficients  drop  to  about  -0.06.  With  both 
test  section  diffusers  (Figures  59  through  61)  in  the  24"  test  section  the  general 
data  trend  is  unaffected  by  the  presence  of  the  cavities. 

The  effect  of  the  cavities  in  the  120"  test  section  is  presented  in 
Figures  62  through  66.  Generally  there  is  seen  to  be  little  effect  of  the  cavities  on 
the  pressure  distributions.  The  one  exception  is  the  effect  of  the  2"  cavity 
downstream  of  that  cavity  where  the  pressure  coefficient  is  reduced  to  about  -0.04. 
With  both  test  section  diffusers  there  is  more  scatter  in  the  data  without  a 
consistent  trend. 

In  conclusion,  all  of  the  surface  pressure  data  suggest  that  the 
presence  of  the  cavities  has  no  consistent  effect.  It  may  be  concluded  that  cavities 
may  be  deployed  without  significantly  modifying  the  resultant  pressure 
distribution. 
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(a)  120"  test  section. 

Figure  44.  Effect  of  the  cavity  configuration  on  the  boundary-layer  thickness 


(b)  120"  test  section  with  the  14°  test-section  diffuser 
Figure  44.  Continued. 


(c)  120"  test  section  with  the  7°  test-section  diffuser. 
Figure  44.  Continued. 


(d)  24"  test  section 
Figure  44.  Continued. 
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(a)  120"  test  section. 

Figure  45.  Effect  of  the  cavity  configuration  on  the  displacement  thickness. 


(b)  120"  test  section  with  the  14°  test-section  diffuser. 
Figure  13.  Continued. 


(c)  120"  test  section  with  the  7°  test-section  diffuser. 
Figure  45.  Continued. 


(d)  24"  test  section 
Figure  45.  Continued. 
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(e)  24"  test  section  with  the  14°  test-section  diffuser. 
Figure  45.  Continued. 


(f)  24"  test  section  with  the  7°  test  section  diffuser. 
Figure  45.  Concluded. 


x,  in. 

(a)  120"  test  section. 

Figure  46.  Effect  of  the  cavity  configuration  on  the  momentum  thickness. 


x,  in. 


(b)  120"  test  section  with  the  14°  test-section  diffuser. 
Figure  46.  Continued. 
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(e)  24"  test  section  with  the  14°  test-section  diffuser. 
Figure  46.  Continued. 


(f)  24"  test  section  with  the  7°  test-section  diffuser. 
Figure  46.  Concluded. 
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(e)  24"  test  section  with  the  14°  test-section  diffuser. 
Figure  47.  Continued. 


(f)  24"  test  section  with  the  7°  test-section  diffuser. 
Figure  47.  Concluded. 


0 


2< 


0.006 


0.004  h 


0.002  h 


0 


■ 

no.  c 

_ n  _ 

;avities,  D,  in. 

i 

--*--1,2 

.  V 

-■A”- 

*■ . 

4,  1 

L  , 

—  iJ 

—  i,  i 

.  0-c>; 
: 

1 

ggjj 

■ 

m 

i  ■  ■  i  i 

100 


104  108 

x,  in. 


112 


(a)  120"  test  section. 

Figure  48.  Effect  of  the  cavity  configuration  on  the  skin-friction  coefficient. 


(b)  120"  test  section  with  the  14°  test-section  diffuser. 
Figure  48.  Continued. 
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(e)  24"  test  section  at  x  ~  5"  with  the  14°  test-section  diffuser. 
Figure  49.  Concluded. 
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(a)  120”  test  section  at  x  ~  108". 

Figure  50.  Effect  of  the  cavity  geometry  on  the  boundary-layer  log-law  profile. 
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(b)  120"  test  section  at  x  «  108"  with  the  14°  test-section  diffuser. 

Figure  50.  Continued. 


y+ 


(c)  120"  test  section  at  x  ~  108"  with  the  7°  test-section  diffuser. 

Figure  50.  Continued. 
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(e)  24"  test  section  at  x  ~  12"  with  the  14°  test-section  diffuser. 
Figure  50.  Continued. 


(f)  24"  test  section  at  x  »  12"  with  the  7°  test-section  diffuser. 
Figure  50.  Concluded. 


Figure  51.  Effect  of  the  cavity  geometry  on  the  boundary-layer  log-law  profile. 
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(b)  120"  test  section  at  x  =  1 12"  with  the  14°  test-section  diffuser. 

Figure  51.  Continued. 


y+ 

(c)  120"  test  section  at  x  =  112"  with  the  7°  test-section  diffuser. 

Figure  51.  Continued. 
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(f)  24"  test  section  at  x  =  16"  with  the  7°  test-section  diffuser. 
Figure  51.  Concluded. 
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(a)  7°  test-section  diffuser  data  compared  with  log-law  equation  using  Cj  =3.5. 

Figure  52.  Comparison  of  120"  test-section  experimental  data  at  x  =  112"  with  the 
log-law  equation  with  Cj  changed  from  5.0. 
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(b)  14°  test-section  diffuser  data  compared  with  log-law  equation  using  Cj  =  3 

Figure  52.  Concluded 


(a)  120"  test  section  at  x  ~  101". 

Figure  53.  Effect  of  the  cavity  configuration  on  the  boundary-layer  profile. 
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(b)  120"  test  section  at  x  ~  101"  with  the  14°  test-section  diffuser. 

Figure  53.  Continued. 
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(c)  120"  test  section  at  x  ~  101"  with  the  7°  test-section  diffuser. 

Figure  53.  Continued. 
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(a)  120"  test  section  at  x  =  1 12  . 

Figure  54.  Effect  of  the  cavity  configuration  on  the  boundary-layer  profile. 
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(b)  120"  test  section  at  x  =  1 12"  with  the  14°  test-section  diffuser. 

Figure  54.  Continued. 
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(e)  24"  test  section  at  x  =  16"  with  the  14°  test-section  diffuser. 
Figure  54.  Continued. 


(f)  24"  test  section  at  x  =  16"  with  the  7°  test-section  diffuser. 
Figure  54.  Concluded. 
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(c)  120"  test  section  at  x  =  1 12  with  the  7°  test-section  diffuser. 
Figure  55.  Continued. 
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(d)  24"  test  section  at  x  =  16". 
Figure  55.  Continued. 
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(e)  24"  test  section  at  x  =  16"  with  the  14°  test-section  diffuser. 
Figure  55.  Continued. 
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(f)  24"  test  section  at  x  =  16"  with  the  7°  test-section  diffuser. 
Figure  55.  Concluded. 
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Figure  58.  Effect  of  the  multiple  cavities  on  the  ceiling  pressure  distribution  (24 

test  section  with  D  =  1"). 


0.5 


Figure  60.  Effect  of  the  multiple  cavities  on  the  ceiling  pressure  distribution  (24" 
test  section,  14°  test-section  diffuser  with  D  =  1 "). 


Figure  61 .  Effect  of  the  multiple  cavities  on  the  ceiling  pressure  distribution  (24" 
test  section,  7°  test-section  diffuser  with  D  =  1 "). 
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Figure  62.  Effect  of  the  single  cavity  size  on  the  ceiling  pressure  distribution 
(120"  test  section  with  D  =  1"  or  2"). 


test  section  with  D  =  1"). 


Figure  64.  Effect  of  the  single  cavity  size  on  the  ceiling  pressure  distribution 
(120"  test  section,  14°  test-section  diffuser  with  D  =  1"  or  2"). 


Figure  65.  Effect  of  the  multiple  cavities  on  the  ceiling  pressure  distribution  (120" 
test  section,  14°  test-section  diffuser  with  D  =  1"). 
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Figure  66.  Effect  of  the  multiple  cavities  on  the  ceiling  pressure  distribution  (120" 
test  section,  7°  test-section  diffuser  with  D  =  1"). 


3.  Cavity  With  a  Door 


As  described  in  the  introduction,  the  effect  of  cavities  are  being  studied  as  a 
means  for  increasing  airfoil  camber  to  provide  a  new  high  lift  device.  One  aspect 
of  this  application  is  the  effect  of  surface  non-uniformities  (which  could  be  caused 
by  cavity  doors)  on  the  boundary  layer  and  on  the  surface  pressure  distribution. 
Tests  were  conducted  which  compared  a  plain,  smooth  ceiling  with  a  ceiling 
which  had  a  two-inch  square  cavity  which  was  closed  by  a  flush  door.  The  door 
was  mounted  by  a  pivot  at  its  upstream  end.  When  the  wind  tunnel  was  not 
operating,  the  closed  door  was  adjusted  to  be  flush  with  the  adjacent  tunnel 
ceiling.  When  the  tunnel  was  turned  on  and  operated  at  a  freestream  velocity  of 
120  feet  per  second,  it  was  noted  that  the  door  deflected  into  the  flow 


147 


approximately  1/16  to  1/8  inch  at  its  downstream  edge.  The  larger  deflections 
occurred  when  one  of  the  test-section  diffusers  was  present  ( (j)  =  7°  or  14°). 

For  these  tunnel  geometries,  the  measured  boundary-layer  thicknesses  are 
presented  in  Figure  67  and  the  pressure  distributions  are  presented  in  Figure  68.  In 
each  figure,  part  (a)  presents  data  from  the  24"  test  section  and  part  (b)  has  data 
from  the  120"  test  section.  The  data  for  the  plain  ceiling  are  presented  by  the  solid 
symbols  and  the  data  for  the  cavity  with  a  door  are  presented  by  open  symbols. 
The  boundary-layer  thicknesses  for  the  24"  test  section  cavity  with  a  diffuser  and  a 
cavity  with  a  closed  door  (Figure  67(a))  are  thicker  than  those  for  the  plain  ceiling. 
The  boundary  layer  thickness  for  the  closed  cavity  with  the  7°  test-section  diffuser 
increased  much  more  rapidly  than  any  of  the  other  configurations.  This  difference 
is  probably  due  to  the  increased  projection  of  the  cavity  door  trailing  edge  into  the 
freestream,  which  was  a  significant  fraction  of  the  8  of  0.25  to  0.45  inch  in  this 
test  section.  The  data  for  the  120"  test  section  (Figure  67(b))  show  a  more 
noticeable  effect  of  the  cavity  door  projection  even  though  the  door  projection  was 
a  small  fraction  of  the  measured  8,  which  was  greater  than  an  inch.  This  indicates 
that  poor  fit  of  a  cavity  door  could  increase  skin  friction  drag.  It  is  seen  in  Figure 
68  that  there  were  differences  in  the  longitudinal  pressure  distribution  for  both  test 
sections.  This  indicates  that  the  surface  non-uniformity  had  an  effect  on  the 
inviscid  freestream  which  determines  the  surface  pressure  distribution. 

A  second  set  of  tests  was  conducted  which  compared  the  flow  in  the 
presence  of  a  two-inch  square  open  cavity  with  a  similar  two  inch  cavity  (  the 
closed  cavity  used  in  Figures  67  and  68)  where  a  door  was  opened  into  the  cavity 
to  form  a  front  cavity  wall  of  an  open-door  cavity.  The  open  cavity  without  a  door 
had  a  smooth  surface  upstream  of  the  cavity.  In  contrast,  the  open  cavity  door  had 
a  small  gap  of  about  1/16  inch  from  the  upstream  end  of  the  ceiling  surface.  These 
results  are  presented  in  Figure  69  (boundary  layer  thicknesses)  and  Figure  70 
(longitudinal  pressure  distribution).  While  these  results  appear  to  be  similar  to 
those  presented  for  the  closed  cavity  in  Figures  67  and  68,  there  are  notable 
differences.  The  boundary-layer  thicknesses  for  all  of  the  open-door  cavity 
configurations  showed  an  increase  with  longitudinal  distance  similar  to  the 
increase  of  the  closed  cavity  which  had  the  largest  increase,  the  7°  test-section 
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x,  in 


(a)  24"  test  section 

Figure  67.  Effects  of  closed  cavity  door  and  test-section  diffusers  on  the  variation 
of  boundary-layer  thickness  as  a  function  of  longitudinal  location. 


x,  in 

(b)  120"  test  section 
Figure  67.  Concluded 
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(a)  24"  test  section 

Figure  69.  Effect  of  open  cavity  door  and  test-section  diffuser  angle  on  the  ceiling 
boundary-layer  thickness  as  a  function  of  longitudinal  location. 


x,  in 

(b)  120"  test  section 
Figure  69.  Concluded 
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(a)  24"  test  section 

Figure  70.  Effect  of  open  cavity  configuration  and  test-section  diffuser  angle  on 
the  ceiling  longitudinal  pressure  distribution. 


(b)  120"  test  section 
Figure  70.  Concluded 
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configuration.  The  reasons  for  this  result  have  not  been  identified.  In  contrast,  it 
is  seen  that  with  the  open  cavities  (Figure  70)  there  were  smaller  differences  in  the 
longitudinal  pressure  distribution  for  either  test  section  length  than  that  shown  by 
the  closed  cavity  and  plain  ceiling.  These  results  suggest  that  the  open  cavities 
produce  less  change  in  the  surface  pressure  distributions  than  do  surface  non¬ 
uniformities  in  an  otherwise  plain  ceiling. 
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IV.  INCOMPRESSIBLE  NAVIER-STOKES  COMPUTATIONS 


There  are  some  simple  flows  related  to  the  present  investigation  which  may 
be  represented  by  empirical  correlations.  One  example  is  a  two-dimensional, 
turbulent  boundary-layer  profiles  using  the  1/7  power  law  (u/U,*,  =(y/S/1/7^). 
More  complete  flow  solutions  require  more  complete  sets  of  equations.  There  is  a 
range  of  equation  sets  of  increasing  complexity  available  for  the  solution  of  fluid 
flows.  All  of  these  equation  sets  may  be  obtained  as  simplifications  of  the  Navier- 
Stokes  equations.  The  particular  equations  selected  for  a  problem  of  interest 
depends  on  the  characteristics  of  that  flow.  The  simplest  of  these  equations  is  the 
linear  potential  Laplace  equation  which  assumes  inviscid,  incompressible, 
irrotational  flow. 


V2O  =  0  (26) 

When  there  are  viscous  effects  more  complex  equation  sets  are  needed.  Two- 
dimensional,  viscous  flows  in  a  conventional,  subsonic  wind  tunnel  test  section 
may  usually  be  accurately  computed  using  the  boundary-layer  equations 
(Equations  11  through  13).  These  equations  describe  the  attached  boundary  layer 
near  the  walls  which  develops  adjacent  to  the  inviscid  outer  flow,  which  is 
representative  of  the  bulk  of  the  freestream  flow. 

The  boundary-layer  concept  of  Prandtl  assumes  that  the  viscous  layer  is 
small  (thin  boundary  layer)  relative  to  the  streamwise  flow  distance  (8/L «1). 
As  a  consequence,  the  Navier-Stokes  equations  may  be  simplified  to  obtain  the 
boundary  layer  approximation  for  steady,  two-dimensional,  incompressible, 
constant  fluid-property  flow  [Ref.  90].  The  present  experimental  investigation 
adds  an  adverse  pressure  gradient  and  one  or  more  cavities  to  produce  a  flow 
which  is  more  complicated  than  those  which  can  be  resolved  by  the  boundary- 
layer  equations.  Especially  difficult  is  the  resolution  of  the  recirculating  flow  in  a 
cavity  and  its  interaction  with  the  freestream  in  an  adverse  pressure  gradient.  For 
example,  in  the  120-inch  test  section  in  the  present  investigation,  the  cavity  flow 
had  a  viscous  layer  whose  thickness  was  nearly  the  same  as  the  cavity  depth.  In 
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contrast,  for  the  24  inch  test  section  the  viscous  boundary  layer  thickness  was 
much  less  than  the  cavity  depth. 

Based  on  the  present  literature  survey,  neither  experimental  nor 
computational  data  have  been  obtained  for  these  flows  prior  to  the  present 
investigation.  The  present  experiment  measured  the  flow  in  the  vicinity  of  several 
cavity  configurations  with  a  freestream  Mach  number  of  the  order  of  0.10.  For 
flows  such  as  these,  the  incompressible  Navier-Stokes  equations  are  a  more 
suitable  representation  than  the  boundary  layer  equations.  The  incompressible 
Navier-Stokes  equations  are  derived  from  the  Reynolds-averaged  Navier-Stokes 
equations  by  assuming  an  incompressible  flow  (i.e.,  M  =  0  and  a  =  °°). 

Since  the  temperature  was  essentially  constant  in  the  present  experimental 
investigation  the  energy  equation  is  not  required.  The  two-dimensional, 
incompressible  Navier-Stokes  equations  in  Cartesian  coordinates  become  the 
following: 


.  .  .  du  dv 

continuity  equation,  —  +  —  = 

dx  dy 


(27) 


x  momentum  equation, 


du  du 
— +  u  — 
dt  dx 


1  dp  fd2u,d2u) 
“ — M)  ■■■-■■+  9 

p  dx  ^  dx  dy  J 


(28) 


y-momentum  equation, 


dv  dv 
— +  u— +  v 
dt  dx 


dv 

dy 


1  dp  ,  J  d2v  ,  d2v^ 

—  _ — |-  — y  H - y 

p  dy  ^dxz  dy 


(29) 


A  numerical  solution  of  these  equations  may  be  obtained  for  example  using  either 
a  finite-element  or  finite-difference  scheme,  a  suitable  grid,  and  the  definition  of 
appropriate  boundary  conditions. 
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A.  INCOMPRESSIBLE  NAVIER-STOKES  CODE  (ins2d) 


The  ins2d  computer  code  [Ref.  91  and  92]  numerically  solves  the  two- 
dimensional,  incompressible  Navier-Stokes  equations  in  a  time-accurate  manner. 
In  a  pseudocompressibility  formulation,  a  time  derivative  of  pressure  is  added  to 
the  continuity  equation, 


9p  a( 3u  3v^| 
— +  (5  —  +  — 
dx  ydx  dy  j 


(30) 


Note  that  |3  is  the  pseudocompressibility  constant  and  x  represents  a  pseudotime 
which  is  not  related  in  any  way  to  physical  time.  Equation  30  is  then  combined 
with  the  momentum  equations  to  obtain  the  following  equations  in  two- 
dimensional  Cartesian  coordinates. 


J;D+i^E-E»)+!;(F-F’)=0  (31) 
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The  convective  terms  are  upwind-differenced  using  a  Roe  [Ref.  93]  flux- 
difference  split  approach  that  has  uniformly  high  accuracy  throughout  the  interior 
grid  points.  This  scheme  was  derived  as  an  approximate  Riemann  solver  for  the 
compressible  flow  equations.  In  this  formulation  the  Reynolds  stress  has  been 
approximated  as  a  function  of  the  strain-rate  tensor,  and  thus  represents  a  sum  of 
the  kinematic  viscosity  and  the  turbulent  eddy  viscosity.  Constant  kinematic 
viscosity  is  assumed  for  simplicity.  The  viscous  fluxes  are  differenced  using 
second-order  accurate  central  differences. 

The  system  of  equations  (equations  31)  is  solved  numerically  using  an 
implicit  line-relaxation  scheme.  Application  of  a  first-order  backward  Euler 
formula  to  equations  31  yields  the  delta-form  equation 


1 

JAx 


1  + 


f— T 


(32) 


Here  the  superscript  n  is  the  pseudotime  iteration  count  and  R  is  the  residual 
vector.  Unsteady  problems  are  solved  with  the  use  of  subiterations  in  pseudotime 
at  each  physical  time  step. 

At  the  inflow  boundary  there  is  one  characteristic  wave  traveling  out  of  the 
computational  domain  the  since  fluid  is  traveling  into  the  domain.  When  the 
inflow  velocity  profile  is  not  known,  there  is  a  vector  of  variables  which  is  held 
constant  and  defined  as 
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(33) 
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For  the  present  calculation  the  vertical  inflow  velocity  component  v  was  set  to 
zero  and  the  horizontal  inflow  velocity  u  was  set  to  unity  at  all  of  the  inflow 
plane  grid  points  except  at  the  floor  and  ceiling  grid  points  where  both  were  zero. 

At  the  outflow  boundary  there  are  two  characteristic  waves  traveling  out  of 
the  computational  domain  since  fluid  is  also  leaving  the  domain.  This 
computation  used  a  specified  static  pressure  at  the  outflow  plane  which  was 
defined  as 
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The  tangential  flow  boundary  conditions  on  the  tunnel  floor  and  ceiling  were 
obtained  by  setting  both  the  u  and  v  velocity  components  to  zero. 


B.  GRID  GENERATION 


The  grid  was  generated  using  several  computers  and  codes.  Several 
FORTRAN  codes  were  written  to  generate  both  the  freestream  flow  grid  between 
the  wind  tunnel  ceiling  and  the  floor  (sjtgg.f  and  rtgg.f)  and  the  cavity  grid  (cgg.f 
and  cgga.f).  Then  these  grids  were  coupled  by  the  PEGSUS  code  [Ref.  94]  which 
uses  an  overset  scheme.  Typically  these  codes  were  run  interactively.  For  the 
basic  wind  tunnel  flow  computations  the  grid  from  code  sjtgg.f  was  used  directly 
in  the  ins2d  code.  These  FORTRAN  programs  were  solved  using  a  Silicon 
Graphics  (SGI)  Challenger  L  workstation.  This  computer  has  4  processors,  a 
MIPS  R4400  microprocessors  with  a  clock  speed  of  150  MHz,  256  MB  random- 
access  memory  and  2GB  random-access  disk  storage. 
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For  a  boundary  layer  flow,  the  largest  velocity  gradient  occurs  adjacent  to 
the  walls  and  normal  to  the  flow  direction.  For  the  present  ins2d  Navier-Stokes 
flow  solver,  it  is  recommended  that  the  finest  spacing  near  the  wall  lead  to  a 
minimum  non-dimensional  inner  layer  distance  normal  to  the  wall  y+  of 
approximately  1.  The  y+  is  reduced  by  reducing  the  grid  spacing  near  the  tunnel 
wall.  In  Anderson,  Tannehill,  and  Pletcher  [Ref.  90]  several  simple 
transformations  are  presented  which  can  be  used  to  cluster  the  grid  in  regions  of 
large  gradients  such  as  boundary  layers.  These  transformations  are  part  of  a 
general  family  proposed  by  Roberts  [see  Ref.  90,  pp.  247-250].  One 
transformation  was  developed  to  refine  the  mesh  near  the  walls  of  a  two- 
dimensional  tunnel  using  equations 

x  =  x 
and 


(P  +  2a)[(P  + 1)  /  ((3  -  l)](y  a)/(1  a)-p  +  2a 

y_h  f  _  (y-a)/(l-a)  ] 

(2a  +  lK  l  +  [(p  +  l)/(p-l)]  [ 


(35) 


where  h  is  tunnel  height,  a  =  1/2  refines  the  grid  equally  near  y  =  0  and  y  =  h,  and 

(  8Y1/2  8 

P  =  1  —  for  0<  —  <1  (36) 

H  l  h)  h 

Coordinates  x  and  y  represent  a  uniform  Cartesian  grid  and  coordinates  x  and  y 
represent  a  dusted  grid.  A  computer  program  (rtgg.f)  was  written  to  generate  the 
grid  for  the  wind  tunnel  used  in  the  experimental  investigation. 

An  example  of  Roberts'  transformation  grid  clustering  is  shown  in  Figure 
71(a)  for  the  region  between  the  tunnel  floor  and  the  tunnel  centerline.  The  region 
from  the  tunnel  centerline  to  the  tunnel  ceiling  is  the  reflection  of  the  clustering 
shown  in  Figure  71.  The  horizontal  scale  is  a  uniformly-spaced  coordinate  T[i 
non-dimensionized  by  the  scale  factor  SF,  and  the  vertical  scale  is  clustered  grid 
scale  yi  norm  which  is  then  converted  from  the  0  to  1  scale  to  the  dimensions  of 
the  tunnel  half  height.  The  grid  clustering  is  determined  by  the  value  of  P  which  is 
shown  to  range  from  2,  for  the  least  clustering,  to  1.001,  for  the  greatest  clustering. 
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The  largest  two  values,  2  and  1.2,  provide  only  a  limited  clustering.  In  the  next 
section  results  from  the  ins2d  code  will  be  presented  for  the  three  finest  scale 
factors  shown  in  Figure  71(a):  1.05,  1.01,  and  1.001.  These  results  will  be  used  to 
evaluate  the  usefulness  of  the  Roberts  grid-transformation  scheme  for  representing 
these  boundary-layer  flows. 


An  alternative  to  Roberts'  transformation-clustering  scheme  was  developed 
by  the  author  using  a  tangential  grid-spacing  scheme.  This  transcendental  function 
was  used  to  transform  a  uniform  spacing  to  a  variable  grid  spacing.  The  fineness 
of  the  grid  spacing  was  determined  by  a  tangential  grid-spacing  scale  factor  SF. 
Initially  the  scheme  divided  a  distance  into  a  uniform  spacing  where 


gy,  _  ^lmax  til 

NC-1 


(37) 


using 

r)!  =  tan-1(-SF)  and  Timax  =tan_1(SF)  (38) 


The  value  of  the  scale  factor  is  increased  when  a  finer  grid  is  desired  at  the  outer 
grid  edges.  For  the  present  investigation,  SF  was  typically  2  in  the  streamwise 
direction  and  ranged  from  16  to  72  normal  to  the  freestream  direction. 


As  an  illustration  of  the  final  steps  only  the  y  coordinate  equations  are  used. 
Similar  calculations  are  done  in  the  x  direction.  For  the  next  step,  a  normalized 
spacing  from  -1  to  1  was  computed  as 


tan  1  rjj 
tan"1  Tlmax 


(39) 


Finally  the  physical  coordinates  were  computed  as 

yi=yorg+^L  (40) 


The  fraction  in  the  second  term  scales  length  from  y4  to  a  range  from  0  to  1 . 
Since  yorg  is  the  minimum  physical  dimension  and  L  is  the  dimension  length, 
the  yi  values  cover  the  entire  length  of  this  direction.  Program  sjtgg.f  was  used 
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to  algebraically  calculate  a  two-dimensional  grid  for  the  wind  tunnel  freestream 
including  the  tunnel  ceiling  and  floor,  with  or  without  the  test  section  diffuser. 
Program  cgg.f  was  used  to  generate  each  cavity  grid.  Both  calculations  clustered 
the  grid  using  a  tangential  spacing  scheme. 

The  normalized  spacings  obtained  from  this  procedure  is  illustrated  for 
several  scale  factors  in  Figure  71  (b).  The  increased  clustering  near  the  walls  is 
readily  apparent.  At  the  highest  scale  factors  the  grid  tends  to  be  sparse  near  the 
center  of  the  physical  space.  For  the  streamwise  direction  the  flow  tends  to  a  zero 
gradient  in  the  potential  flow  region  away  from  the  boundary  layer.  As  a 
consequence  the  coarser  spacing  is  adequate  near  the  tunnel  centerline.  The  grid 
size  variation  is  well  behaved  for  the  present  calculations  because  it  changes  by 
less  than  1 .2  between  adjacent  grid  cells.  This  grid  clustering  scheme  is  evaluated 
in  the  next  section  to  determine  its  suitability  for  the  present  investigation. 


h/SF 

(a)  Roberts  grid  transformation. 

Figure  7 1 .  Clustered  grid  spacing  as  a  function  of  a  uniform  spacing  for  several 

scale  factors. 
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(b)  Tangential  grid  transformation. 
Figure  71.  Concluded 


The  tangential  grid  spacing  is  used  for  both  directions  in  each  cavity  and  for 
both  tunnel  directions.  When  the  tunnel  has  one  or  more  cavities,  the  tangential 
grid  spacing  in  the  stream  wise  direction  must  be  modified  to  satisfy  an  adjacent 
grid  spacing  variation  of  1.2  or  less.  This  requirement  was  satisfied  as  part  of  the 
program  cgga.f.  First  this  program  generates  the  grid  for  the  specified  cavity. 
Then  the  tunnel  grid  generated  in  sjtgg.f  is  read  and  modified.  In  the  streamwise 
direction  the  grid  spacing  is  reduced  immediately  upstream  and  downstream  of  the 
cavity.  Generally  the  two  grids,  freestream  and  cavity,  are  compatibily  spaced  in 
the  vertical  direction  and  require  no  modification.  If  modification  is  needed  to 
satisfy  the  1.2  adjacent  spacing  criterion  the  program  cgga.f  accomplishes  this 
change  in  the  freestream  grid.  If  more  than  one  cavity  is  present  then  cgga.f 
continues  on  to  add  these  cavities  to  complete  the  grid. 
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C.  CFD  RESULTS 


After  the  grid  generation  was  completed,  the  CFD  code  computed  the  fluid 
flow  solution  and  analysis  of  the  results.  After  the  grids  were  prepared,  two 
additional  input  files  were  prepared  which  provided:  (1)  the  boundary  conditions 
and  (2)  the  ins2d  input  data  for  the  desired  code  options.  The  present  input 
boundary  conditions  were  (a)  an  inflow  boundary  with  a  constant  velocity  and 
characteristic  relation  for  pressure,  (b)  an  outflow  boundary  using  characteristic 
relations  for  velocity  and  constant  static  pressure,  (c)  the  ceiling  with  a  no-slip 
wall  with  the  wall  normal  vector  pointing  in  the  negative  computational  direction, 
and  (d)  the  floor  with  a  no-slip  wall  with  the  wall  normal  vector  pointing  in  the 
positive  computational  direction.  About  40  items  of  data  make  up  the  ins2d.in 
file.  These  data  are  defined  in  a  user  file  described  by  the  ins2d  code  and  will  not 
be  discussed  in  detail  in  this  paper. 

The  ins2d  code  was  run  on  a  Cray  C-90  supercomputer  with  8  processors, 
256  Mw  of  random  access  memory  (ram),  8  Gw  of  disk  storage.  The  floating¬ 
point  processor  speed  was  1  GFLOP  which  is  about  20  faster  than  the  Challenge  L 
workstation  used  for  grid  generation.  A  typical  case  with  a  201  by  321  grid  used 
about  7.75  Mw  of  ram  and  required  about  2.2  seconds/iteration  or  345  x  10"7 
seconds/iteration/grid  point.  The  convergence  criterion  is  based  on  the  magnitude 
of  the  change  in  the  residual  term  from  the  previous  iteration.  It  was  found  that  if 
a  case  converged  it  took  40  to  150  iterations.  If  a  case  took  more  iterations  it  did 
not  converge.  Some  cases  were  run  for  as  many  as  20,000  iterations  without 
convergence. 

For  data  analysis,  the  desired  ins2d  results  were  then  selected  from  the 
solution  file  using  a  postprocessing  program  called  plot3d  on  the  Cray  C-90. 
These  results  were  then  transferred  as  a  group  of  ASCII  data  files  using  a  file 
transfer  protocol  to  an  Apple  Macintosh  Ilci  personal  computer.  This  computer 
used  a  Motorola  68030  microprocessor  which  had  a  clock  speed  of  40  MHz  and  an 
8  MB  random-access  memory  and  an  80  MB  hard  disk.  In  most  cases  these  data 
were  read  by  a  plotting  program  called  Kaleidagraph. 
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1.  Wind  Tunnel  Flow  Computation 


Several  grid  geometries  were  examined  for  the  basic  24"  test  section 
(identified  as  COOOO)  and  for  the  120"  test  section  (identified  as  C8000)  used  in 
the  experimental  investigation.  For  this  phase  of  the  computational  investigation 
both  the  Robert's  transformation  and  the  tangential  transformation  were  used  to 
generate  grids  to  represent  the  flow  in  the  wind  tunnel  for  both  the  24-inch  test 
section  and  the  120-inch  test  section.  Also  examined  for  both  test  section  lengths 
were  the  number  of  grid  points  in  the  y  direction,  normal  to  the  freestream.  For 
the  24-inch  test  section,  the  effect  of  the  longitudinal  length  of  the  x  grid  was 
varied. 


a.  24  Inch  Test  Section 


The  grid  was  used  to  represent  the  24  inch  wind  tunnel  test  section 
(COOOO)  from  the  floor  to  the  ceiling  and  between  the  longitudinal  stations  of  -16 
to  24  inches.  A  few  preliminary  cases  were  run  for  a  shorter  longitudinal  grid 
which  ran  from  0  to  24  inches.  It  was  found  that  this  provided  boundary  layer 
thicknesses  which  were  too  thin  in  relation  to  the  experimental  data.  The  grid 
clustering  for  Roberts'  transformation  is  determined  by  the  value  used  for  the 
stretching  parameter  (3  (equation  36).  For  the  24"  test  section  the  boundary  layer 
thickness  5  ranged  from  0.22  to  0.45  inch  depending  on  the  longitudinal  survey 
location.  This  lead  to  a  recommended  stretching  parameter  of  about  1.01.  In  the 
ins2d  code  the  resultant  grid  resulted  in  a  minimum  y+  value  of  10.0  which  is 
larger  than  the  recommended  value  of  about  1. 

Using  201  grid  points  in  the  y  direction,  normal  to  the  freestream, 
three  values  of  (3  were  used  which  resulted  in  a  reduction  in  the  minimum  y+  from 
14.0  to  1.8.  An  additional  grid  with  401  y  grid  points  was  also  used  to  obtain  a  y+ 
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value  near  the  desired  value  of  1.  The  following  table  outlines  the  Robert's 
transformation  cases  which  were  used  to  evaluate  these  grids. 


C0000  y  grid  points 

Robert's  grid  (3 

min.  y+ 

201 

1.05 

14.0 

201 

1.01 

10.0 

201 

1.001 

1.8 

401 

1.001 

0.9 

The  ins2d  code  iterates  until  either  the  maximum  number  of  iterations  are 
completed  or  the  computation  converges.  None  of  these  cases  converged  within 
5000  iterations.  This  undesirable  result  indicates  Lax's  equilibrium  theorem  [Ref. 
90]  is  not  being  satisfied.  This  theorem  states  that  given  a  properly-posed  initial- 
value  problem  and  a  finite-difference  approximation  that  satisfies  the  consistency 
condition,  stability  is  the  necessary  and  sufficient  condition  for  convergence. 

The  computed  boundary-layer  thicknesses  (8)  are  presented  in 

Figure  72  along  with  the  experimental  data  and  the  1/7  power  law  result.  The 
latter  two  values  are  shown  to  be  in  good  agreement.  The  computational  result 
from  the  coarsest  grids  produced  thicknesses  greater  than  experimental  data.  The 
computational  result  from  the  two  finest  grids  where  y+  is  near  1  show  thicknesses 
which  are  less  than  the  experimental  data.  These  inadequate  results  reflect  the  lack 
of  convergence  in  these  computations  and  suggest  that  Roberts'  transformation 
grids  do  not  provide  a  satisfactory  grid  for  the  present  problem.  Even  though  the 
thicknesses  do  not  show  good  agreement  with  the  experimental  data,  the 
boundary-layer  profiles  at  longitudinal  locations  5.5  (Figure  73(a))  and  16  inches 
(Figure  73(b))  show  agreement  between  the  experimental  and  the  computational 
results.  This  suggests  that  the  profile  is  a  less  useful  parameter  for  evaluating  the 
adequacy  of  the  computational  results. 
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data  P  sf  y+  y  grid 


Figure  72.  Boundary-layer  thickness  as  a  function  of  longitudinal  tunnel  location 
computed  by  the  ins2d  for  24"  tunnel  using  Roberts’  transformation  grid. 
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(a)  x  =  5.5". 

Figure  73.  Boundary-layer  profile  computed  by  the  ins2d  for  24"  tunnel  using 

Roberts'  transformation  grid. 
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(b)  x  =  16". 

Figure  73.  Concluded. 

The  tangential-grid  scheme  was  also  used  for  the  24  inch  test  section 
tunnel  (C0000).  Scale  factors  of  32  and  64  were  selected  to  achieve  y+  values 
which  approached  1.  The  number  of  y  grid  points  was  increased  from  201  to  401 
to  obtain  a  minimum  value  of  1 .2.  In  contrast  to  the  Robert's  transformation  all  of 
these  cases  converged  in  from  48  to  76  iterations.  These  rapid  solutions  enabled 
the  cases  to  be  run  in  the  debug  queue  on  the  Cray  C-90  computer  in  less  than 
three  minutes.  This  rapid  convergence  provides  confidence  that  these  grids 
provide  consistent,  stable  computational  results. 


Tangential  grid 
scale  factor 

C0000  y  grid  points 
y+  /  (iterations  for  convergence) 

201 

401 

32 

4.8/53 

2.4/76 

64 

2.4  /  48 

1.2/75 

- 1/7  power  law 

- ins2d  p= 1.001 
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Examples  of  the  computed  boundary  layer  thicknesses  using  the  tangential-grid 
scheme  are  presented  in  Figures  74(a)  for  201  y  grid-points  and  in  Figure  74(b) 
for  401  y  grid-points.  Both  figures  present  data  for  scale  factors  of  32  and  64.  In 
both  figures  there  is  excellent  agreement  between  the  computational  results,  the 
experimental  data,  and  the  1/7  power  law  at  longitudinal  locations  at  5.5  inches 
and  12  inches.  The  reasons  for  the  disagreement  at  the  16  inch  station  is 
unknown.  It  is  possible  there  may  be  an  influence  of  the  diffuser  downstream  of 
the  test  section  on  the  downstream  end  of  the  test  section.  This  diffuser  was  not 
modelled  in  these  calculations.  Figure  75  presents  comparisons  of  computational 
results,  the  1/7  power  law,  and  the  experimental  data  for  the  boundary  layer 
profiles.  All  of  these  results  show  agreement.  These  comparisons  for  the  24  inch 
test  section  indicate  that  the  tangential  grid  scheme  provides  a  useful  grid  which  is 
validated  by  the  experimental  data.  It  should  be  noted  that  even  though  both  of  the 
grid  schemes  tried  here  satisfied  the  minimum  desired  y+  value  of  1  only  the 
tangential  grid  scheme  provided  good  agreement  with  experimental  data.  In  this 
case  there  was  little  difference  among  grids  which  had  y+  values  less  than  5. 


x,  in. 


(a)  201  grid  points  in  y  direction. 

Figure  74.  Boundary-layer  thickness  as  a  function  of  longitudinal  tunnel  location 
computed  by  the  ins2d  for  24"  tunnel  at  x  =  16"  using  the  tangential  grid. 
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(b)  401  grid  points  in  y  direction. 
Figure  74.  Concluded. 


(a)  Tangential  scale  factor  of  32  and  201  grid  lines  vertically. 

Figure  75.  Boundary-layer  profile  computed  by  the  ins2d  for  24"  tunnel  at  x  =  16" 

using  the  tangential-transformation  grid. 
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(b)  Tangential  scale  factor  of  32  and  401  grid  lines  vertically 
Figure  75.  Continued. 


(c)  Tangential  scale  factor  of  64  and  201  grid  lines  vertically. 
Figure  75.  Concluded. 


b.  120  Inch  Test  Section 


The  grid  was  used  to  represent  the  120  inch  wind  tunnel  test  section 
(C8000)  from  the  floor  to  the  ceiling  and  between  the  longitudinal  stations  from  0 
to  120  inches.  The  grid  clustering  for  Roberts'  transformation  is  determined  by  the 
value  used  for  the  stretching  parameter  (3  (equation  36).  For  the  120"  test  section, 
the  boundary  layer  thickness  (8)  ranged  from  1.1  to  1.5  inches  depending  on  the 
longitudinal  survey  location.  This  leads  to  a  recommended  stretching  parameter  of 
about  1.05.  In  the  ins2d  code  the  resultant  grid  yields  a  minimum  y+  value  of  7.0, 
which  is  larger  than  the  recommended  value  of  about  1.  Using  201  grid  points  in 
the  y  direction,  normal  to  the  freestream,  three  values  of  (3  were  used  which 
resulted  in  a  reduction  in  the  minimum  y+  from  8.6  to  0.7.  None  of  these  cases 
resulted  in  converged  calculations  even  after  10,000  or  more  iterations.  As 
discussed  earlier,  this  is  an  undesirable  feature  of  these  grids.  The  following  table 
outlines  Roberts'  transformation  cases  which  were  used  to  evaluate  these  grids. 


C8000  Robert's  grid  (3 

min.  y+ 

iterations 

1.01 

8.6 

17,628 

1.001 

1.3 

10,126 

1.0005 

0.7 

10,000 

The  computed  boundary  layer  thicknesses  (8)  are  presented  in  Figure 
76  along  with  the  experimental  data  and  the  1/7  power  law  result.  The 
experimental  values  are  much  less  than  those  from  the  1/7  power  law,  as  noted 
earlier  in  the  discussion  of  the  experimental  data.  The  computational  result  from 
the  coarsest  grid  produced  thicknesses  which  were  less  than  one-half  of  that 
measured  in  the  experiment.  The  computational  result  from  the  two  finest  grids 
where  y+  is  near  1  provide  thicknesses  which  are  less  than  one-third  of  that 
measured  in  the  experiment.  These  inadequate  results  reflect  the  lack  of 
convergence  in  these  computations  and  suggest  that  the  Robert's  transformation 
grids  are  not  satisfactory  for  the  120  inch  test  section. 
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The  tangential-grid  scheme  was  also  used  for  the  120-inch  test- 
section  tunnel  (C8000).  Scale  factors  from  8  to  72  were  selected  to  achieve  a 
range  of  y+  values  which  approached  1.  The  number  of  y  grid  points  ranged  from 
201  to  601  to  evaluate  grid-density  effects.  In  contrast  to  the  tangential  grid 
transformation  for  the  24  inch  test  section,  not  all  of  these  cases  converged  in  less 
than  100  iterations.  The  y+  and  iterations  associated  with  the  tangential-grid 
scheme  are  presented  in  the  following  table: 


Tangential  grid 
scale  factor 

C8000  y  grid  points 
y  +  /  iterations 

201 

401 

601 

8 

8.3  /  8,600  nc 

6.2/ 10,115  nc 

16 

6.5  / 10,249  nc 

3.5/7,859  nc 

32 

3.1 /65  c 

1.7/ 9,334  nc 

64 

1.5 /48  c 

0.8/91  c 

0.6/1 1,271  nc 

72 

1.4 /46  c 

0.7  /  86  c 

... 

Here  nc  means  'not  converged  solution'  and  c  means  'converged  solution'.  It  is 
noted  that  while  all  of  the  converged  cases  required  less  than  100  iterations,  the 
unconverged  cases  were  run  for  about  8000  or  more  iterations.  It  is  noted  that 
with  grids  using  scale  factors  of  8  and  16,  the  solutions  did  not  converge 
Increasing  the  grids  from  401  to  601  y  grid  points  (using  a  scale  factor  of  64) 
caused  the  calculations  to  diverge  even  though  the  y+  value  was  reduced  from  0.8 
to  0.6. 


Calculated  boundary-layer  thicknesses  are  compared  in  Figure  77. 
Figure  77(a)  presents  the  results  for  201  y  grid  points.  Only  the  unconverged  scale 
factor  of  8  results  differ  significantly  from  the  experimental  data.  Even  the 
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unconverged  scale  factor  16  results  provide  good  agreement  with  the  experimental 
data.  As  a  result,  good  agreement  is  achieved  with  y+  of  6.5  down  to  1.4.  Figure 
77(b)  presents  the  results  for  401  y  grid  points.  None  of  the  unconverged  results 
agree  with  the  experimental  data.  The  scale  factor  8  result  overestimates  the 
boundary  layer  thickness  and  the  scale  factors  of  16  and  32  underestimate  the 
experimental  data.  Only  the  two  converged  results  for  scale  factors  of  64  and  72 
achieve  reasonable  agreement  with  the  experimental  data.  In  this  y  grid  size 
convergence  was  achieved  for  y+  of  1.7  or  less.  The  final  comparison  in  Figure 
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Figure  76.  Boundary-layer  thickness  computed  by  the  ins2d  for  120"  tunnel  using 
Roberts'  transformation  grid,  201  y  points. 
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(a)  201  grid  lines  vertically. 

Figure  77.  Boundary-layer  profile  computed  by  the  ins2d  for  120"  tunnel  using 

the  tangential-transformation  grid. 
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(b)  401  grid  lines  vertically. 
Figure  77.  Continued. 
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(c)  Tangential  scale  factor  of  64. 

Figure  77.  Concluded. 

77(c)  for  scale  factor  64  grids  show  that  increasing  the  number  of  y  grid  points 
from  201  to  401  to  601  reduces  the  y+  from  1.5  to  0.8  to  0.6  for  the  finest  grid. 
Yet  the  worst  result  was  from  the  finest  grid,  which  is  also  the  only  unconverged 
result  shown  in  this  plot. 

In  summary,  the  computational  results  from  both  of  the  test  section 
lengths  used  in  the  experimental  investigation,  show  that,  first,  a  grid  which 
provides  a  converged  solution  is  needed  before  a  reasonable  comparison  with 
experimental  data  can  be  expected.  The  value  of  y+  is  of  less  importance  as 
shown  by  the  Robert's  grid  where  a  y+  of  0.7  gave  poor  results;  yet  a  tangential 
grid  with  a  y+  of  7  gave  good  results.  The  appropriate  grid  clustering,  size,  and 
density  are  not  always  apparent  before  comparing  results  with  experimental  data. 
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2.  Driven  Cavity  Computations 


A  frequently  used  model  problem  for  cavity  flow  is  the  driven  cavity.  This 
cavity  has  three  stationary  walls  and  one  wall  which  moves  at  a  constant  velocity. 
The  velocity  magnitude  is  determined  by  the  Reynolds  number  of  the  calculation. 
Ghia  et  al  [Ref.  67]  conducted  a  very  detailed  analysis  using  a  vorticity-stream- 
function  formulation  for  the  incompressible  Navier-Stokes  equations.  Uniform 
grids  were  used  in  these  calculations.  For  Reynolds  numbers  from  100  to  3200  the 
grid  was  129  by  129;  and  for  Reynolds  numbers  of  5000  or  greater  the  grid  was 
257  by  257.  An  example  of  their  streamline  flow  pattern  was  presented  in  Figure 
10  for  a  Reynolds  number  of  10,000.  In  this  paper,  flows  were  computed  for 
Reynolds  numbers  from  1000  to  10,000  and  then  the  results  were  compared  with 
those  from  several  other  investigators  [Ref.  68,  95,and  96].  The  Ghia  u  velocity 
results  along  a  vertical  line  which  passes  through  the  center  of  the  cavity  compared 
with  the  other  computations  in  Figure  78.  The  origins  in  Figure  78  are  displaced 
to  separate  the  various  curves.  At  low  Reynolds  numbers  the  boundary  layers  are 
very  thick.  As  Reynolds  number  reaches  and  exceeds  5000,  8  asymptotes  toward 
a  converged  magnitude.  Away  from  the  cavity  walls,  the  velocities  tend  toward  a 
linear  variation  especially  for  the  high  Reynolds  number  cases  where  the  boundary 
layer  thickness  is  a  small  fraction  of  the  cavity  depth.  In  most  of  these 
comparisons,  the  Ghia  results  are  very  consistent  with  those  of  the  other  three 
investigators.  The  most  notable  difference  occurs  when  the  Ghia  results  are 
compared  with  those  from  Ref.  95,  which  used  the  coarsest  grid  (50  x  50)  of  the 
four  sets  of  results  shown. 

Comparable  calculations  done  using  the  ins2d  incompressible  Navier- 
Stokes  code  [Ref.  91  and  92]  for  Reynolds  numbers  of  1000  and  10,000  are 
presented  in  Figure  79  with  the  results  from  Ghia  et  al  [Ref.  67].  In  Reference  91 
it  was  found  that  a  clustered  81  by  81  grid  agreed  with  the  Ghia  et  al  [Ref.  67] 
results.  For  the  present  investigation,  an  ins2d  calculation  was  made  using  an  81 
by  8 1  grid  with  tangential  spacing.  The  grid  was  generated  using  the  cavity  grid 


177 


-0.4  -0.2  0.0 


u 

g.c. 

0-0  0.2  0.4 


0.6  0.8  1.0 


y 

1.0 

0.8 

0.6 

0.4 

0.2 

0.0 


Figure  78.  Comparison  of  u  velocity  along  vertical  lines  through  the  geometric 

center  of  the  cavity  [Ref.  67]. 
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generation  code  cgg.f  as  described  previously.  The  u  (Fig.  79(a))  and  v  (Fig. 
79(b))  velocity  components  on  vertical  or  horizontal  lines,  respectively,  are  non- 
dimensionalized  by  the  moving  wall  velocity.  These  comparisons  show  very  good 
agreement  for  both  velocity  components.  This  calculation  is  consistent  with 
previous  ins2d  results  and  confirms  that  consistent  grids  and  boundary  conditions 
are  being  used  in  the  present  investigation. 


(a)  u  velocities. 

Figure  79.  Comparison  between  ins2d[Ref.  91  and  92]  and  Ghia  et  al  [Ref.  67] 
Navier-Stokes  computational  results  for  Re  of  1000  and  10,000. 
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(b)  v  velocities. 
Figure  79.  Concluded. 


3.  Open  Cavity  Flow  Computations 


The  ins2d  results  of  the  previous  two  sections  provided  consistent  results 
for  the  two  SJSU  tunnel  test  section  configuration's  (24"  and  120")  boundary 
layers  and  for  a  driven  cavity  problem.  These  data  provided  the  confidence 
needed  to  progress  to  the  computations  of  the  wind  tunnel  with  a  cavity  and  with  a 
test  section  diffuser. 

Initial  computational  efforts  focused  on  the  24"  test  section  with  one  cavity 
(D  =  1").  The  tunnel  grid  from  the  tunnel  grid  generation  program  stgg.f  was 
used  as  a  starting  point.  The  single  cavity  was  added  using  the  cavity  grid 
generation  program  cgga.f  .  This  program  used  an  81  by  81  cavity  grid  with 
tangential  spacing,  using  a  scale  factor  of  16.  Then  the  tunnel  grid  with  a  scale 
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factor  of  32  vertically  was  augmented  from  201  high  by  495  long  to  satisfy  the 
longitudinal  spacing  criteria  after  the  cavity  grid  was  introduced.  This  criteria  uses 
a  maximum  grid  spacing  increase  of  no  more  than  20  percent  between  adjacent 
grid  cells.  The  augmentation  increased  the  grid  to  201  high  by  495  long  for  the 
basic  24"  test  section.  The  resultant  experimental  and  computational  boundary 
layer  thicknesses  are  presented  in  Figure  80  for  both  the  tunnel  without  a  cavity 
and  the  tunnel  with  the  1"  cavity.  These  results  show  a  small  increase  (~  0.03")  in 
5  when  the  cavity  was  present.  In  contrast  the  experimental  data  show  a  crossover 
in  8  ,  so  that  at  x  »  5  the  basic  tunnel  has  a  slightly  thinner  8,  while  at  x  =  16  the 
basic  tunnel  has  an  increased  8.  In  the  region  immediately  downstream  of  the 
cavity  there  is  only  a  small  difference  in  the  8  for  these  two  configurations.  The 
boundary  layer  profiles  were  very  similar  to  the  one  shown  in  Figure  75(a). 

Similar  computations  were  done  for  the  120"  test  section  both  with  and 
without  the  7°  test  section  diffuser.  In  this  case  the  same  81  by  81  tangentially- 
spaced  cavity  grid  was  used.  The  freestream  grid  was  increased  from  321  to  561 
by  cgg.f  to  satisfy  grid  spacing  change  rate  requirements.  The  ins2d  code  did  not 
converge  to  a  solution  for  this  case.  The  freestream  grid  was  split  into  two  grids 
(201  by  321  each)  and  merged  using  the  PEGSUS  overset  scheme  [Ref.  94].  The 
first  grid  extended  from  x  =  0"  to  96"  and  the  second  grid  from  x  =  96"  to  120". 
This  provided  a  denser  grid  in  the  test  section  region  where  the  boundary  layer 
thicknesses  were  measured  and  where  the  cavity  and  test  section  diffuser  were 
located.  For  the  cavity  cases,  the  cavity  grid  was  overset  as  a  third  grid  using 
PEGSUS.  These  computations  converged  in  about  200  to  240  iterations 
depending  on  the  specific  configuration. 

The  experimental  and  computational  boundary  layer  thickness  results  are 
presented  for  the  test  section  region  in  Figure  81.  The  experimental  boundary- 
layer  thickness  was  slightly  greater  with  the  1 "  cavity  than  without  the  cavity.  The 
computed  8's  are  nearly  the  same  as  the  cavity  case  where  the  experimental  data 
are  a  small  increment  (=  0.09")  thicker.  For  the  case  with  the  test  section  diffuser 
(Figure  82)  the  experimental  data  for  the  cavity  case  has  a  slightly  thinner  8  (A8 
~  0.03")  while  the  computational  data  has  a  little  larger  difference  (A8  ~  0.05"). 
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While  the  computational  results  shown  in  Figure  80  through  82  have  some 
disagreement  with  a  few  experimental  data  points,  the  trends  tend  to  be  consistent 
with  the  experimental  data  trends  for  the  zero  pressure  gradient  in  both  test 
sections;  the  presence  of  the  cavity  slightly  increases  the  boundary  layer  thickness. 
For  the  120"  test-section  case  with  the  7°  test-section  diffuser,  the  cavity  slightly 
decreases  the  boundary-layer  thickness  8.  Overall,  both  experimental  and 
computational  results  show  only  small  changes  in  the  boundary  layer 
characteristics  due  to  the  presence  of  a  1"  cavity. 


Figure  80.  Boundary-layer  thickness  as  a  function  of  longitudinal  tunnel  location 

computed  by  ins2d  for  24"  tunnel. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


One  airfoil  high-lift  concept  uses  multiple  cavities  on  its  upper  surface  to 
increase  camber.  The  usefulness  of  cavities  in  an  adverse  pressure  gradient  is 
determined  by  their  effect  on  the  viscous  flow  in  the  boundary  layer  downstream 
of  the  cavities.  Maintenance  of  attached  flow  requires  that  the  flow  over  the 
cavities  and  downstream  of  them  retains  enough  momentum  to  overcome  the 
kinetic  energy  loss  due  to  an  adverse  pressure  gradient,  shear-layer  flow  gradients, 
and  viscous  dissipation.  For  this  concept's  feasibility  to  be  demonstrated  the  effect 
of  cavities  on  boundary-layer  characteristics  especially  in  an  adverse  pressure 
gradient  must  be  determined.  The  present  investigation  concentrated  on  the  effect 
of  cavity  flow  on  the  attached  surface  boundary-layer  characteristics  in  adverse 
pressure  gradients. 

The  purpose  of  the  present  investigation  was  to  experimentally  and 
computationally  determine  the  effect  of  1,  2,  or  4  two-dimensional  span  wise 
cavities  with  square  cross-sections  on  the  boundary-layer  characteristics  in  both  a 
constant-pressure  flow  and  for  two  adverse  pressure-gradient  flows.  The 
experimental  investigation  was  conducted  in  the  San  Jose  State  University  (SJSU) 
12"  by  12"  low-speed  tunnel.  Two  test-section  lengths,  24"  and  120",  were  used  to 
obtain  a  boundary-layer  thickness  relative  to  the  cavity  depth  5/D  which  was 
either  less  than  0.5  ( 5/D  <  0.5  ),  where  unsteady  flow  oscillations  may  exist  in  the 
cavity  (or  cavities),  or  greater  than  1.0  ( 5/D  >  1.0  )  where  there  is  predominately  a 
steady,  standing  vortex  flow  in  the  cavity  (or  cavities).  Adverse  pressure  gradients 
(dp/dx  >  0)  were  obtained  using  7°  and  14°  inserts  on  the  floor  of  the  test  section. 
Total  pressure  profiles  in  the  boundary  layer  were  measured  at  several  longitudinal 
locations  both  upstream  and  downstream  of  the  cavities.  Computational  results 
were  obtained  using  a  numerical  solution  to  the  incompressible  Navier-Stokes 
equations.  Additional  tests  were  conducted  with  a  single  cavity  which  was  either 
open  or  closed  to  determine  the  effects  of  surface  discontinuities  on  the  boundary- 
layer  development. 

The  measured  profiles  demonstrated  boundary  layer  similarity  for  the  basic 
test  section  and  for  the  7°  test-section  diffuser.  There  was  a  significant  boundary- 
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layer  profile  variation  from  the  1/7  power  law  for  the  14°  test-section  diffuser 
which  was  typical  of  an  adverse  pressure  gradient.  The  single  and  double  cavities 
showed  only  a  small  deviation  from  the  profile  measured  with  no  cavity.  The 
four-cavity  configuration  profile,  in  most  cases,  represented  a  boundary  layer 
profile  similar  to  those  found  in  an  adverse  pressure  gradient.  All  of  the  surface 
pressure  data  suggest  that  the  presence  of  the  cavities  had  no  noticeable  effect.  It 
may  be  concluded  that  cavities  may  be  deployed  with  only  a  small  change  to  the 
boundary-layer  profile  and  without  significantly  modifying  the  resultant  pressure 
distribution.  This  important  conclusion  shows  that  this  high-lift  concept  may  be 
feasible  because  the  multiple  cavities  make  only  small  changes  in  the  boundary- 
layer  characteristics. 

Another  aspect  of  this  application  is  the  effect  of  surface  non-uniformities 
which  could  be  caused  by  cavity  doors  on  the  boundary  layer  and  on  the  surface 
pressure  distribution.  As  expected,  a  poor  fit  of  a  cavity  door  could  increase 
boundary-layer  thickness  and  the  skin-friction  drag.  There  were  differences  in  the 
longitudinal  pressure  distribution  for  both  test  sections.  These  results  indicated 
that  the  surface  uniformity  for  a  closed  cavity  had  an  effect  on  the  freestream 
characteristics,  which  determine  the  surface  pressure  distribution.  It  was  also 
found  that  the  open,  cavities  produced  less  change  in  the  surface  pressure 
distribution  than  did  surface  non-uniformity  in  an  otherwise  plain  ceiling.  This  is 
a  further  indication  of  the  feasibility  of  the  use  of  cavities. 

In  summary,  the  computational  results  for  both  of  the  test-section  lengths 
used  in  the  experimental  investigation,  show  that  a  grid  which  provides  a 
converged  solution  is  needed  before  a  reasonable  comparison  with  experimental 
data  can  be  expected.  A  small  value  of  y+  is  of  less  importance  than  convergence. 
This  was  shown  by  the  unconverged  Roberts’  transformation  grid  with  a  y+  of  0.9 
which  gave  poor  results,  while  a  converged  tangential  grid  with  a  y+  of  7  gave 
good  results.  Computations  for  the  120"  test  section  with  either  a  one  inch  cavity 
and/or  with  a  7°  test-section  diffuser  showed  only  small  changes  in  the  boundary- 
layer  characteristics.  The  appropriate  grid  clustering,  size,  and  density  were  not 
always  apparent  before  comparing  computational  results  with  experimental  data. 
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Stated  another  way,  it  was  found  that  the  experimental  data  were  needed  to  guide 
the  computational  effort. 

This  experimental  and  computational  investigation  has  shown  that  multiple 
open  cavities  in  either  zero  or  adverse  pressure  gradients  make  only  small  changes 
in  the  boundary-layer  characteristics.  It  was  also  found  that  a  small  backward¬ 
facing  step  made  larger,  adverse  changes  to  the  boundary  layer  than  those  made  by 
the  cavities.  It  is  concluded  that  multiple  cavities  may  be  a  feasible  high-lift 
concept.  It  is  recommended  that  the  cavity  concept  be  implemented  in  an  airfoil 
and  that  it  be  tested  in  a  wind  tunnel  to  quantify  its  high-lift  and  drag 
characteristics. 
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APPENDIX  A 


TEST  RUN  SCHEDULE 


The  run  schedule  for  the  experimental  investigation  is  presented  in  this 
appendix.  It  lists  the  following  information: 

Run  Sequential  number  assigned  to  a  group  of  data 


Configuration  A  four  digit  code  which  describes  the  tunnel  configuration 
code  The  first  digit  identifies  the  test  section  length: 

0  or  1  refers  to  the  24"  long  test  section 
8  refers  to  the  120"  long  test  section 
The  second  digit  identifies  the  diffuser  angle: 

0  means  basic  test  section  with  no  diffuser 
7  means  the  7°  test  section  diffuser 
14  means  the  14°  test  section  diffuser 
The  third  digit  identifies  the  cavity  depth  in  inches: 

0  means  no  cavity 

1  means  1  inch  deep  cavity(-ies) 

2  means  2  inch  deep  cavity 

3  means  2  inch  deep  cavity  with  a  movable  door  which 

can  be  used  to  close  the  cavity  or  open  the  cavity 
to  the  freestream  flow  (i.e.,  c-closed  or  o-open). 
The  fourth  digit  identifies  the  number  of  cavities  in  the 
streamwise  direction  as  either  0, 1, 2,  or  4  cavities. 
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Run  type 

boundary  layer 
probe  position 
b.l.  probe  height 

b.l.  rake  installed 

Notes 


This  identifies  the  primary  data  taken  as  either  CP,  surface 
Cp  distribution,  or  BL,  boundary  layer  profiles. 

This  identifies  the  x  location  of  the  boundary  layer  probe  tip. 

This  identifies  the  vertical  dimension  of  the  boundary  layer 
probe. 

This  identifies  whether  the  b.l.  rake  is  installed  at  the  x=16" 
location. 

Notes  related  to  the  run  describing  unexpected  items  or 
information  unique  to  a  particular  run  which  is  not  described 
in  the  other  columns. 


An  entry  is  made  in  the  run  schedule  only  when  there  is  a  change  in  a 
parameter.  If  there  is  no  entry  then  the  parameter  is  unchanged  from  the  previous 
run. 
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202 


Configuration  Run  boundary  b.l.  probe  b.l.  rake  Notes 

code  Type  layer  probe  head  installed 


203 


205 


255  _ _ 0.034 _ New  Pt,bl  probe  w/  0.034"  pro 

256  _ 0.022 _ Tried  a  modified  0.022"  probe 

257  Tried  a  second  0.022”  probe 


Configuration  Run  boundary  b.l.  probe  b.l.  rake  Notes 

code  Type  layer  probe  head  installed 
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APPENDIX  B 


SURFACE  PRESSURE  TAP  LOCATIONS 


The  ceiling  surface  and  cavity  pressure  tap  locations  used  to  obtain  pressure 
distributions  in  the  present  investigation  are  listed  in  this  appendix.  It  includes  the 
following  tables: 


Table 

number 

of 

ceiling 

taps 

number 
of  cavity 
taps 

Surface 

1 

16 

N/A 

7°  test  section  diffuser 

2 

16 

N/A 

14°  test  section  diffuser 

3 

19 

N/A 

basic  ceiling  (i.e.,  no  cavities) 

4 

20 

15 

one  2"  deep  cavity  and  adjacent  surfaces 

5 

16 

8 

one  2"  deep  cavity  w/  trap  door  and  adjacent  surfaces 

6 

17 

11 

one  1"  deep  cavity  and  adjacent  surfaces 

7 

17 

21 

two  1"  deep  cavities  and  adjacent  surfaces 

8 

16 

44 

four  1"  deep  cavities  and  adjacent  surfaces 

The  basic  test  section  of  this  wind  tunnel  uses  only  the  basic  ceiling  whose  pressure 
tap  locations  are  given  in  Table  3  and  a  plain  uninstrumented  tunnel  floor.  When 
an  adverse  pressure  gradient  is  desired  one  of  the  test  section  diffusers  whose 
pressure  tap  locations  are  given  in  either  Table  1  or  Table  2  is  added  on  the  tunnel 
floor  along  with  a  modified  contraction  at  the  test  section  entrance.  When  a  cavity 
configuration  is  desired  the  basic  ceiling  is  replaced  by  the  appropriate  ceiling 
given  in  Tables  4  through  8.  All  of  the  x  station  dimensions  in  these  Tables  are 
given  for  the  24"  long  test  section.  The  120"  long  test  section  is  achieved  by 
adding  a  96"  insert  with  a  variable  floor  location  upstream  of  the  original  test 
section.  As  a  consequence  the  x  stations  for  the  120"  long  test  section  may  be 
obtained  by  adding  96"  to  the  x  station  values  given  in  these  Tables. 
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Table  1.  7°  test  section  diffuser  pressure  tap  locations: 


Tap 

x,  inches 

y,  inches 

1 

2.031 

-4 

2 

3.015 

-4 

3 

4.031 

-4 

4 

5.063 

-4 

5 

7.148 

-4.128 

6 

7.985 

-4.221 

7 

8.978 

-4.294 

8 

9.97 

-4.441  1 

9 

10.947 

-4.550 

10 

11.924 

-4.658 

11 

12.917 

-4.769 

12 

13.909 

-4.879 

13 

15.925 

-5.103  1 

14 

17.942 

-5.327 

15 

19.927 

-5.547 

16 

21.881 

-5.765 

Table  2.  14°  test  section  diffuser  pressure  tap  locations: 


Tap 

x,  inches 

y,  inches 

1 

2.000 

-1.500 

2 

3.000 

-1.500  1 

3 

4.015 

-1.500 

1  4 

5.000 

-1.500 

5 

7.395 

-1.849 

6 

7.910 

-1.978 

7  1 

8.881 

-2.220 

8 

9.821 

-2.455 

9 

10.791 

-2.958 

10 

11.791 

-2.948 

11 

12.762 

-3.191 

12 

13.732 

-3.433 

13 

15.703 

-3.926 

14 

r  F7i>44  * 

-4.411 

15 

20.524 

-5.131 

16 

21.524 

-5.381 
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Table  3.  Basic  test  section  ceiling. 


x,  inches 

y, inches 

2.094 

6 

3.01 

6 

4.063 

6 

5.083 

6 

6.063 

6 

7.073 

6 

8.063 

6 

9.063 

6 

10.073 

6 

11.083 

6 

12.083 

6 

13.073 

6 

14.083 

6 

15.063 

6 

16.063 

6 

17.063 

6 

18.073 

6 

20.063 

6 

22.063 

6 

Table  4.  One  2"  deep  cavity  and  the  adjacent  ceiling  surface. 


x,  inches 

2 

6 

3.015 

6 

4 

6 

5.015 

6 

6.041 

6 

6.52 

6 

6.97 

6 

9.5313 

6 

9.98 

6 

10.48 

6 

11.015 

6 

12.01 

6 

13.01 

6 

14.015 

6 

15.01 

6" 

16.02 

6 

17.02 

6 

18.02 

6 

IH9 

6 

22.03 

6 
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Table  5.  One  2"  deep  cavity  with  a  trapdoor  and  the  adjacent  ceiling  surface. 


x,  inches 

2.03 

6 

2.97 

6 

4 

6 

4.98 

6 

8.06 

6 

10 

6 

11.03 

6 

12.02 

6 

13.08 

6 

14.063 

6 

15.07 

6 

16.063 

6 

17.06 

6 

18.06 

6 

20.063 

6 

22.063 

6 

Table  6.  One  1"  deep  cavity  and  the  adjacent  ceiling  surface. 


x,  inches 


2.031 


3 


6.24 


8.531 


10.016 


11.031 


12 


13.016 


14.016 


15.031 


16.062 


17.031 


18.047 


20.062 


22.031 
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Table  7.  Two  1"  deep  cavities  and  the  adjacent  ceiling  surface. 


x, inches 

y,  inches 

1.969 

6 

2.969 

6 

3.938 

6 

4.969 

6 

5.75 

6 

7.75 

6 

10 

6 

10.969 

6 

11.969 

6 

12.984 

6 

13.984 

6 

14.984 

6 

16 

6 

16.969 

6 

17.984 

6 

19.969 

6 

21.969 

6 

x,  inches  y, inches 


Cavity 

number 


1 


Table  8.  Four  1"  deep  cavities  and  the  adjacent  ceiling  surface. 


■M 

x,  inches 

y,  inches 

ns 

x,  inches 

y,  inches 

i 

1.969 

6  : 

i 

6 

6.48 

2 

2.938 

6 

2 

6 

6.87 

3 

3.95 

6 

3 

6.1 

6.95 

4 

4.938 

6 

4 

6.25 

6.94 

5 

5.72 

6 

5 

6.5 

6.935 

6 

7.67 

6 

6 

6.75 

6.93 

7 

9.71 

6 

7 

6.91 

6.92 

8 

11.69 

6 

8 

6.94 

6.89 

9 

13.2 

6 

9 

6.94 

6.75 

10 

13.94 

6 

io 

6.94 

6.5 

11 

14.94 

6 

11 

6.94 

6.25 

12 

15.97 

6 

12 

7.95 

6.46 

13 

16.97 

6 

13 

7.95 

6.86 

14 

17.95 

6 

14 

8.05 

6.93 

15 

19.95 

6 

15 

8.2 

6.93 

16 

21.94 

6 

16 

8.47 

6.935 
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APPENDIX  C 


SURFACE  PRESSURES 


The  pressure  distributions  measured  on  the  wind  tunnel  surfaces  in  the 
cavities  and  on  the  adjacent  freestream  surfaces  are  presented  in  this  appendix. 
The  data  from  each  cavity  are  presented  in  three  contiguous  plots:  (1)  front  cavity 
wall;  (2)  cavity  floor;  and  (3)  rear  cavity  wall.  The  freestream  surface  plots  are 
presented  for  all  of  the  pressures  measured  on  the  y  =  6  inches  surface.  The  data 
are  faired  only  for  the  surfaces  upstream  and  downsteam  of  the  cavity  (or  cavities). 
The  following  table  relates  the  configurations  and  their  descriptions  with  the  figure 
numbers: 


Config¬ 

uration 

Test 

Section 

Length 

Number  of 

Cavities 

Cavity 

Depth 

Surface 

Pressures 

Cavity 

Pressures 

0011 

24  in. 

1 

1  in. 

Cl 

C2 

0711 

24  in. 

1 

1  in. 

C3 

C4 

1411 

24  in. 

1 

1  in. 

C5 

C6 

8011 

120  in. 

1 

lin. 

C7 

C8 

8711 

120  in. 

1 

1  in. 

C9 

CIO 

8411 

120  in. 

1 

1  in. 

Cll 

C12 

0012 

24  in. 

2 

1  in. 

C13 

C14 

0014 

24  in. 

4 

1  in. 

C15 

C16 
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Figure  Cl.  Pressure  distribution  on  the  tunnel  surfaces  adjacent  to  the  0011 


Figure  C2.  Pressure  distribution  in  the  0011  cavity 


Figure  C3.  Pressure  distribution  on  the  tunnel  surfaces  adjacent  to  the  0711 


Figure  C4.  Pressure  distribution  in  the  071 1  cavity 


CP 


Figure  C5.  Pressure  distribution  on  the  tunnel  surfaces  adjacent  to  the  141 1 

cavity. 

Cavity 
0.3 

0.2 

0.1 

0 

-0.1 

Y:  6  7  (6.938  in.)  7  6 

X:  (6.50in.)  6.5  7.5  (7.50  in.) 

Figure  C6.  Pressure  distribution  in  the  141 1  cavity. 
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Figure  Cl.  Pressure  distribution  on  the  tunnel  surfaces  adjacent  to  the  8011 


Figure  C8.  Pressure  distribution  in  the  8011  cavity 


Figure  C9.  Pressure  distribution  on  the  tunnel  surfaces  adjacent  to  the  871 1 


Figure  CIO.  Pressure  distribution  in  the  871 1  cavity 
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Figure  Cl 2.  Pressure  distribution  in  the  8411  cavity. 


Figure  Cl 3.  Pressure  distribution  on  the  tunnel  surfaces  adjacent  to  the  0012 

cavities. 
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Figure  C14.  Pressure  distribution  in  the  0012  cavities 


Figure  Cl 5.  Pressure  distribution  on  the  tunnel  surfaces  adjacent  to  the  0014 

cavities. 
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(a)  Pressure  distribution  in  the  first  two  cavities. 
Figure  Cl 6.  Pressure  distribution  in  the  0014  cavities 


(b)  Pressure  distribution  in  the  second  two  cavities 
Figure  Cl 6.  Concluded. 


Figure  Cl 9.  Pressure  distribution  on  the  tunnel  surfaces  adjacent  to  the  1412 

cavities. 
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Figure  C21.  Pressure  distribution  on  the  tunnel  surfaces  adjacent  to  the  1414 

cavities. 


(a)  Pressure  distribution  in  the  first  two  cavities. 
Figure  C22.  Pressure  distribution  in  the  1414  cavities. 
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(b)  Pressure  distribution  in  the  second  two  cavities. 
Figure  C22.  Concluded. 


x,  in. 

Figure  C23.  Pressure  distribution  on  the  tunnel  surfaces  adjacent  to  the  8012 

cavities. 


231 


Figure  C24.  Pressure  distribution  in  the  8012  cavities 
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(a)  Pressure  distribution  in  the  first  two  cavities. 
Figure  C26.  Pressure  distribution  in  the  8014  cavities. 
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(b)  Pressure  distribution  in  the  second  two  cavities. 
Figure  C26.  Concluded. 
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Figure  C27.  Pressure  distribution  on  the  tunnel  surfaces  adjacent  to  the  8714 

cavities. 
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(a)  Pressure  distribution  in  the  first  two  cavities. 
Figure  C28.  Pressure  distribution  in  the  8714  cavities. 
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Figure  C30.  Pressure  distribution  in  the  8412  cavities 
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(a)  Pressure  distribution  in  the  first  two  cavities. 
Figure  C32.  Pressure  distribution  in  the  8414  cavities. 
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(b)  Pressure  distribution  in  the  second  two  cavities. 
Figure  C32.  Concluded. 
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Figure  C33.  Pressure  distribution  on  the  tunnel  surfaces  adjacent  to  the  0021 


Figure  C34.  Pressure  distribution  in  the  0021  cavity 


Figure  C35.  Pressure  distribution  on  the  tunnel  surfaces  adjacent  to  the  8021 


Figure  C36.  Pressure  distribution  in  the  8021  cavity 


Figure  C37.  Pressure  distribution  on  the  tunnel  surfaces  adjacent  to  the  8421 

cavity. 
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Figure  C38.  Pressure  distribution  in  the  8421  cavity. 
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Figure  C39.  Pressure  distribution  on  the  tunnel  surfaces  adjacent  to  the  1421 

cavity. 
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Figure  C40.  Pressure  distribution  in  the  1421  cavity. 


241 


c 


p 


0.3 

0.2 

0.1 

0 

-0.1 

0  5  10  15  20  25 

x,  in. 


:  — O —  Upstream,  < 
r— 0—  Downstreai 
:  -  -  Upstream,  < 

3pen 

n,  open  j 
closed 

:  —  I 

)ownstreai 

n,  closed  ! 

. 

i 

■  ■  ■  ■  *  ■  ■  ■  ■  *  ■  ■  ■  ■  1 

Figure  C41.  Pressure  distribution  on  the  tunnel  surfaces  adjacent  to  the  0031 

cavity. 
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Figure  C42.  Pressure  distribution  in  the  0031  cavity. 
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Figure  C43.  Pressure  distribution  on  the  tunnel  surfaces  adjacent  to  the  0731 


Figure  C44.  Pressure  distribution  in  the  00731  cavity 


Figure  C45.  Pressure  distribution  on  the  tunnel  surfaces  adjacent  to  the  1431 

cavity. 
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Figure  C46.  Pressure  distribution  in  the  1431  cavity. 
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Figure  C47.  Pressure  distribution  on  the  tunnel  surfaces  adjacent  to  the  8031 


Figure  C48.  Pressure  distribution  in  the  8031  cavity 
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Figure  C49.  Pressure  distribution  on  the  tunnel  surfaces  adjacent  to  the  8431 

cavity. 
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Figure  C50.  Pressure  distribution  in  the  8431  cavity. 
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APPENDIX  D 


BOUNDARY  LAYER  PROFILES 


The  boundary  layer  data  from  this  investigation  is  presented  for  each  wind 
tunnel  configuration  in  this  appendix.  There  were  four  to  six  longitudinal 
locations  for  each  configuration  where  boundary  layer  profiles  were  measured. 
Each  figure  identifies  the  x  location  for  the  data  presented.  The  displacement 
thickness  was  computed  from  the  data  and  then  used  to  non-dimensionlize  the 
distance  y  which  was  measured  from  the  wind  tunnel  ceiling.  A  pitot  probe 
measured  the  total-pressure  variation  and  the  static  pressure  was  obtained  from  a 
port  located  on  the  ceiling  at  the  appropriate  longitudinal  location.  As  expected, 
the  static  pressures  were  nearly  constant  and  the  total  pressures  varied  from  a  few 
pounds  per  square  foot  greater  than  the  static  pressure  near  the  surface  to  the 
freestream  total  pressure  at  the  boundary  layer  edge.  The  following  table  relates 
the  configurations  and  their  description  with  the  figure  numbers: 


Config¬ 

uration 

Test 

Section 

Length 

Number  of 

Cavities 

Cavity 

Depth 

Figure 

Number 

Notes 

0000 

24  in. 

0 

N/A 

D1 

8000 

120  in. 

0 

N/A 

D2D 

0700 

24  in. 

0 

N/A 

D3 

8700 

120  in. 

0 

N/A 

D4D 

1400 

24  in. 

0 

N/A 

D5 

8400 

120  in. 

0 

N/A 

D6 

Config¬ 

uration 

T.S.  Length 

No.  of 

Cavities 

1 

Figure  No. 

Notes 

0011 

24  in. 

1 

N/A 

D7 

0021 

24  in. 

1 

2  in. 

D8 

0711 

24  in. 

1 

1  in. 

D9 

1411 

24  in. 

1 

lin. 

D10 

1421 

24  in. 

1 

2  in. 

Dll 

8011 

120  in. 

1 

1  in. 

D12 

8021 

120  in. 

1 

2  in. 

D13 

8711 

120  in. 

1 

1  in. 

D14 

8411 

120  in. 

1 

1  in. 

D15 

8421 

120  in. 

1 

2  in. 

D16 

0012 

24  in. 

2 

1  in. 

D17 

0014 

24  in. 

4 

lin. 

D18 

0714 

24  in. 

4 

lin. 

D19 

1412 

24  in. 

2 

1  in. 

D20 

1414 

24  in. 

4 

1  in. 

D21 

8012 

120  in. 

2 

1  in. 

D22 

8014 

120  in. 

4 

1  in. 

D23 

8714 

120  in. 

4 

lin. 

D24 
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Config¬ 

uration 

T.S.  Length 

No.  of 

Cavities 

Figure  No. 

Notes 

8412 

120  in. 

2 

1  in. 

D25 

8414 

120  in. 

4 

1  in. 

D26 

0031c 

24  in. 

1 

2  in. 

D27 

003  lo 

24  in. 

1 

2  in. 

D28 

0731c 

24  in. 

1 

2  in. 

D29 

073  lo 

24  in. 

1 

2  in. 

D30 

1431c 

24  in. 

1 

2  in. 

D31 

1431o 

24  in. 

1 

2  in. 

D32 

8031c 

120  in. 

1 

2  in. 

D33 

803  lo 

120  in. 

1 

2  in. 

D34 

8431c 

120  in. 

1 

2  in. 

D35 

2  in. 


8431o 


120  in. 


1 


D36 


U/Ue 


Figure  Dl.  Boundary  layer  profiles  for  0000  configuration. 
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Figure  D2.  Boundary  layer  profiles  for  8000  configuration. 
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Figure  D3.  Boundary  layer  profiles  for  0700  configuration. 


Figure  D4.  Boundary  layer  profiles  for  8700  configuration. 


Figure  D6.  Boundary  layer  profiles  for  8400  configuration. 


Figure  D9.  Boundary  layer  profiles  for  0711  configuration. 


Figure  DIO.  Boundary  layer  profiles  for  141 1  configuration. 


Figure  Dll.  Boundary  layer  profiles  for  1421  configuration. 


Figure  D12.  Boundary  layer  profiles  for  8011  configuration. 


Figure  D13.  Boundary  layer  profiles  for  8021  configuration. 


Figure  D14.  Boundary  layer  profiles  for  8711  configuration. 


Figure  D15.  Boundary  layer  profiles  for  8411  configuration. 


Figure  D16.  Boundary  layer  profiles  for  8421  configuration. 


Figure  D19.  Boundary  layer  profiles  for  0714  configuration. 


Figure  D20.  Boundary  layer  profiles  for  1412  configuration. 


Figure  D21.  Boundary  layer  profiles  for  1414  configuration. 


Figure  D22.  Boundary  layer  profiles  for  8012  configuration. 


Figure  D23.  Boundary  layer  profiles  for  8014  configuration. 


Figure  D24.  Boundary  layer  profiles  for  8714  configuration. 


Figure  D25.  Boundary  layer  profiles  for  8412  configuration. 


Figure  D26.  Boundary  layer  profiles  for  8414  configuration. 


Figure  D27.  Boundary  layer  profiles  for  0031c  configuration. 


Figure  D28.  Boundary  layer  profiles  for  003  lo  configuration. 
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Figure  D29.  Boundary  layer  profiles  for  0731c  configuration. 
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Figure  D30.  Boundary  layer  profiles  for  073  lo  configuration. 
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Figure  D31.  Boundary  layer  profiles  for  1431c  configuration. 


Figure  D32.  Boundary  layer  profiles  for  1431  ©configuration. 


Figure  D33.  Boundary  layer  profiles  for  8031c  configuration. 


Figure  D34.  Boundary  layer  profiles  for  803  lo  configuration. 
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